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provided  by  Mr.  Robert  E.  McCarty  (AFWAL/FIER)  as  the  Air  Force 
Project  Engineer.  The  work  described  herein  represents  a  continu¬ 
ation  of  previous  developments  performed  in-house  at  the  University 
of  Dayton  Research  Institute,  and  on  Air  Force  Contract 
F33615-76-C-3103. 
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Mr.  Robert  E.  McCarty.  The  analytical  development  performed  by 
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SECTION  1 


INTRODUCTION 


Finite  element  methods  and  other  numerical  techniques  for 
performing  complex  structural  analyses  have  matured  to  the  point 
that  they  may  be  used  with  confidence  in  the  development  and 
final  qualification  of  complex  structural  designs.  Numerical 
structural  analysis  is  routinely  used  in  a  number  of  industries 
to  ensure  maximum  safety  and  reliability,  and  practical  appli¬ 
cations  of  enormous  size  have  been  accomplished  Cl]  for  the 
qualification  of  designs  for  which  full-scale  testing  is  imprac¬ 
tical  or  impossible. 

As  the  size  and/or  complexity  of  an  application  increases, 
however,  the  advantages  of  a  computer  simulation  are  sometimes 
reduced  due  to  the  time  and  cost  associated  with  the  preparation 
and  checking  of  input  data.  For  the  geometric  description  of  a 
finite  element  model,  the  tabulation  of  nodal  point  coordinates 
and  the  connections  of  elements  to  the  nodes  accounts  for  the 
bulk  of  the  numerical  data.  In  a  large  finite  element  model, 
manual  preparation  of  this  data  alone  may  consume  weeks  or  even 
months . 

This  report  describes  a  system  of  data  preprocessing  programs 
which  generate  and  manipulate  modeling  input  for  three-dimensional 


finite  element  structural  analysis.  The  preprocessing  system  is 
oriented  toward  the  construction  of  models  for  applications 
involving  thick  shell  and  solid  geometries,  although  axisymmetric 
and  two-dimensional  models  may  also  be  prepared.  Most  of  the 
facilities  of  the  preprocessor  are  concerned  with  the  geometric 
description  of  a  structure,  since  it  represents  the  largest  and 
most  complex  portion  of  the  problem  definition.  However,  the 
remaining  parts  of  the  model,  such  as  constraints,  loading,  and 
material  properties,  can  also  be  prepared  automatically.  Facil¬ 
ities  are  also  provided  within  the  present  system  of  programs 
for  interactive  plotting,  optimization  of  a  completed  model  for 
solution  efficiency,  and  communication  with  other  computer 
programs  which  perform  related  data  processing  functions. 

1.1  OVERVIEW  OF  PREPROCESSOR  FUNCTIONS" 

The  present  preprocessing  system  is  composed  of  a  series  of 
compatible  computer  programs,  which  perform  three  primary  classes 
of  functions: 

-  data  entry  and  translation, 

-  model  editing,  refinement  and  display,  and 

-  data  output  and  reformatting. 

These  functions  are  outlined  briefly  below,  and  in  detail  in 
Sections  2,  3  and  4  respectively. 
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The  motivation  for  developing  the  preprocessor  described  here 
is  the  current  need  for  a  means  of  constructing  finite  element 
analysis  models  for  aircraft  windshield- type  structures,  which 
can  be  broadly  characterized  as  solid  or  thick  shell  construc¬ 
tions.  As  a  result,  most  input  functions  are  strongly  oriented 
toward  problems  involving  curved,  three-dimensional  geometries 
which  might  be  described  by  surface  shapes  and  thickness  distri¬ 
butions,  by  lofting  coordinates,  or  by  mathematical  expressions 
in  curvilinear  coordinates. 

Data  entry  functions  performed  within  the  preprocessor 
provide  a  convenient  means  of  defining  the  geometry  of  this 
thick  shell  /  solid  class  of  structures;  analytically-defined 

shapes,  lofting  input,  coarse  surface  grid  data,  and  surface 
patch  input  are  supported,  by' means  of  file,  keyboard  and/or 
digitizing  pad  input.  "Superelement"  input,  which  defines 
a  geometric  shape  in  terms  of  the  coarsest  possible  finite 
element  mesh,  is  also  possible  for  general  two-  and  three- 
dimensional  regions.  For  models  prepared  using  other  data 
preprocessing  systems  (such  as  PATRAN-G  [2]  or  IMPRESS  [3]), 
data  translation  routines  are  available  for  converting  most 
finite  element  data  from  external  formats  to  the  data  forms 
understood  by  this  system. 

Every  data  entry  (or  translation)  function  results  in 
data  files,  which  are  stored  in  a  certain  internal  format 
for  use  in  the  preprocessor.  This  format  is  designed  for 
easy  translation  to  other  forms  as  needed  to  interface  with 
other  preprocessors  and/or  specific  analysis  programs. 


1.3 


Once  a  finite  element  model  is  stored  in  the  internal 
format  of  the  preprocessor,  any  number  of  operations  may  be 
performed  to  edit,  refine  or  otherwise  manipulate  the  model 
geometry  and  properties.  Two  or  more  models  may  be  merged 
to  form  a  single  geometric  model  at- this  stage  as  well. 

During  this  phase  of  preprocessing,  a  model  can  always  be 
saved  and  accessed  later  for  further  processing,  so  that  the 
modelling  process  need  not  be  accomplished  in  a  single  session 
at  the  computer  terminal. 

Final  version(s)  of  a  completed  finite  element  model  are 
translated  into  forms  understood  by  other  computer  programs 
in  the  data  output  and  reformatting  process.  The  finite 
element  data  can  be  formatted'  for  a  specific  analysis  program, 
or  rewritten  in  a  form  which  can  be  transferred  to  other 
preprocessor  formats  as  needed. 

The  interrelation  of  these  three  phases  of  preprocessing  is 
discussed  in  detail  in  the  following  subsection,  which  describes 
many  of  the  possible  data  paths  through  the  total  system.  Since 
the  preprocessing  system  is  constructed  from  a  group  of  separate 
programs  which  may  be  executed  in  any  order,  the  possibility  of 
introducing  specially-written  small  programs  to  perform  unusual 
or  unsupported  functions  also  exists.  This  open-ended  organiza¬ 
tion  can  be  a  great  advantage  for  the  experienced  user  who  is 
handy  with  FORTRAN,  while  less  experienced  (or  less  computer- 
oriented)  individuals  can  exercise  the  capabilities  of  the 
system  "as- is”. 
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1.2  PATHS  THROUGH  THE  SYSTEM 

The  organization  of  the  preprocessor  is  shown  in  its  most 
general  form  in  Figure  1.1.  Geometric  and  other  data  are  entered 
via  keyboard,  data  file  or  digitizing  pad,  using  one  of  the  data 
entry  or  interface  facilities  of  the  system.  At  the  opposite 
end,  a  completed  model  is  extracted '  from  the  system  format  and 
placed  in  a  usable  form  using  additional  data  translators.  In 
between,  the  main  preprocessor  (PREP)  may  be  used  to  modify, 
refine,  list,  plot  or  otherwise  manipulate  the  data. 

The  choice  of  a  data  entry  program  is  invariably  governed  by 
the  form  in  which  the  information  describing  the  model  is  most 
readily  accessible.  Table  1.1  summarizes  the  types  of  input 
which  are  accepted  by  each  of  the  data  entry  utilities  :  CREATE, 
IJKGEN,  CORGEN,  AGRID,  and  SPATCH. 

In  some  situations,  the  geometric  description  of  a  model  may 
already  exist  in  the  form  of  a  finite  element  mesh.  When  this  is 
the  case,  one  of  the  data  interfaces  TRNSFR,  IMPRINT,  or  NEUTRAL  . 
can  be  used  to  convert  the  information  to  the  preprocessor 
internal  format.  For  data  formats  not  recognized  by  the  system, 
a  set  of  general-purpose  data  translating  subroutines  are  pro¬ 
vided,  which  will  permit  the  conversion  of  external  data  with  a 
minimum  of  user  programming.  These  utility  routines  are  de¬ 
scribed  in  Subsection  3.3. 
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Figure  1.1.  Preprocessor  Organization. 


TABLE  1.1 


SUMMARY  OF  INPUT  DATA  TYPES  FOR  DATA  ENTRY  MODULES 


The  central  preprocessing  program,  PREP,  uses  a  single  file 
format  for  data  input  and  output.  Each  of  the  data  entry  and 
interface  utilities  creates  finite  element  data  files  in  the 
PREP  standard  format;  once  expressed  in  this  form,  the  finite 
element  mesh  can  be  reprocessed  numerous  times  using  PREP,  and 
may  be  combined  with  other  models  stored  in  the  same  form.  Data 
stored  in  the  preprocessor  internal  format  may  be  retained  on 
auxiliary  storage  and  modified  as  many  times  as  desired.  Printed 
output  and  geometry  plots  are  also  available  at  this  stage. 

Data  translators  which  interpret  the  output  of  the  prepro¬ 
cessor  are  of  two  types:  the  REFMT  processor  generates  finite 
element  data  which  is  readable  directly  by  the  finite  element 
analysis  program  MAGNA[4],  and  NEUTRAL  produces  a  formatted  data 
file  which  may  be  translated  for  use  with  other  preprocessors 
and/or  analysis  codes.  Since  NEUTRAL  can  also  be  used  to  trans¬ 
late  the  neutral  file  back  to  PREP  format,  it  is  also  a  useful 
means  of  transferring  modelling  data  from  one  computer  system  to 
another.  Obviously,  the  entire  output  translation  step  may  be 
replaced  if  necessary  by  one's  own  conversion  programs.  The 
PREP  data  sequence,  as  it  occurs  in  files  output  by  the  PREP 
processor,  is  described  fully  in  Section  6  for  such  applications. 


SECTION  2 


DATA  ENTRY  FUNCTIONS 


The  initial  step  in  generating  a  model  using  the  present 
preprocessing  system  is  concerned  with  the  general  descrip¬ 
tion  of  the  structural  geometry.  In  this  phase,  the  objec¬ 
tive  is  to  create  a  geometric  data  base  which  adequately 
describes  the  structure  or  substructure  of  interest.  The  model 
so  defined  will  be  generally  too  coarse  for  the  actual  finite 
element  stress  solution,  but  will  be  sufficient  to  identify  the 
structural  geometry  which  may  be  edited,  refined  and  merged 
with  other  model  files  to  result  in  a  final  model  which  will 
be  suitable  for  analysis. 

The  data  entry  modules  provided  by  the  preprocessing 
system  are  applicable  to  rather  arbitrary  structural  geometries. 
The  modules  discussed  in  this  section  detail  the  provisions  for 
utilizing  various  forms  of  model  data  input  to  define  a  'coarse' 
finite  element  mesh.  It  is  anticipated  that  the  user  will 
employ  the  editing  capabilities  of  the  verification  and  modifi¬ 
cation  modules  to  refine  the  coarse  mesh  to  a  mesh  suitable  for 
analysis. 

The  CREATE  data  entry  program  provides  the  simplest  possible 
form  of  input,  in  which  nodal  coordinates  and  element  connections 
are  defined  explicitly  in  free  format  at  the  keyboard.  This  mode 
of  input  can  be  useful  when  large  portions  of  a  model  can  be 
defined  in  terms  of  uncomplicated  shapes. 


The  data  generator  IJKGEN  allows  the  definition  of  a  model 
in  analytical  form,  using  known  equations  of  a  surface  or  other 
shape  to  set  up  a  topologically  regular  mesh  of  nodes  and 
elements.  Built  in  generation  options,  and/or  user-written  sub¬ 
routines,  may  be  used  in  generating  the  model. 

The  CORGEN  processor  assists  the  user  in  combining  pre¬ 
viously  defined  coordinate  data,  lofting  data  existing  either  in 
digitizing  form  or  numerical  form,  and  additional  data  that  may 
be  input  directly  from  the  keyboard.  This  program  makes  use  of 
an  intermediate  nodal  group  concept  to  allow  implementation  of 
interpolation  schemes  to  provide  for  regular  mesh  generation. 

AGRID  is  a  surface-fitting  utility  program  which  permits  the 
definition  of  a  regular  surface  mesh  from  scattered  data  which  is 
not  suitable  for  input  in  any  other  form.  This  module  uses  a 
special  smoothing  technique  which  computes  surfaces  passing 
through  a  number  of  specified  points  based  on  a  'least  curvature' 
criterion. 

The  SPATCH  utility,  in  conjunction  with  user-written  input 
routines,  facilitates  the  translation  of  geometric  data  stored 
in  the  form  of  bicubic  '  surface  patches'  into  the  forms  recog¬ 
nized  by  the  preprocessor. 

Each  of  t^ese  data  entry  functions  is  described  in  further 
detail  in  the  following  Subsections.  For  modelling  data  which 
is  incompatible  with  the  data  definition  utilities  discussed 
below,  the  reader  is  referred  to  Subsection  3.3,  which  describes 
a  set  of  general  interface  routines  which  can  be  adapted  to  the 
translation  of  most  types  of  data  into  the  preprocessor  formats. 


2.1  CREATE  -  Coarse  Mesh  Input 

The  CREATE  program  is  oriented  toward  problems  which  involve 
geometric  shapes  that  might  be  defined  by  a  very  coarse  grid  of 
finite  elements,  even  though  a  much  finer  mesh  would  be  required 
for  stress  analysis.  Such  a  simple  model  {e.g.,  a  tapered  plate 
defined  using  a  single  eight-node  brick)  can  be  translated  into 
the  preprocessor  format  using  CREATE,  and  refined  for  analysis 
with  PREP.  CREATE  is  fully  interactive,  prompting  the  user  for 
each  of  the  inputs  it  requires.  There  are  currently  limits  of 
500  nodes  and  100  finite  elements. 

The  access  procedures  for  CREATE  are  detailed  in  Appendix  C 
for  all  computer  systems  currently  supporting  CREATE.  Once 
initiated,  the  program  prompts  the  user  with  the  following 
pattern  of  requests: 

i.)  is  this  a  re-edit  session?, 

ii.)  define  a  2-D  or  3-D  coordinate  system  to  be  used, 

iii.)  input  the  nodal  coordinates, 

iv.)  edit  the  nodal  coordinates?, 

v.)  input  the  element  nodal  connectivity, 

vi.)  edit  the  element  connectivity?, 

vii.)  edit  nodal  data  or  element  connectivity? 

Sample  program  executions  are  provided  in  Section  7.1  for  CREATE. 

A  feature  is  provided  which  allows  for  re-editing  a  previously 
CREATE'd  file  to  eliminate  mistakes.  The  user  is  initially  asked 
if  the  current  session  is  a  re-edit  session.  If  the  user  responds 
that  it  is  a  re-edit  session,  there  must  be  a  file  present 
containing  previous  output  from  the  CREATE  program.  Should  the 
user  select  a  re-edit  session  there  will  only  be  one  edit  pass 
made  through  the  data. 

Data  entered  to  CREATE  must  be  either  two-dimensional  or 
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three-dimensional  in  form.  For  two-dimensional  data,  two  coor¬ 
dinates  (X,  Y)  must  be  defined  for  each  node;  three-dimensional 
coordinates  require  a  third  (Z)  value.  The  user  may  define  nodes 
that  are  not  connected  to  the  model  being  defined,  but  all  nodes 
specified  in  the  element  nodal  connectivity  must  be  defined. 
Unused  nodes  may  be  eliminated  in  the  PREP  module. 

Elements  are  defined  in  CREATE  with  reference  to  the  27-node 
solid  and  9-node  planar  elements  shown  in  Figures  2.1  and  2.2. 

If  the  user  has  selected  a  2-D  coordinate  system,  only  four  to 
nine  nodes  are  required  to  define  each  planar  element.  The 
elements  illustrated  in  Figure  2.2  are  typical  2-D  elements  the 
user  may  define.  When  specifying  elements  using  a  2-D  coordinate 
system,  the  user  must  input  the  element  number,  the  number  of 
nodes  to  be  used  to  define  this  element  (from  four  to  nine), 
followed  by  the  nodes  themselves.  Three  dimensional  models  are 
defined  in  the  same  manner  as  for  2-D,  except  that  each  finite 
element  may  contain  from  four  to  twenty-seven  connected  nodes. 
Since  each  element  may  have  a  different  number  of  nodes,  the  user 
must  enter  the  element  number,  the  number  of  nodes  to  be  used  to 
define  the  element  and  the  node  numbers  involved. 

The  completed  model  data  is  written  to  the  default  file 
specified  in  the  access  procedures  for  CREATE.  This  data  file 
is  unformatted  and  is  suitable  for  input  directly  to  the  PREP 
preprocessing  module  for  modification  and  reformatting. 


2.2  IJKGEN  -  ANALYTICAL  SURFACE  GEOMETRY  INPUT 


The  IJKGEN  mesh  generator  is  a  stand-alone  program  which 
creates  finite  element  data  from  analytical  descriptions  of  the 
geometry  in  question.  Figure  2.3  demonstrates  the  types  of 
geometries  which  lend  themselves  to  the  generation  scheme  used 
i..  IJKGEN;  in  each  instance,  the  shape  can  be  imagined  as  varying 


(at  least  piecewise)  continuously  with  the  "generator  indices" 

I,  J  and  K. 

A  very  simple  example  of  the  method  of  generation  in  IJKGEN 
appears  in  Figure  2.4.  Since  the  three  coordinates  R,  0,  Z  are 
orthogonal,  we  might  express  them  in  terms  of  generator  indices 
I,  J,  K  as 

R  =  R  (I) 

e  =  e  (j)  (2.1) 

Z  =  Z  (K) 

in  which  I  varies  from  one  to  the  number  of  nodes  in  the  radial 
direction,  and  so  on.  Given  the  limits  on  R,  0  and  Z,  IJKGEN 
would  generate  a  mesh  by 

-  looping  over  all  combinations  of  (I,J,K), 

-  for  each  (I,J,K)  combination,  "interpolating"  for  (R,0,Z), 

-  transforming  each  (R,0,Z)  to  Cartesian  (X,Y,Z)  coordinates. 

Connectivity  is  generated  for  the  model  using  the  fact  that  the 
shape  is  topologically  regular;  that  is,  the  connectivity  for 
the  cylinder  is  precisely  the  same  as  for  a  rectangular  parallele 
piped  having  the  same  number  of  nodes  in  each  direction  as  the 
cyl  inder . 

The  basic  mode  of  operation  of  IJKGEN  uses  a  number  of  built- 
in  options  for  geometric  shapes,  including  rectangular  flat 
plates,  cylindrical  shells  and  spherical  shells.  Furthermore, 
two  built-in  options  are  provided  for  interpolating  the  curvi¬ 
linear  coordinates  in  terms  of  I,  J  and  K:  (1)  uniform  mesh 
spacing;  and  (2)  logarithmic  grading,  with  the  ratio  of  the  sizes 
of  the  first  and  last  elements  in  each  direction  specified  by  the 
user.  Figures  2.5  and  2.6  illustrate  the  default  options  for 
coordinate  systems  and  mesh  spacing  in  IJKGEN. 


When  the  default  options  in  IJKGEN  are  used,  the  data  required 
as  input  includes 

-  the  coordinate  system  type  (  Cartesian  ,cyl  indrical , spherical )  , 

-  the  number  of  elements  in  each  direction, 

-  the  mesh  spacing  option,  and 

-  limits  on  the  curvilinear  coordinate  values. 

The  generated  mesh  is  always  composed  of  27-node  solid  elements, 
which  may  be  reduced  to  lower-order  elements  in  PREP  if  necessary 
(see  Section  4.1,  the  MASK  function).  To  ensure  that  27-node 
elements  which  are  to  be  MASKed  to  create  16-node  elements  (for 
example)  are  properly  oriented,  IJKGEN  also  allows  the  selection 
of  a  "thickness  direction"  prior  to  generating  the  element 
connectivity  data. 

Figure  2.7  shows  a  sample  execution  of  IJKGEN,  using  the 
default  options.  In  this  case,  the  shape  is  rectangular,  the 
mesh  is  graded  in  all  three  directions,  and  no  thickness  direc¬ 
tion  for  the  model  is  specified.  Figure  2.8  shows  a  plot  (ob¬ 
tained  from  PREP)  of  the  finite  element  model  as  output  by  IJKGEN. 

Naturally,  the  simple  geometric  shapes  included  in  IJKGEN  as 
built-in  functions  represent  the  exception,  not  the  rule,  for 
practical  applications.  When  the  shape  of  a  model  is  known  in  an 
analytical  form  other  than  those  included  in  IJKGEN,  three  user- 
written  subroutines  may  be  introduced  to  describe  the  geometry  in 
equation  form: 

-  SUBROUTINE  U INPUT 

-  SUBROUTINE  SURFAC  ( I ,J ,K, ALPHA, BETA, ZETA) 

-  SUBROUTINE  CRDTRN  ( ALPHA, BETA, ZETA, X,Y, Z) 

The  subroutine  UINPUT  is  called  at  the  beginning  of  execution,  to 
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tmtmmmumtmmmttniiiini 
nuitittt  teem  ijkgem  * sssssttss 
umtHiKmttuttmmmtttmntti 

I  J  K  d  E  N  -  GENERATION  Of  GEOMETRIC  MESH  DATA  FOR  SOLID, 
THICK  SHELL  OR  PLATE  FINITE  ELEMENT  MODELS,  USING  AN  INTEGER  - 
COORDINATE  INDEXING  SCHEME.  OPTIONAL  USER  ROUTINES  ARE  - 

(1)  SURF AC  ( I, J.K, ALPHA, BETA, ZETA)  -  DEFINE  MESH  GEOMETRY 

(2)  CRDTRN  (ALPHA, BETA, ZETA,X,V,Z)  -  COORD.  TRANSFORMATION 
BUILT-IN  OPTIONS  INCLUDE  RECTANGULAR,  CYLINDRICAL  OR  SPHERICAL 
COORDINATES,  AND  UNIFORM  OR  PROPORTIONALLY  GRADED  MESH  SPACING 

(3)  UlrSMJT  -  USER  PARAMETER  INPUT  ROUTINE  (INITIALIZE  DATA 

IN  BLANK  COMMON! 


lUUmtt  USER  SUBROUTINE  'SURF AC'  NOT  GIUEN  IIIIIIWl 

BUILT-IN  MESH  DIUISION  OPTIONS  ARE  AS  FOLLOWS  - 
(!>  -  UNIFORM  MESH  IN  EACH  DIRECTION 

(2)  -  GRADED  MESH  (SPECIFY  RATIO  Of  FIRST/LAST  ELEMENT  SI2E 

ENTER  OPTION  (  1  .  a  > . .2 

ENTER  THE  RATIO  OF  FIRST  /  LAST  ELEMENT  LENGTHS  FOR  EACH 
COORDINATE  DIRECTION  (ALPHA,  BETA.  ZETA)  (R-l  FOR  UNIFORM) 

ENTER  LENGTH  RATIOS  (R1,R2,R3)  . .345 

txnttl XXX  USER  SUBROUTINE  'CRDTRN'  NOT  GIUEN  XXXXXXXXt* 

BUILT-IN  COORDINATE  SYSTEM  TRANSFORMATION  OPTIONS  ARE  - 

<1!  RECTANGULAR,  (2)  CYLINDRICAL,  (3)  SPHERICAL 

ENTER  COORDINATE  SYSTEM  OPTION  (1.2,3)  -1 

tt  PLEASE  NOTE  THE  FOLLOWING  CONUENTlONS  FOR  RECTANGULAR  SYSTEM  SI 

'ALPHA'  •  X  'BETA'  •  V  'ZETA'  •  Z 
A  RIGHT-HANDED  SYSTEM  IS  ASSUMED. 

ENTER  LIMITING  SURFACE  COORDINATE  UALUES  - 


1.  ALPHA (MIN) 

3.  BETA  (MIN) 
S.  ZETA  (MIN) 


2.  ALPHA (MAX) 

4.  BETA  (MAX) 

5.  ZETA  (MAX) 


-S.  S.  t.  13.  -2.  1 


ENTER  THE  NUMBER  OF  CLEMENTS  TO  BE  GEN¬ 
ERATED  IN  THE  ALPHA,  BETA  AND  ZETA  CO¬ 
ORDINATE  DIRECTIONS,  RESPECTIVELY  ....  3  3  3 


»«  BEGIN  GENERATION  PHASE  It* 

J2J2  <£  TO  BE  GENERATED  *  343 

NUMBER  OF  ELEMENTS  TO  BE  GENERATED*  27 

JX,JMi2!f?,_®rtecTI0!L?r  THE  MODEL  (IF  ONE  EXISTS.  AS  IN  A 
THICK  SHELL)  fKiST  BE  IDENTIFIED  TO  ORIENT  ELEMENTS  PROPERLY? 

OPTIONS  ARE  ( t  > ALPHA,  (2)BCTA,  (3)ZETA,  OR  (4  KUNIPPORTANO . 

ENTER  THICKNESS  DIRECTION  CODE  (1.2.3. 4)  -4 

*M  DATA  GENERATION  COMPLETE  III 

mill  IJKQEM  TERMINATED  Qltlt 
STOP 

•2SSOO  MAXIMUM  EXECUTION  FL. 

_  >147  CP  SECONDS  EXECUTION  TIME. 


Figure  2.7.  Sample  Execution  of  IJKGEN. 
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Figure  2.8.  Plot  of  IJKGEN  Sample  Model  from  PREP. 


allow  initialization  of  parameters  (e.g.,  in  COMMON  blocks)  which 
might  be  needed  in  the  other  two  routines.  Subroutine  SURFAC 
computes  the  values  of  the  curvilinear  coordinates  (ALPHA, BETA, 
ZETA)  as  functions  of  (I,J,K),  and  therefore  allows  the  intro¬ 
duction  of  non-orthogonal  coordinates  and  variable  mesh  spacing. 
The  CRDTRN  routine  accepts  curvilinear  coordinate  values  (ALPHA, 
BETA, ZETA)  and  returns  Cartesian  coordinates  X,  Y,  Z;  in  effect, 
CRDTRN  defines  the  coordinate  system  in  which  IJKGEN  generates 
the  finite  element  mesh. 

Figure  2.9  shows  a  model  of  the  B-l  bomber  left 
windshield  generated  with  IJXGEN  [6].  In  this  case,  UINPUT  has 
been  used  to  read  in  "key  points"  for  the  windshield  geometry, 
which  is  cylindrical  with  irregular  boundaries;  a  two-element 
mesh  which  uses  only  the  key  points  and  a  few  intermediate  points 
is  shown  in  Figure  2.10.  SURFAC  is  used  to  generate  intermediate 
points  on  the  windshield  by  linear  interpolation,  and  CRDTRN 
performs  the  final  transformation  to  Cartesian  coordinates. 


Figure  2.10.  Coarse  B-l  Windshield  Mesh  using  Key  Point 
Coordinates. 
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2.3  CORGFU  -  General  Surface  Geometry  Input 

CORrnri  is  a  general  purpose  finite  eienent  preprocessor 
program  designed  to  interpret  surface  loftinn  tvpe  data,  either 
from  sets  of  coordinate  data  or  fron  digitized  data  in  the  for^ 
of  a  geometric  description.  The  pronran  will  also  accept  direct 
creation  of  nodel  geometry  and  eienent  connectivity  fron  keyboard 
input.  reonetric  data  input  nav  cone  fron  a  variety  of  nedia 
such  as  disc  files,  keyboard  or  digitizer  tablet,  or  a  conbination 
of  these.  Reference  points  entered  by  anv  of  these  neans  nay  he 
used  to  generate  data  alone  lines,  circular  arcs  or  nore  neneral 
curved  paths. 

Data  entry  is  performed  according  to  "croups"  of  points  where 
data  fron  an  individual  lofting  plane  (cross-section)  will  cor¬ 
respond  to  a  sinnle  group  as  illustrated  in  Finure  2.11.  For  each 
group  of  data  input,  the  user  nay  select  a  different  input  nediun 
(e.g  file,  keyboard  or  digitizer),  a  separate  coordinate  systen 
transfoma tion  (Cartesian,  cylindrical,  or  spherical),  independent 
coordinate  translations  and  other  parameters  which  define  the 
placement  of  the  group  within  the  global  coordinate  systen.  Input 
fron  the  digitizing  tablet,  for  instance,  is  always  given  in 
planar  form  (2-D)  but  nay  be  arbitrarily  located  in  three  dimen¬ 
sional  space  hv  specifying  the  location  of  a  positional 
digitizing  plane  using  any  three  non-collinear  points  lying  in 
the  plane. 

Ultimately,  the  geonetric  data  presented  to  CORCEP  must  be 
arranged  in  nodal  groups  containing  similar  numbers  of  points. 

Data  which  does  not  conform  to  this  scheme  may  be  interpolated  or 
otherwise  modified  after  being  input  to  CORGFP.  The  commands 
available  in  CORCFT1  are  generally  of  three  types:  (1)  data  input 
commands;  (2)  generation/intorpolation  commands;  and  (3)  element 
generation  commands.  Input  conmands  define  an  input  medium,  or 
cause  data  to  be  read  fron  an  al ready-def ined  medium.  Generation 
and  interpolation  commands  initiate  internal  generation  of  nodes, 


usually  for  the  purpose  of  "fillinn  out"  one  or  more  nodal  nrouns. 
Element  qeneration  commands  cause  the  surface  model  to  be  completed 
by  defininq  surface  element  connection  data.  - 

Often  the  available  data  is  not  suitable  for  qeneratinq  a  finite 
element  mesh  directly,  due  to  poor  point  spacinq  or  density.  In  such 
situations,  interpolation  of  point  data  may  be  performed  in  CORCTM, 
usino  the  method  described  in  Reference  [5],  Interpolation  of  new 
nodes  between  existinq  nodes,  or  placement  of  entire  new  node  qroups 
between  nreviously-def ined  cross-sections  produces  smoothly  interpol¬ 
ated  data,  even  in  situations  where  spline  fittinq  may  be  inadequate. 
Fiaure  2.12  illustrates  the  salient  differences  between  these 
me  thods . 

Whenever  the  user  specifies  the  input  of  nodal  coordinate 
data,  prompts  will  request  necessary  information  for  forminq 
qroups  from  that  data.  For  example,  if  the  user  were  to  specify 
file  inout  of  nodes,  the  user  would  be  asked  to  snecifv  the 
number  of  nodes  per  qroup  to  be  input  as  well  as  the  number  of  • 

qroups  that  should  be  input  concurrently  containinq  that  number 
of  nodes  per  qroup,  before  returninq  control  to  the  user. 

Simple  nodellinq  requirements  will  entail  no  special  interaction 
by  the  user  to  alter  or  select  qroups,  since  CORCEM  is  structured 
to  assume  the  user  has  entered  nodes  and  qroups  in  the  order 
they  are  to  be  utilized  in  element  qeneration. 

Element  qeneration  in  COROEH  defaults  to  utilize  all  croups 
currently  defined,  in  the  order  they  were  defined.  This  may  be 
altered  by  the  user  in  special  cases  where  new  qroups  are  created 
or  qroups  were  not  input  in  the  correct  order  for  the  nenerator. 
Additionally,  the  user  should  be  aware  that  when  a  nroun  is 
specified  for  in  terno  la  t  ion ,  the  new  croup  will  replace  the  old 
qroup.  This  can  result  in  sienificant  disarray  of  node  nunberinq 
sequences  in  the  completed  model.  Provisions  are  made  for  renum¬ 
ber  inq  the  nodes  in  the  nREP  module  to  eliminate  problems  with 
non-optimal  node  nunberinq.  Should  the  user  create  qroups  that 
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arc  not  built  with  eauivalent  nunhers  of  nodes,  CORC.nn  i  1 1 
automatically  invoke  the  interpolation  routines  to  honoqenizo 
all  qroups  specified  for  elenent  qeneration.  Element  qeneration 
has  two  options  available  for  positioninq  midside  nodes.  The 
proqram  default  is  to  always  aenerato  interpolated  nodes  betwaen 
existino  nodes  such  that  midside  nodes  are  properlv  located. 

Occasionally,  it  is  necessary  to  extensively  aJ  ter  input 
qeonetry  data  to  achieve  an  appropriate  model  for  analysis. 
Provisions  are  available  for  user  creation  of  new  croups  by 
snecifvinq  nodes  to  be  used  for  interpolation  and  by  user  selec¬ 
tion  of  qroups  to  be  utilized  in  elenent  oeneration.  Addition¬ 
ally,  a  number  of  options  exist  for  alterinq  node,  croup  and 
elenent  data  to  meet  a  variety  of  nodellinq  requirements.  In 
the  event  the  user  wishes  to  utilize  the  elenent  qenerator  but 
does  not  wish  to  have  new  nodes  interpolated  then  the  croups 
must  be  defined  such  that  an  odd  number  of  nodes  are  located  in 
each  qroup  and  that  an  odd  number  of  croups  are  selected  for 
elenent  qeneration. 

2.3.1  Interactive  Proqram  Execution 

The  interactive  capabilities  of  C0RGFT1  are  summarized  in 
Table  2.1.  The  proqram  is  essentially  command  driven,  but 
oriented  such  that  the  user  is  displaved  only  a  menu  of  those 
options  or  commands  immediately  relevant.  Txtensive  help 
facilities  assist  the  user  in  understand inc  what  options  are 
available  and  provides  information  relevant  to  the  current  stace 
of  model  development. 

Initial  input  to  COROEN  consists  of  establishinc  the  proper 
coordinate  system  transformations,  inplenentinq  anv  coordinate 
translations,  resettinq  default  nodal  attributes,  and  any  other 
parameters  that  affect  the  input  of  nodal  coordinate  data. 

Fiqure  2.13  lists  the  default  options  of  CORCEN  and  what  commands 
are  used  to  alter  them.  Once  the  user  has  adeouately  defined 


TABLE  2.1 


CORGEN  INTERACTIVE  PROCESSING  CAPABILITIES 


Node 

Input 

File  X 

Keyboard  X 

Digitizer  X 

Editing 

File 

Keyboard  X 

Digitizer  X 

Interpolation  X 

Auto-generation 

Auto-interpolation  X 

Generation  (other)  X 

Save/Restore  X 


Group 


X 


X 

X 

X 


X 


Element 


X 

X 


X 


X 


X 


Other  features: 

o  2-D  lofting  plane  specification  in  arbitrary  space 
o  Default  values  for  thickness  and  coordinates 
o  Coordinate  transformations  and  translations 


Default  Options  for  CORGEN  Program  Execution 


Function 


Default 


Command 


Auto  Interpolation 
Auto  Mesh  Generation 
Coordinate  System 
Digitizer  Pick  Tolerance 
Lofting  Plane 


Translation  Factors 
Thickness 


Yes 

Yes 

Cartesian 
5%  horizontal  axis 
(0,0,1) 

X  vertical  screen 
Y  horizontal  screen 
Z  into  screen 
x=0 ,  y=0 ,  z=0 
1.0 


INTErpolate 

ELEMent 

COORdinate 

TOLErance 


PLANe 

COORdinate 

THICk 


Figure  2.13.  CORGEN  Default  Options  and  Commands 


Menu  Paths 


A.  File  input  of  nodal  coordinates,  automatic  interpolation, 
and  mesh  generation: 

NODE,  FILE,  READ,  FINI,  STOP. 

B.  Keyboard  input  of  coordinates,  order  groups  for 
interpolation  then  selective  element  generation: 

NODE,  KEYB ,  READ? 

GROUP,  ALTER? 

ELEM,  GENER? 

STOP. 


Figure  2.14.  CORGEN  Command/Menu  Paths. 
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the  input  coordinate  system,  a  mechanism  should  be  selected  for 
nodal  coordinate  input.  In  defining  the  nodes  to  COP.OEN,  the 
user  must  be  coqnizant  of  the  need  to  incorporate  nodes  into 
'nodal  qroups'  .  All  nodal  input  sequences  will  question  the 
user  as  to  how  many  nodes  are  to  be  placed  in  the  next  input 
qroup.  There  should  qenerallv  be  an  equal  number  of  nodes  in 
all  croups  or  the  interpolation  functions  will  automatically  be 
invoked  to  reconfiqure  the  qroups  so  that  all  are  equal. 

The  menu  paths  of  C^R^ET!  are  detailed  in  riqure  2.1*.  The 
user  is  encouraqed  to  follow  a  str ainhtforward  approach  to  rodel 
qeneration  by  first  defininq  nodes  in  nodal  qroups,  then  selection 
any  needed  interpolation  and  finally  qeneratinq  elements.  Many 
of  the  features  in  CORGEN  are  automatic  and  will  be  invisible  to 
the  user  if  nodal  qroups  are  defined  properly  at  the  outset. 

Each  master  command  will  be  discussed  below.  All  commands  are 
enabled  by  typinq  the  one  to  four  letter  abbreviation  enclosed 
in  parentheses  in  the  commands  listinq.  The  command  listinqs 
are  alwavs  available  by  typinq  COMMands.  Master  commands  are 
those  commands  which  are  available  for  the  COMMAND,.:  prompt. 


2.3.2  COMMANDS  LISTING 

The  user  has  available  a  number  of  command  options  at  all 
times.  The  proqran  is  structured  to  display  only  those  commands 
which  are  immediately  neccessary  to  what  the  user  is  doinq  based 
on  previous  commands  entered.  The  master  commands  are  listed  in 
Fiqure  2.15.  The  master  commands  will  provide  the  user  with 
more  specific  options  to  accomplish  particular  tasks  bv  renues- 
tinq  more  specific  input  or  data.  The  COMMands  command  will 
always  display  the  current  optional  commands  available  to  the 
user.  Master  commands  are  always  accessible  bv  the  user, 
unless  specific  numeric  data  has  been  requested  such  as  nodal 
coordinate  values.  The  proqran  default  is  to  provide  a  listinq 
of  the  available  commands  only  when  the  user  types  COMMands. 
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CORGEN  MASTER  COMMANDS 


COMMANDS. 


( COMM) AND S  LISTING 
(HELP) 

(NODE)  DATA 
(ELEM)ENT  DATA 
(GROU)P  DATA 
(SAVE)  MODEL 
(REST)ORE  MODEL 
(COOR)DINATE  SYSTEM 
(INTE)RPOLATE  DATA 
(DISP)LAY  VALUES 
(FINI)SH  MODEL 
(THIC)KNESS 
(STOP) 


Figure  2.15.  Master  Commands  for  CORGEN. 
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2.3.3  COORDINATE  SYSTEM 

Three  autonatic  coordinate  transfomations  are  built  into 
CORGEM.  The  Cartesian  (rectilinear)  coordinate  systen  ir 

rnt’lt ,  c^-r''5.  iorito  cnnvorr>  ioors  *"■*.?  *  ?  ! r,v  o  c  i  r  \  "  ot* 

r'ata  input  an  (KrPPI,n,!F?A)  where  R  is  the  radius  of  a  sr>here  and 
THETA  and  PHI  represent  annlos  in  decrees,  locatinq  the  points  on 
a  spherical  surface.  Cylindrical  coordinates  conversion  is  also 
available  for  nodel  data  defined  by  (R, THETA, Z)  coordinates, 
where  R  is  the  radius  fron  the  center  of  a  cylinder,  THETA  is 
the  ancle  in  decrees  of  the  point  about  the  circular  axis  and  Z 
defines  the  loncitudinal  distance  alonq  the  cvlinder  axis.  The 
three  coordinate  systems  are  illustrated  in  Finure  2.16.  '’’he 

COORdinate  command  selects  the  current  coordinate  svsten  until  a 
new  one  is  specified. 

The  user  is  queried  if  coordinate  translations  of  data  from 
the  oriqin  are  desired  whenever  the  COORdinate  comnand  is  executed 
The  translation  values  will  remain  in  effect  until  chanced  or 
disabled  by  aqain  executinc  the  COORdinate  comnand.  Coordinate 
translations  are  specified  as  positive  or  neqative  macnitudes  to 
direct  the  node  in  the  X,  Y  and  Z  directions  respectively/ .  The 
translations  are  applied  to  input  coordinate  data  after  the  coor¬ 
dinate  transformations. 

2.3.4  DIGITIZING  DATA  COMMAND 

DIGItize  is  an  input  command  desiqned  for  dicitizinc 
nodes  and  element  connectivity  utilizinc  loftinn  planes 
oriented  arbitrarily  in  three  dimensional  space.  nata  may 
be  input  to  CORCEM  in  the  form  of  dicitized  points  from 
sources  such  as  loftinn  dravinqs. 

The  oxritizc  comnand  is  currenty  implemented  to  operate 


Y 


a.  Right-handed  Cartesian  Coordinate  System  Definition  (X,  Y,  Z) 


Y 


b.  Right-handed  Spherical  Coordinate  System  Definition  (R,<£,0) 


Y 


c.  Right-handed  Cylindrical  Coordinate  System  Definition  (R,  0,  Z) 


CORGEN  COORDINATE  SYSTEMS 

Figure  2.16.  Right  Handed  Cartesian  (X,Y,Z)  Coordinate  System; 

Right  Handed  Spherical  (R,4»f0)  Coordinate  System; 
Right  Handed  Cylindrical  (R,9,Z)  Coordinate  System. 
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with  'Tektronix  qranhics  tablets,  teminal  (cross-hair)  screen 
cursors,  or  terminals  that  enulate  Tektronix  qraphics  terminals. 
The  command  uses  both  command  prompt  and  question/answer 
techniques  to  establish  the  dinitizinq  area.  The  user  will  be 
pronpted  initially  for  input  to  identify  the  type  of  eraphics 
terminal  beinq  used;  the  type  of  diqitizinn  to  be  done  (screen 
cursor  or  tablet);  the  location  of  a  menu  area  if  one  is  desired; 
the  definition  of  a  loftinn  plane;  and  the  orientation  of  the 
drawinn  to  be  dioitized  on  the  tablet.  Once  these  initial  iters 
are  defined  the  user  is  required  to  bonin  diqitizinq  coordinate 
points.  After  the  user  has  completed  input  of  points  for  the 
first  loftinq  plane,  a  menu  of  command  options  is  displayed  to 
allow  the  user  to  redefine  the  loftinq  plane;  add,  move  or  delete 
nodes;  create  or  edit  elements;  or  to  reset  default  values  of 
parameters  such  as  layer  thicknesses  or  coordinate  translations. 
Fiqure  2.17  illustrates  the  initialization  information  required. 

Uhen  the  user  selects  the  tablet  input  option,  the  prooran 
must  request  information  to  locate  the  model  on  the  tablet  and  to 
set  a  user  coordinate  system  independent  from  that  on  the  tablet. 
In  this  case  the  user  must  first  diqitize  the  point  that  will 
represent  the  oriqin  (0,0)  of  the  first  drawinq  to  be  input. 

This  point  is  followed  by  a  second  diqitized  point  which  should 
define  the  maximum  point  alonq  the  USPR’S  horizonta'  axis. 

The  user  should  note  that  the  user’ s  horizontal  and  vertical 
axes  do  not  need  to  alien  with  the  tablet  (or  terminal)  hori¬ 
zontal  axes  as  the  proqram  will  provide  the  appropriate  coor¬ 
dinate  translations. 

An  optional  ’menu'  area  is  allowed  to  facilitate  proqram 
execution.  This  menu  consists  of  lf>  blocks  arranqed  in  a  4  by  4 
matrix  as  illustrated  in  Fiqure  2. in.  Uhen  the  menu  area  option 
is  selected  it  allows  the  user  to  nuickly  exit  the  diqitizinn 
mode  or  to  easily  undo  any  operation  the  user  may  have  not 
really  wanted  to  perform.  In  addition,  the  menu  area  can  be  used 
to  desianate  any  one  of  up  to  nine  different  predefined  thickness 


PLOTTING  PACKAGES: 

1  -  DI-3000 

2  -  TEK  PLOTIO 

3  -  NO  PACKAGE 

ENTER  PACKAGE  TYPE  ( 1-3 2 
TERMINAL  TYPES: 

1  -  NONGRAPHICS 

2  -  TEK  4006 

3  -  TEK  4010,  4012,  4052 

4  -  TEK  4014,  4016 

5  -  TEK  4014  EG 

6  -  TEK  4114 

ENTER  TERMINAL  TYPE...:  5 
DIGITIZER  TYPES: 

1  -  TERMINAL  W/ CURSOR 

2  -  TERMINAL  W/TABLET 

3  -  TABLET  STAND-ALONE 

ENTER  DIGITIZER  TYPE  (1-3)....:  2 


ESTABLISH  A  TABLET  MENU  AREA  FOR 
CONTROLLING  DIGITIZING?  ( Y,N)  .  .  :  Y 
THE  USER  MUST  DIGITIZE  THE  LOWER 
LEFT  AND  UPPER  RIGHT  CORNERS  OF 
THE  TABLET  MENU  AREA,  THEN  ENTER 
THE  NUMBER  OF  HORIZONTAL  AND 
VERTICAL  MENU  OPTIONS  AVAILABLE. 
DIGITIZE  THE  LOWER  LEFT  MENU  AREA..: 

DIGITIZE  THE  UPPER  RIGHT  MENU  AREA..: 


TABLET  -  SET  USER  ORIGIN  AND  AXES 
TO  ORIENT  THE  MESH  AXES  ON  THE 
TABLET.  THE  USER  MUST  DIGITIZE 
TWO  POINTS  -  THE  HORIZONTAL 
ORIGIN  OF  THE  MODEL,  THEN  THE 
MAXIMUM  HORIZONTAL  VALUE  ON  THE 
HORIZONTAL  AXIS.  THE  USER-ASSIGNED 
COORDINATES  FOR  THESE  TWO  POINTS 
MUST  THEN  BE  ENTERED  TO  PROPERLY 
SCALE  THE  DIGITIZED  DATA.  DIGITIZE 
THE  ORIGIN  OF  THE  HORIZONTAL  AXIS.: 

ENTER  THE  COORDINATES  FOR  THIS  POINT(X,Y).: 
DIGITIZE  THE  MAXIMUM  HORIZONTAL  POINT: 

ENTER  THE  COORDINATES  FOR  THIS  POINT(X,Y).: 


Figure  2.17.  DIGItize  Command  Initialization 
Questions . 


0.  ,0. 
100. ,0. 


SCALE  IS:  1  USER  UNIT  =  65.300  UNITS. 


THE  CURRENT  LOFT  PLANE  IS  DEFINED 
BY  THE  FOLLOWING  VECTORS.: 


VI  = 

( 

1.0000, 

0.0000, 

0.0000) 

V2  = 

( 

0.0000, 

0.0000, 

0.0000) 

V3  = 

( 

0.0000, 

1.0000, 

0.0000) 

CHANGE  THIS  LOFT  PLANE?  (  Y ,  N)  .  .  .  .  :  N 

DO  YOU  WISH  TO  INPUT  THICKNESS 
VALUES? (Y,N)....:  Y 

WHICH  THICKNESS  OPTION? 

1  -  INPUT  A  THICKNESS  VALUE  FOR  EACH  POINT 

2  -  SET  ONE  THICKNESS  FOR  ALL  POINTS 

3  -  DEFINE  LAYERED  THICKNESS  INDICES 

ASSOCIATED  WITH  EACH  POINT  OR  GROUP 
ENTER  THICKNESS  OPTION  (1,3)....:  2 

ENTER  THE  OVERALL  THICKNESS  (>0)....:  0.75 

THE  PICK  TOLERANCE  ■  4.00 

DO  YOU  WANT  TO  CHANGE  IT?  ( Y,N) . :  N 
USER  MUST  DIGITIZE  NODES  BEFORE  ANY 
OTHER  COMMAND  MAY  BE  EXECUTED. 

DIGITIZE  NODES.  WHEN  FINISHED 
DIGITIZE  THE  ’R'  CHARACTER  IN 
THE  MENU  BOX,  OR  DIGITIZE  THE 
SAME  POINT  3  TIMES.  DIGITIZE 
THE  'D'  TO  DELETE  THE  LAST  NODE 
INPUT. 

DISPLAYED  POINTS  ARE  THOSE 
ALREADY  PLOTTED  BY  USER. 

BEGIN  DIGITIZING  POINTS: 

READY?  — >  Y 


Figure  2.17.  (Continued) 
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values  to  each  node  as  it  is  dioitizeri,  in 

only  a  single  thickness  value  for  all  node 

blocks  are  defined  to  be: 

lieu  of 

s  input. 

pr  ed  e  f  i  n  i  nn 

The  If, 

FLOCK  LAPP  L 

COWAN  n 

1-9 

0 

R 

r 

P 

• 

Pinnon 

thickness 

ind  ex 

( predefined ) 

none 

return 

return 

delete 

last 

node 

none 

mcnr.n 

none 

none 

return 

return 

d  el  etc 

1  ast 

element 

none 

rrovnop 

not  available  -  user  oust  dioitize 

three  tines  to  exit  this  connand . 

the  sane 

po  int 

DELNOn 

none 

none 

return 

return 

restore 

last 

node 

none 

rr  it.  lm 

none 

none 

return 

return 

restore 

none 

1  ast 
element 


LOF*”  —  not  used  — 


Note  that  the  "0"  and  "R"  blocks  appear  twice  each  in  the 
menu  area,  as  shown  in  Figure  2.18.  The  menu  area  may  be 
located  anywhere  within  the  available  digitizing  area,  and 
may  be  set  to  any  size  which  is  convenient  for  the  user.  The 
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specified  location  of  the  menu  area  is  qiven  by  first  d  in  i- 
tizinq  the  lover  left  corner  of  the  nenu  block,  then  diqi- 
tizina  the  upper  riqht  corner  of  the  nenu  block.  The  pronran 
will  display  the  nenu  block  when  it  enters  the  PICNOP  connand . 
Each  tine  the  PSPLAY  connand  is  executed  the  nenu  area  vill  be 
redrawn.  To  select  any  i  ten  in  the  nenu  area  the  user  need 
only  place  the  dinitizinn  stylus  within  the  block  surroundinq 
the  character  desired. 

A  default  loftinq  plane  is  defined  which  places  the  qlobal 
(3-D)  X  axis  alonq  the  local  (2-D)  Y  axis  (vertical  screen); 
the  qlobal  Y  axis  lies  alonn  the  local  X  axis  (horizontal  screen 
or  tablet);  and  the  nlobal  E  axis  is  positive  into  the  terninal 
screen  ( away  fron  the  user).  This  is  illustrated  in  Fiqure  2.10 
The  user  nay  alter  the  loftinq  plane  by  specifyinq  any  three 
noncollinear  points  that  lie  within  the  desired  plane.  The 
LOFT  connand  contains  all  further  discussion  of  this  ootion. 

The  user  is  qiven  the  option  of  specifyinq  thickness  value 
to  be  associated  with  each  point  input.  Three  options  exist  for 
specifyinq  the  thickness  values.  First,  the  user  nay  input  a 
thickness  value  for  each  point  AS  THEY  ARC  BEING  DIOITIZCP. 

The  proqran  will  pronpt  the  user  for  each  coordinate  thickness 
value  after  it  has  been  diqitized.  A  second  option  allows  the 
user  to  specify  a  sinqle  thickness  value  to  be  assiqned  to  all 
nodes  as  they  are  dinitized.  This  sinqle  thickness  value  vill 
renain  in  effect  until  channed  by  the  user  vi th  the  THICK 
connand.  The  third  option  al lows  the  user  to  define  up  to  nine 
different  thickness  values.  These  values  are  indexed  in  the 
order  thev  are  input  to  ihe  proqran  (e.q.  first  value  =  index  1) 
and  are  referenced  by  their  index  value.  Once  a  point  is  input 
by  the  user,  an  index  (nenu  area  nunbers  1-9)  nust  then  be 
diqitized  to  define  the  thickness  at  that  node. 

To  connlete  the  initialization,  the  proqran  displays  the 
current  pick  tolerance  (resolution)  in  user  units  and  queries 


the  user  if  a  chanqe  is  desired .  A  smaller  pick  tolerance 
will  require  the  user  to  be  more  exact  in  selectinq  nodes  or 
elements  for  editinq  or  deletinq.  It  is  recommended  that  the 
user  select  a  smaller  tolerance  if  the  coordinates  to  be  input 
lie  relatively  close  to  each  other  (i.e.  within  5%  total  distance). 

Once  the  initialization  is  completeo  t^e  proqran  will 
automatically  execute  the  Pirnop  command  .  The  PlCNOn  command 
nav  be  exited  at  any  time  into  the  command  mode.  Finure  2.20 
lists  the  various  commands  available  in  Digitize  and  a  list 
of  the  F!nLn  command  output.  Fach  command  is  discussed  below. 

nimriop  allows  the  user  to  specify  the  input  of  points. 

The  user  also  must  define  a  thickness  value  for  each  point  input. 

As  discussed  above  there  are  several  options  that  allow  the  user 
to  select  a  constant  thickness  value  for  all  points;  to  specifv 
via  the  keyboard  a  different  thickness  for  each  point  input;  or 
a  means  of  specifvinq  two  to  nine  thickness  and  diqitizinq  the 
index  of  one  of  those  to  associate  with  the  current  point.  mhe  • 

user  must  have  selected  the  menu  option  in  order  to  use  the  indexed 
thickness  values  options.  Once  the  user  has  defined  enouoh  nodes 
he  may  exit  by  diqitizinq  the  same  point  three  times  in  succession 
or  by  pickinq  the  1 R'  or  ' C'  blocks  in  the  menu  area.  If,  at 
any  time,  the  user  wishes  to  delete  the  last  node  entered, 
simply  diqitize  the  ' D'  block  in  the  menu  area.  Repeated  use 
of  the  • D'  block  will  continue  to  delete  nodes  from  the  hiqhest 
numbered  node  down. 

The  DIOFLF  command  allows  the  user  to  specify  element 
connectivity  either  by  kevboard  input  or  by  pickinq  currently 
defined  nodes  with  the  diqitizer.  At  least  four  nodes  must  be 
specified  to  create  an  element  and  as  many  as  nine  nodes  mav 
be  qiven.  Should  the  user  desire  to  utilize  the  dioitizer  to 
delete  elements,  then  at  least  five  nodes  must  be  qiven  where 
the  fifth  node  corresponds  to  local  node  number  9.  Piqure  2.21 
illustrates  the  order  nodes  are  to  be  qiven  to  define  an  element 
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ENTER  COMMAND:  HELP 


CONTROL  COMMANDS 


COMMAND 

FUNCTION 

DIGNOD 

- 

DIGITIZE  NODES  (REQUIRED  FIRST) 

DIGELE 

- 

DIGITIZE  ELEMENT  CONNECTIVITY 

DELNOD 

- 

DELETE  A  NODE 

DELELE 

- 

DELETE  AN  ELEMENT 

DSP LAY 

- 

DISPLAY  ALL  DIGITIZED  POINTS 

HELP 

- 

PRINTS  THIS  LIST 

MENU 

- 

CHANGE  MENU  LOCATION 

MOVNOD 

- 

MOVE  A  NODE  TO  A  NEW  POSITION 

LOFT 

- 

DEFINE  A  NEW  LOFTING  PLANE 

SAVE 

- 

SAVE  MODEL  DATA 

STOP 

STOP  PROGRAM 

Figur*  2.20.  DIGItize  Command  HELP  Feature 
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either  by  keyboard  or  by  digitizer.  All  four  corner  nodes  are 
required  to  define  each  element,  while  the  remaininq  five  nodes 
are  optional.  The  user  must  specify  the  hiohest  local  node 
number  to  be  defined  as  the  total  number  of  nodes  to  be  aiven 
for  an  element  connectivity.  For  example,  if  the  user  wishes  to 
qive  5  node  nunbers  to  define  an  element,  but  have  then  corre¬ 
spond  to  local  node  locations  1-^  and  9,  the  connectivity  list 
will  be  1 , 2, 3 , 4 , 0 , 0 , 0 , 0, 5  where  all  intermediate  nodes  that  are 
not  defined  must  be  set  to  0 .  If  only  4  nodes  are  to  be  aiven 
the  user  nay  -just  specify  those  4  and  the  rest  will  be  set  to  0. 
To  exit  from  the  DICELC  command  the  user  may  diqitize  the  same 
point  three  tines  or  pick  the  '  R'  or  ' C '  blocks  in  the  menu 
area  (if  the  diqitize  option  is  selected). 

2.3.5  DISPLAY  VALUES 


The  user  nay  execute  the  DISPlay  connand  to  obtain  information 
concerninq  the  status  of  the  model  beinq  defined  and  user  defin¬ 
able  parameters  that  affect  model  data  definition.  Fiqure  2.22 
illustrates  the  output  provided  by  the  DISPlay  command.  All  model 
parameter  values  are  calculated  in  response  to  the  command,  so 
repeated  use  of  the  command  can  be  wasteful.  CORGEN  will  automat¬ 
ically  provide  this  information  before  proqran  termination. 


2.3.6  CLEMENT  DATA 

The  ELEMent  command  provides  the  user  with  capabilities  for 
inplenentinq  the  nesh  qeneration  features  of  CORCr*!.  Options 
exist  for  selectinq  alternate  groups  of  nodes  to  be  utilized  in 
the  element  generator  or  for  the  user  to  enter  new  elements 
directly  from  the  keyboard.  Editing  features  are  provided  which 
allow  the  user  to  alter  nodes  specified  in  element  connectivity 
or  to  replace  or  delete  elements.  Elements  can  be  listed  by  the 
selection  of  a  range.  Figure  2.23  illustrates  the  command  options 
available  with  the  ELEMent  command.  The  nine  node  biquadratic 
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MODEL  SUMMARY  DATA: 

NUMBER  OF  NODES  DEFINED . : 

HIGHEST  NUMBERED  NODE . : 

NUMBER  OF  UNDEFINED  NODES . : 

NUMBER  OF  GROUPS  DEFINED . : 

HIGHEST  GROUP  DEFINED . : 

NUMBER  OF  UNDEFINED  GROUPS.  . . : 

NUMBER  OF  ELEMENTS  DEFINED .  .  .  .  ; 

HIGHEST  ELEMENT  DEFINED . : 

NUMBER  OF  UNDEFINED  ELEMENTS . . : 

NUMBER  OF  ELEMENT  REFERENCES 
TO  UNDEFINED  NODES . : 

m* 


Figure  2.22.  DISPlay  Command  Output  listing 
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I 


ELEMENT  -  specify  element  generation  editing 

ALTER  -  change  element  connectivity 
GENERATE  -  generate  elements  using  selected  groups 
FILE  -  input  element  connectivity  from  file 
KEYBOARD  -  input  element  connectivity  from  keyboard 
HELP  -  help  user 

INTERPOLATE  -  refine  mesh  by  interpolation 
LIST  -  list  elements 
DELETE  -  delete  elements 

TABLE  -  give  summary  of  element  statistics 


Figure  2.23.  ELEMent  Command  Options. 


element  illustrated  in  Figure  2.21  is  the  element  type  qenerated 
by  the  ELEMent  command.  The  nronram  allows  the  user  to  specify  a 

element  connectivity  directly,  however,  the  user  must  specify  — 3 

elements  that  correspond  to  the  local  node  nunberinq  qiven  in 
Fiaure  2.21  for  the  quadratic  surface  element. 

I 

If  the  INTFrpolate  option  has  not  been  invoked  at  the  time 
elements  are  qenerated,  CORCEN  permits  all  existinn  nodal  data  to 
be  interpolated  (refined)  in  both  surface  directions)  prior  to 
qeneratinq  elements.  Invokinq  the  interpolation  option  at  this 
time  produces  twice  as  many  elements  in  each  direction  (refinement 
may  also  be  postponed  until  PREP  is  executed;  see  Subsection  4.1). 

2.3.7  FINISH  MODEL 

Once  the  model  nodal  coordinate  data  has  been  input  properlv 
to  form  appropriate  qroups,  the  user  may  select  the  Finish 
command  to  complete  all  element  qeneration  and  preprocessor  •■'"i 

functions  to  conclude  the  proqram.  This  command  will  access  the 
current  qroups  list  whether  established  by  default  or  selected 
by  the  user,  and  invoke  any  necessary  interpolate  -  r '•endures 
before  oxoeutinn  the  element  venerator.  The  command  writes  an 
unformatted  data  file  which  the  user  may  input  directly  to  the 
EXPAND  nodule  to  obtain  the  3-D  solid  element  model  needed  for 
input  to  the  PREP  pronram.  This  command  will  provide  a  listinq 
of  the  model  parameters  illustrated  in  Fiqure  2.22  once  the 
model  has  been  completed.  Included  in  this  listinq  will  be  the 
file  number  contain  i-nn  the  current  model  data.  The  FTOP  command 
has  the  same  effect  as  Finish. 

2.3.0  CROUP  DATA 

Croups  are  emploved  by  CORCEN  as  an  intermediate  organisa¬ 
tional  structure  to  facilitate  element  qeneration.  The  CRonp 
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connand  options,  illustrated  in  Finure  2.24,  allow  the  user  to 
alter  the  node  ranqes  defined  in  the  groups,  to  list  defined 
q roups,  to  select  alternative  qrouns  or  order  of  qrouns  for 
element  qeneration  and  to  define  new  qrouns  for  use  in  element 
neneration  or  new  node  interpolation.  The  concept  of  the  qroup 
consistinc  of  cross-sections  of  sone  model,  nrecludes  the 
definition  of  anv  node  in  more  than  one  qroup  specified  for 
element  qeneration.  The  user  nav,  however,  create  new  croups 
with  existinq  nodes  for  the  purposes  of  interrola tine  new  nodes. 
The  new  qrouns  may  contain  qreater  or  fewer  nodes  than  the 
orininal  croups  and  will  replace  the  oriqinal  croup  sent  to 
the  interpolation. 

CORfFP  provides  interactive  prompts  to  insure  that  the  user 
properlv  defines  new  qroups  or  makes  valid  chances  to  existinc 
croups.  Croups  are  defined  bv  a  beqinninc  node,  an  endinc  node 
and  a  qroup  desiqnator.  The  croup  desiqnator  assists  the  user 
in  keepinq  track  of  how  the  qrouns  were  created.  Table  2.2 
lists  the  different  qroup  designators  that  are  assicned  and 
defines  their  neaninqs.  The  current  options  for  CORCCP  assume 
that  the  user  is  able  to  always  input  or  cenerate  nodes  such 
that  functionally  related  nodes  are  adiacent  to  one  another. 
Special  cases  nay  arise  where  the  user  must  specify  random  lists 
of  nodes  for  interpolation  to  yield  a  new  qroup  with  the  proper 
characteristics  for  the  model  to  be  defined. 

2.3.9  IITTCRPOLATF  PATA 

CORCHN  provides  the  user  v/ith  a  means  of  creatine  new 
coordinates  based  on  the  model  data  provided.  The  iMTrrrolate 
command  allows  users  to  input  either  a  croun  of  nodes  or  a  list 
of  nodes  to  use  for  the  creation  of  new  nodes.  The  input  for 
the  interpolation  procedures  requests  the  user  to  specify  what 
nodes  or  qroup  to  use  and  how  many  nodes  the  new  croup  beinc 
created  should  have.  CORCF’I  will  interpolate  the  q  iven  nodes 
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GROUP  -  select  group  specification  and  editing 


ALTER  -  change  group  composition 

INTERPOLATE  -  generate  nodes  with  selected  group 
KEYBOARD  -  keyboard  input  of  group  composition 
HELP  -  help  user 
LIST  -  list  groups 
DELETE  -  delete  groups 

TABLE  -  give  summary  of  group  statistics 


Figure  2.24.  GROUP  Command  Options. 
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TABLE  2.2 


CORGEN  GROUP  DESIGNATORS 


Groups  will  be  labelled  depending  on  the  type  of  input 
used  to  define  the  nodes  to  the  program.  The  following 
group  designators  are  used: 


A 

C 

D 

F 

I 

K 

L 

S 

T 


-  arc  generator 

-  curve  generator 

-  deleted  group 

-  file  input 

-  interpolated 

-  keyboard  input 

-  line  generator 

-  screen  digitizer 

-  tablet  digitizer 


and  provide  a  new  qroup  of  nodes  if  no  qrouo  was  specified  or 
it  will  ovennrite  the  qroup  that  was  qiven  as  input.  The  two 
end  point  nodes  of  the  qroup  will  be  included  in  the  new  arotm 
to  insure  the  model  boundaries  are  not  affected.  r’hile  the  new 
qroup  will  replace  the  old  qroup,  new  node  nunbers  must  be 
qenerated  for  the  nodes  created  by  interpolation.  These  new 
node  numbers  in  no  wav  correspond  to  the  old  node  nunbers.  nhe 
actual  nodal  coordinates  used  to  do  the  interpola tion  are  in  no 
wav  affected  by  the  qeneration  of  new  nodes. 

2.3.in  none  DATA 

COROEH  is  structured  to  provide  the  user  with  several  node 
input  features  which  operate  independently  and  concurrently . 

The  user  nay  access  a  FI  LF  input  of  nodal  data,  then  enter  the 
PTYHoard  option  to  add  additional  coordinate  data.  In  addition, 
the  dinitizinq  command  nay  be  executed  at  any  tine  to  add  loft inn 
type  data  coordinates.  The  node  input  routines  will  prompt  the 
user  to  declare  how  many  nodes  are  to  be  included  in  the  next 
qroup.  If  appropriate,  the  user  also  will  be  asked  to  specify 
how  many  qroups  are  to  be  created  before  proqran  control  will  be 
returned  to  the  user.  The  input  routine  will  then  create  or 
input  enouqh  nodes  to  equal  the  number  of  nodes  per  qroup  tines 
the  number  of  qroups.  The  qroups  created  nay  then  be  utilized 
to  create  other  nroups  of  nodes  or  to  specify  how  nodes  should 
be  orqanized  for  element  qeneration.  Fioure  2.25  lists  the 
options  for  the  I10PF  command. 

All  nodal  coordinate  data  is  stored  as  four  values.  Ppon 
input  to  CORCinri  all  data  is  converted  to  rectilinear  coordinates 
and  stored  as  (X,Y,2)  triples  for  each  coordinate  location.  In 
addition,  a  thickness  parameter  is  associated  with  each  coor¬ 
dinate  set  for  use  in  FXPAPPinq  the  surface  model  to  a  solid 
model  prior  to  use  by  the  PRFP  preprocessor.  The  user  nav 
specify  a  default  thickness  value  which  will  be  anponded  to  all 


ALTER  -  change  nodal  attributes 

GENERATE  -  create  new  nodes  using  automatic  generators 

FILE  -  input  nodes  from  a  file 

KEYBOARD  -  input  nodes  from  keyboard 

HELP  -  help  user 

LIST  -  list  nodes 

DELETE  -  delete  nodes 

TABLE  -•  provide  a  table  of  nodal  statistics 
INTERPOLATE  -  create  new  nodes  via  interpolation 


Figure  2.25.  NODE  Command  Options. 


the  coordinate  data  input  until  the  user  resets  it.  The  AL^Er 
option  also  permits  the  user  considerable  flexibility  in  nakinq 
chanqes  to  attributes  (e.q.  -  X, Y , Z ,  THICK )  of  a  ranqe,  list  or 
qroun  of  nodes.  Hew  nodes  nay  be  input  to  CORFEH  via  the  ALTFr 
node , . however ,  all  nodes  entered  in  this  fashion  will  be  placed 
in  a  sinqle  new  qroup.  A  new  qroup  will  he  created  each  tine 
the  user  leaves  the  sinqle  node  edit  option  and  one  or  more 
previously  undefined  nodes  were  entered. 

KFYPoard  input  of  nodes  allows  the  user  to  input  nodes  in 
coordinate  triples  plus  a  thickness  value.  Any  coordinate 
transformation  can  be  specified  for  the  data  as  well  as  a 
separate  coordinate  translation.  Loftinq  planes  cannot  be  used 
with  the  keyboard  input.  All  data  for  KEYHoard  input  is  read 
format-free,  meaninq  the  user  may  type  the  data  bv  separatinc 
each  value  with  a  comma  or  blank. 


FILE  input  provides  the  user  with  a  convenient  means  to 
input  coordinate  data  already  available  on  cards,  tapes  or 
disk  files.  The  user  may  input  data  as  coordinate  triples 
plus  a  thickness  value  or  may  select  to  specify  a  default 
thickness  and  just  input  the  coordinate  triples.  Aqain,  the 
coordinate  transformations  may  be  selected  for  spherical  or 
cylindrical  coordinates  if  necessary.  The  user  must  be 
careful  to  specify  how  many  croups  are  to  be  constructed  by 
the  current  REAP  option.  The  command  will  request  th  ?  user 
to  specify  the  number  of  nodes  that  are  to  be  included  in 
each  qroup  be  inq  built  then  the  number  of  croups  to  be 
constructed.  The  number  of  nodes  on  the  file  must  be  at 
least  the  number  of  nodes  per  qroup  times  the  number  of  nroups 
to  build. 

The  Digitize  command  allows  the  user  to  create  dicitizinc 
planes  arbitrarily  located  in  three-dimensional  space.  The 
dicitizinc  planes  are  defined  either  by  nrovidinc  the  equation 


of  the  plane,  three  non-collinear  points  or  three  nodes 
currently  defined  that  lie  v/ithin  the  plane.  Once  the  lofting 
plane  has  been  defined  the  user  nay  beqin  digitizing  coor¬ 
dinates  into  croups,  where  each  loft  plane  will  generally 
correspond  to  a  separate  group.  The  conmand  assunes  a 
Cartesian  coordinate  systen  for  all  dinitizino,  although  the 
user  nay  select  coordinate  translations  to  perfom  on  the 
data  prior  to  the  trans fornation  fron  the  digitizing  plane  to 
global  coordinates.  The  user  may  exit  the  digitizing  node  in 
two  ways,  one  is  to  select  the  nenu  option  and  digitize  the 

EXIT  option  in  the  nenu  area.  The  second  nanner  is  to 
digitize  the  sane  point  three  consecutive  tines.  In  the  latter 
case  the  last  three  points  will  all  be  lost. 


2.3.11  RESTORE  DATA 

Data  that  has  been  previously  generated  by  COROEM,  ACRID  or 
SPATCH  in  the  unformatted,  surface  model  file  structure  nay  be 
submitted  to  the  CORCEN  program  with  the  RESTore  command.  This 
command  will  cause  the  current  model  in  the  CORC.ETT  systen  to  be 
deleted  and  the  new  model  to  replace  it.  Once  the  command  has 
been  issued  the  user  will  be  prompted  to  input  the  appropriate 
data  file  specifier  (see  Appendix  C  for  valid  file  names).  CORCEIT 
will  then  input  the  new  nodes  and  elements  and  create  groups  for 
user  convenience. 


2.3.12  SAVE  DATA 

Data  input  to  CORCEN  or  created  by  COROCn  nay  be  saved  at  any 
tine  by  the  user  with  the  SAVE  connand.  This  command  will  cause 
all  currently  defined  nodal  coordinate  data,  element  connectivity 
data  and  groups  specification  data  to  be  written  to  an  unformatted 
surface  model  data  file. 
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2.3.13 


STOP  PROGRAM 


The  STOP  command  will  terminate  CORGEN  program  execution. 
Program  processing  flags  will  be  checked  to  insure  the  user  has 
not  overlooked  any  of  the  crucial  operations  necessary  for  model 
generation.  If  any  phases  of  model  development  have  not  been 
accomplished  the  program  will  perform  any  operations  that  are 
necessary  to  produce  a  complete  model  file. 
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2.4  AGRID  -  Arbitrary  Arrays  of  Points  Mesh  Input 


The  AGRID  data  entry  program  is  designed  for  use  with  data 
that  is  unsuitable  for  use  with  the  remaining  data  input  modules, 
usually  due  to  a  lack  of  order  within  the  geometric  data  itself. 
AGRID  can  accept  a  virtually  arbitrary  collection  of  point  coor¬ 
dinates  which  lie  on  a  surface  of  interest,  and  generate  a  regular 
grid  of  mesh  points  which  define  a  smooth  surface  passing  through 
the  given  data. 


The  method  of  surface- fitting  used  within  AGRID  is  based  upon 
a  minimum-curvature  criterion,  by  which  the  given  data  (in  the 
form  of  point  coordinates  X,  Y ,  Z)  is  used  to  define  a  continuous 
surface  Z  =  Z(X,Y)  which  exhibits  a  minimal  amount  of  oscillation 
between  the  known  points.  In  particular,  the  final  surface  is 
represented  by  a  grid  of  piecewise  bicubic  functions  whose  coef¬ 
ficients  are  determined  by  minimizing  the  curvature  functional 


+  Z 


ryy 


+  (z  z  -z2  )1 

,xx  ,yy  ,xy  J 


dx  dy 


(2.2) 


which  reflects  both  the  mean  and  Gaussian  curvatures  of  the 
surface.  The  minimization  is  carried  out  by  solving  a  small 
finite  element  problem,  in  which  each  bicubic  segment  of  the 
surface  corresponds  to  a  single  finite  element. 


Figures  2.26  through  2.28  illustrate  the  AGRID  solution  pro¬ 
cedure.  The  known  points  lying  on  the  surface  to  be  defined  are 
enclosed  within  a  rectangular  region  which  will  form  a  temporary 
grid  for  interpolation  (Figure  2.26).  A  mesh  of  'interpolation 
elements'  is  defined  by  specifying  the  series  of  X  (or  Y)  coor¬ 
dinates  at  which  mesh  points  are  to  be  positioned,  as  shown  in 
Figure  2.27.  The  mesh  size  and  refinement  are  arbitrary,  but 
more  divisions  would  normally  be  used  in  those  regions  containing 


the  largest  number  of  prescribed  coordinates  or  exhibiting  the 
most  changes  in  curvature.  The  coordinates  of  known  points  on 
the  surface  are  then  specified;  these  may  be  located  anywhere 
within  the  interpolation  grid,  but  not  outside  it.  Finally,  a 
set  of  boundaries  is  defined  which  outline  the  true  limits  of 
the  surface  mesh  to  be  generated,  and  the  number  of  surface 
elements  to  be  generated  along  each  direction  of  the  final  mesh 
is  specified.  The  projection  of  these  boundary  curves  on  the 
(X,Y)  plane  is  a  straight-sided  quadrilateral,  as  shown  in 
Figure  2.28.  The  result  of  this  procedure  is  a  mesh  of  surface 
elements  (see,  for  example,  Figure  2.2)  and  coordinates,  suitable 
for  input  to  EXPAND. 

The  necessary  input  to  AGRID  is  summarized  in  Table  2.3. 

This  data  is  normally  saved  on  a  file  (local  file  INGEOM ,  on  CDC 
machines)  arranged  in  free  format,  with  spaces  or  commas 
delimiting  separate  items  of  data  within  a  record.  The  output 
from  AGRID,  in  addition  to  a  minor  amount  of  printed  output,  is 
in  the  form  of  the  shell  surface  geometry  file  described  in 
Subsection  6.1  and  in  Appendix  A. 

It  should  be  observed  that  AGRID  is  based  upon  an  empirical 
means  of  interpolation,  due  to  the  generality  of  the  data  which 
can  be  accepted.  In  return  for  this  generality,  a  compromise  in 
the  quality  of  the  interpolation  can  result  if  the  amount  of  data 
given  is  insufficient  to  define  the  surface  with  a  fair  degree  of 
accuracy.  For  example,  the  specification  of  three  data  points  to 
AGRID  will  always  result  in  the  generation  of  a  planar  surface, 
even  if  the  intended  result  is  a  curvilinear  surface;  four  points 
will  generate  a  warped  surface  with  small  mean  curvatures.  When 
the  local  changes  in  curvature  are  significant,  sufficient  data 
must  be  specified  to  make  them  evident  to  the  generator.  If  this 
advice  is  followed,  AGRID  is  capable  of  producing  very  nice 
results  from  the  most  disorganized  data. 


TABLE  2.3  -  AGRID  Data  Input 


RECORD  VARIABLE  DESCRIPTION 


1 


2 

3 

4 


NXINT  Number  of  X-stations  in  rectangular  interpolation 

grid 

NYINT  Number  of  Y-stations  in  rectangular  interpolation 

grid 

NSPEC  Number  of  prescribed  coordinate  values  to  be  given 


X(NXINT)  Array  of  X-station  values,  in  ascending  order 


Y(NYINT)  Array  of  Y-station  values,  in  ascending  order 


XVAL  X-Coordinate  for  prescribed  point 

YVAL  Y-Coordinate  for  prescribed  point 

I TYPE  Type  of  value  prescribed 

(=1  for  Z,  =2  for  dZ/dX,  =3  for  dZ/dY) 

VALUE  Value  of  Z,  dZ/dX,  or  dZ/dY  at  point  (XVAL, YVAL) 


<<<  Repeat  record  4  for  each  prescribed  point  (NSPEC  times)  >>> 

5  XOUT ( 4 )  X-Coordinates  of  corner  points  of  output  grid 

6  Y0UT(4)  Y-Coordinates  of  corner  points  of  output  grid 


7  NEL1 

NEL2 


Number  of  elements  in  output  grid,  along  edge  1-2 
Number  of  elements  in  output  grid,  along  edge  2-3 


8 


THICK 


Uniform  thickness  value  for  use  in  EXPAND 


2.5  SPATCH  -  Surface  Patch  Geometry  Input 


Recently,  many  interactive  surface  design  systems  have 
adopted  the  bicubic  patch  as  a  basis  for  geometric  description. 
The  SPATCH  utility  is  a  special-purpose  conversion  program  which, 
with  the  aid  of  a  few  simple  user-written  subroutines,  allows 
the  translation  of  bicubic  surface  patch  data  into  the 
preprocessor  data  format. 

With  the  type  of  surface  patch  considered  in  SPATCH,  all 
three  Cartesian  coordinates  X,  Y,  Z  are  described  as  parametric 
functions  of  two  natural  coordinates  (u,v), 

X  =  X  ( u , v) 

Y  -  Y  ( u , v )  (2.2) 

Z  *  Z  (  u ,  v ) 

where  u  and  v  each  vary  between  zero  and  one.  The  most  common 
form  of  storing  the  surface  patch  data  is  in  arrays  having  the 
form  17] 


xoo 

X01 

x 

vOO 

X 

vOl 

x,  „ 

X,  , 

x 

x 

10 

11 

vio 

vll 

X 

uOO 

XuCl 

X 

uvOO 

Xuv0l 

X  .  . 
ulO 

xull 

XuvlO 

Xuvll 

for  each  of  the  coordinates  X,  Y,  Z.  The  indices  (o,l)  refer  to 
the  values  of  u  and  v  respectively,  and  subscripts  u,  v  denote 
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parametric  derivatives  of  the  coordinate  in  question.  For  - 

example,  Xvlo  is  the  value  of  dx/dv  at  u=l  and  v=o. 

SPATCH  permits  surfaces  described  in  the  form  of  a  collection 
of  such  patches  to  be  translated  into  an  assemblage  of  biquadratic 
surface  elements  (see  Subsection  2.1,  Figure  2.2)  suitable  for  input 
to  EXPAND.  The  input  to  SPATCH  is  performed  through  a  user-written 
subroutine,  since  the  actual  storage  of  the  surface  patch  data  may 
taken  on  a  number  of  different  forms.  The  user  routine  has  the  form 

SUBROUTINE  UPATCH  (PATCH,  THICK,  NU,  NV,  ICOUNT,  IEND) 

DIMENSION  PATC H ( 4 , 4 , 3 ) 


<  code  to  read  surface  patch  data  > 


RETURN 


The  exact  meanings  of  all  formal  parameters  and  other  special 
requirements  are  discussed  in  Appendix  E.  Each  surface  patch 
(described  by  the  coefficients  in  array  PATCH)  may  be  divided 
into  one  or  more  nine-node  surface  elements,  as  determined  by 
the  parameters  NU  and  NV. 


Two  additional  user-written  subroutines  are  provided  for  in 
SPATCH,  to  permit  opening  and/or  positioning  of  the  patch  data 
file  prior  to  entering  UPATCH,  and  closing  of  the  file  upon 
completion.  These  two  routines  will  generally  be  necessary  only 
in  the  VAX  version  of  SPATCH,  or  in  the  event  that  other,  extra¬ 
neous  data  precedes  the  surface  patch  data  on  the  file  to  be 
accessed.  These  two  file  control  subroutines,  UOPEN  and  UCLOSE, 
are  also  described  in  Appendix  E. 
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SECTION  3 


PREPROCESSOR  INTERFACE  MODULES 


The  translation  of  existing  finite  element  model  data  into 
the  standard  preprocessor  format  is  accomplished  using  special- 
purpose  data  translation  utilities,  each  one  corresponding  to  a 
specific  pair  of  data  formats.  Conversion  of  modelling  data 
into  the  standard  PREP  format  is  currently  possible  for  two  data 
file  types:  * 

-  MAGNA  [4]  finita  element  data  deck,*  and 

-  IMPRESS  [3]  preprocessor  database. 

These  two  data  conversion  utilities  are  described  in  subsections 
3.1  and  3.2,  respectively. 

In  addition  to  the  special-purpose  translation  programs 
mentioned  above,  a  set  of  general-purpose  utility  subroutines 
exist  for  the  purpose  of  converting  finite  element  data  in  other 
external  formats.  These  general  interface  routines  are  the 
subject  of  subsection  3.3. 


♦Conversions  to  and  from  the  PREP  file  format  can  also  be 
performed  using  the  NEUTRAL  interface  program  (Section  5) . 
NEUTRAL  is  used  primarily  for  archival  of  data  in 
formatted  form. 


3.1  TRNSFR  -  MAGNA  INPUT  DATA  TRANSLATOR 

The  TRNSFR  data  translator  accepts  as  input  the  standard  data 
deck  recognized  by  the  MAGNA  finite- element  program  [4].  Only  the 
nodal  coordinates  and  element  connectivity  data,  which  represents 
the  bulk  of  the  model  information,  are  converted  to  the  pre¬ 
processor  internal  format.  All  of  the  standard  MAGNA  element 
types  are  supported  by  TRNSFR;  these  include: 

-  truss  (bar)  elements,  with  two  nodes  per  element, 

-  beam  elements,  with  two  or  three  nodes  per  element, 

“  plan©  stress,  plane  strain  and  shear  panel  elements 

having  from  four  to  nine  nodes  per  element, 

-  axisymmetric  elements  with  four  to  nine  nodes, 

-  the  MAGNA  thin  shell  elements,  with  eight  or  sixteen 
vertex  nodes,  and 

-  three-dimensional  solids,  with  from  eight  to  27  nodes. 

TRNSFR  is  capable  of  accepting  all  of  the  incremental  node  point 
and  element  generating  sequences  possible  with  the  MAGNA  program. 

Execution  of  TRNSFR  is  quite  simple  (see  Appendix  C).  On  CDC 
computers,  the  MAGNA  input  deck  must  be  attached  as  a  local  file 
prior  to  execution,  while  on  the  VAX  11/780,  the  data  file  name  is 
prompted  by  the  program.  No  further  input  or  action  is  required 
of  the  user ,  other  than  saving  the  generated  preprocessor  data 
file. 


3.2  IMPRINT  -  IMPRESS  PREPROCESSOR  INTERFACE 


IMPRESS  is  the  main  preprocessing  program  of  a  finite  element 
data  generation  system  consisting  of  three  distinct  programs 
(IMPRESS,  INTRPOS ,  and  IMPACT),  all  of  which  are  described  in 
Reference  [3].  The  modelling  features  available  in  IMPRESS 
include  a  flexible  area-outlining  approach  to  geometry  definition, 
and  the  automatic  calculation  of  intersection  curves  for  planar, 
cylindrical  and  spherical  shapes.  Both  two-  and  three-dimensional 
geometric  models  may  be  prepared  with  the  program. 

Finite  element  data  originally  generated  using  IMPRESS  can  be 
converted  for  use  with  the  present  preprocessor  by  means  of  the 
IMPRINT  translator.  IMPRINT  is  a  collection  of  conversion 
routines  which  can  be  used  as  the  output  overlay  of  the  IMPRESS 
data  base  access  program  INTRPOS.  The  IMPRINT  overlay  is  merged 
with  INTRPOS  using  the  control  statement  sequence  described  in 
Appendix  I  of  Reference  C 3 3 .  When  INTRPOS  is  executed,  the  data 
conversion  procedure  is  invoked  by  entering 

PREPARE  RAW 

Upon  completion  of  the  data  conversion,  the  END  directive  can  be 
used  to  exit  INTRPOS.  The  converted  finite  element  data  is 
written  to  the  local  data  file  TAPE9,  which  should  be  saved 
following  execution  of  the  program. 


Table  3.1  summarizes  the  IMPRESS  data  types,  and  indicates 
those  which  are  supported  by  the  IMPRINT  utility. 


TABLE  3.1  -  IMPRESS  Data  Types 


IMPRESS  Data  Type 


Supported  by  IMPRINT 


Nodal  Coordinates 
Skewed  Coordinate  Systems 
Boundary  Conditions 
Master  /  Slave  Constraints 
Prescribed  Displacements 
Element  Connectivity 
Element  Property  Codes 
Beam  Section  Properties 
Concentrated  Springs  /  Masses 
Nodal  Concentrated  Forces 
Multiple  Loading  Sets 


YES 

NO 

YES 

NO 

NO 

YES 

YES 

NO 

NO 

YES 

YES 


3.3  GENERAL-PURPOSE  INTERFACE  ROUTINES 

A  collection  of  general-purpose  translation  routines  has  been 
developed  to  facilitate  the  conversion  of  foreign  data  formats  to 
the  standard  internal  format  of  the  present  preprocessor.  To  use 
these  interface  routines,  the  user  must  write  a  small  driver 
program  which  reads  the  data  to  be  converted,  and  sends  it  to  the 
translation  subroutines  for  consistency  checking  and  output. 

Seven  utility  routines  are  available: 

1.  INICNV  -  Initialization  routine 

2.  NODCNV  -  Nodal  data  conversion 

3.  ELECNV  -  Element  data  conversion 

4.  BCDCNV  -  Boundary  condition  data  conversion 

5.  NLDCNV  -  Nodal  loads  data  conversion 

6.  ELDCNV  -  Element  distributed  loads  data  conversion 

7.  TRMCNV  -  Termination  routine. 

These  interface  routines  are  described  in  detail  in  Tables  3.2 
through  3.8. 

It  is  important  to  note  that  logical  unit  50  must  be  reserved 
for  output  by  the  interface  routines,  and  should  not  be  modified 
by  the  calling  program.  On  CDC  machines,  the  file  TAPE50  must  be 
declared  on  the  PROGRAM  card,  with  a  minimum  buffer  length  of  512 
words.  For  the  VAX  version  of  the  conversion  routines,  the  output 


TABLE  3 . 2 


INICNV  -  Initialization  Subroutine 


PURPOSE  -  Define  the  number  of  data  entries  to  be  made  in  the 

file,  and  initialize  internal  parameter  values. 

ACCESS  -  CALL  INICNV  ( NUMNP , NUMEL , NUMBC , NUMNL , NUMEF ) 

PARAMETERS  - 

NUMNP  (input)  -  Number  of  nodal  points  in  the  model 
NUMEL  (input)  -  Number  of  finite  elements 
NUMBC  (input)  -  Number  of  boundary  condition 

specifications  - 

NUMNL  (input)  -  Number  of  nodal  loads  specifications 
NUMEF  (input)  -  Number  of  element  loads  specifications 

NOTES  -  -  (1)  This  subroutine  must  be  called  before  any  of  the 
other  data  translation  routines  are  called. 

(2)  Nodal  points  must  be  numbered  from  1  through  NUMNP. 

(3)  Each  boundary  condition  specification  may  consist  of 
a  range  of  nodes  (first,  last,  increment)  all  having 
the  same  constraints  applied. 

(4)  Each  nodal  load  specification  may  consist  of  a  range 
of  nodes  ( first,  last,  increment)  and  three 
components  of  force  (X,  Y,  Z)  which  apply  to  all 
nodes  in  the  sequence. 

(5)  Each  element  load  specification  may  consist  of  a 
range  of  elements  ( first,  last,  increment)  and  a 
loading  type  and  value  which  apply  to  all  elements 
in  the  range. 


TABLE  3 . 3 


NODCNV  -  Nodal  Coordinate  Conversion  Subroutine 


PURPOSE  -  Process  coordinate  data  for  a  single  node  point  of  the 
model . 

ACCESS  -  CALL  NODCNV  (NODE,X,Y,Z) 

PARAMETERS  -  .... 

(input)  -  Node  point  sequence  number 
(input)  -  Cartesian  X-Coordinate 
(input)  -  Cartesian  Y-Coordinate 
(input)  -  Cartesian  Z-Coordinate 

-  (1)  Node  points  must  be  transmitted  to  NODCNV  sequen¬ 
tially,  with  all  node  numbers  between  1  and  NUMNP. 
(2)  Node  coordinates  must  be  defined  prior  to  defining 
elements,  boundary  conditions  or  loading  data. 


NOTES 


TABLE  3.  4 


ELECNV  -  Element  Data  Conversion  Subroutine 


PURPOSE  -  Processes  a  single  finite  element  of  the  model. 

ACCESS  -  CALL  ELECNV  ( IELNUM , MATLNO , I  DIMEN, MAXNOD, NCON) 

PARAMETERS  - 

IELNUM  (input)  -  Element  number  (arbitrary) 

MATLNO  (input)  -  Material  number  (optional) 

IDIMEN  (input)  -  Dimensionality  of  the  element  (1,2,3) 
MAXNOD  (input)  -  Maximum  local  node  number 
NCON  (input)  -  List  of  connected  nodes 

NOTES  -  (1)  If  the  material  number  is  to  be  ignored,  set  MATLNO=0 

(2)  MAXNOD  should  lie  within  the  following  ranges: 

MAXNOD  =*2  -  One-dimensional  elements; 

MAXNOD  =4-9  -  Two-dimensional  elements;  and 

MAXNOD  =8-27  -  Three-dimensional  elements. 

(3)  The  length  of  the  connection  array  NCON  is  determined 
by  MAXNOD.  If  variable-number-of-nodes  elements  are 
used,  intermediate  entries  for  which  no  node  exists 
must  be  entered  as  zeroes. 

(4)  The  conventions  for  node  numbering  in  all  element 
types  are  shown  in  Figures  2.1  and  2.2. 


TABLE  3.5 


BCDCNV  -  Boundary  Condition  Conversion  Subroutine 


PURPOSE  -  Processes  a  single  boundary  condition  specification. 

ACCESS  -  CALL  BCDCNV  (ISTART, IEND, INCR,  IX, IY,  IZ) 

PARAMETERS  - 

ISTART  (input)  -  First  node  in  the  sequence 
IEND  (input)  -  Last  node  in  the  sequence 
INCR  (input)  -  Node  number  increment 

IX  (input)  -  X-Direction  constraint  code  (0,1) 

IY  (input)  -  Y-Direction  constraint  code  (0,1) 

IZ  (input)  -  Z-Direction  constraint  code  (0,1) 

NOTES  -  (1)  The  specified  constraints  will  be  applied  at  nodes 

ISTART,  (ISTART+INCR) ,  ( ISTART+2 *INCR) ,  ...  ,  IEND. 

If  a  single  node  is  to  be  constrained,  IEND  and  INCR 
may  be  set  to  zero;  however,  they  must  appear  in  the 
calling  statement. 

(2)  IX,  IY,  and  IZ  are  defined  as  zero  if  motion  is 

permitted  in  the  X,  Y,  or  Z  directions,  respectively 
Constraints  are  specified  by  IX,  IY,  or  IZ  =  1  as 
appropriate . 


TABLE  3.6 


NLDCNV  -  Nodal  Loads  Conversion  Subroutine 


PURPOSE  -  Processes  a  single  nodal  loading  specification. 

ACCESS  -  CALL  NLDCNV  (ISTART, IEND, INCR, FX, FY, FZ, ICASE) 

PARAMETERS  - 

ISTART  (input)  -  First  node  number  in  sequence 
IEND  (input)  -  Last  node  number  in  sequence 
INCR  (input)  -  Node  number  increment 

FX  (input)  -  X-Direction  force  component 

FY  (input)  -  Y-Direction  force  component 

FZ  (input)  -  Z-Direction  force  component 

ICASE  (input)  -  Load  case  or  group  number 

NOTES  -  (1)  The  range  of  nodes  to  receive  the  specified  load  is 
defined  exectly  as  in  the  boundary  condition  data; 
when  only  a  single  node  is  to  be  loaded,  IEND  and 
INCR  may  be  set  to  zero. 

(2)  Any  integer  value  may  be  ured  for  the  load  case 
number. 


TABLE  3.7 


ELDCNV  -  Element  Loads  Conversion  Subroutine 


PURPOSE  -  Processes  a  single  distributed  element  load  specifi¬ 
cation. 


ACCESS  -  CALL  ELDCNV  ( I  START, I  END , INCR, ITYPE , ICASE, FORCE) 
PARAMETERS  - 


NOTES 


I START 

I  END 

INCR 

ITYPE 

ICASE 

FORCE 


(input) 
( input ) 
( input) 
( input) 
( input) 
( input) 


-  First  element  number  in  sequence 

-  Last  element  number  in  sequence 

-  Element  number  increment 

-  Loading  type/-  direction/  surface  code 

-  Load  case  or  group  number 

-  Loading  magnitude 


(1)  Element  ranges  are  defined  similarly  to  the  nodal 
ranges  for  boundary  conditions  and  nodal  loads.  If 
a  single  element  is  to  be  loaded,  I END  and  INCR  may 
be  set  to  zero,  but  they  must  appear  in  the  calling 
statement. 

(2)  ITYPE  is  defined  as  follows: 


ITYPE  »  -1,  -2,  -3  signifies  a  body  force  (force 
per  unit  volume)  in  the  X,  Y,  or  Z 
direction,  respectively. 

ITYPE  =1,  2,  ...  ,  6  denotes  a  surface  pressure 

on  element  faces  1,  2,  ...,  6,  respectively, 

for  three-dimensional  elements  (IDIMEN=3). 

ITYPE  =1,  2,  3,  4  represents  a  line  load  (force 

per  unit  length)  on  the  corresponding  edge 
of  a  two-dimensional  element  (IDIMEN=2). 


TABLE  3.8 


TRMCNV  -  Termination  Subroutine 

PURPOSE  -  Checks  for  incomplete  blocks  of  data,  and  closes  the 
generated  data  file. 

ACCESS  -  CALL  TRMCNV 

PARAMETERS  -  (none) 

NOTES  -  (1)  In  the  CDC  version  of  the  .termination  routine,  the 
output  data  file  is  rewound,  but  not  returned. 


file  generation  is  user-transparent,  with  the  converted  data  being 
written  automatically  to  the  file  UNFMT.DAT.  A  labelled  COMMON 
block  (COMMON  /DATCNV / )  is  also  used  in  the  interface  subroutines, 
and  should  not  be  modified  by  the  user. 
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MODEL  VERIFICATION  AND  MODIFICATION 


Capabilities  for  developing  a  completed  finite  element  model 
from  an  initial  geometric  description  are  contained  in  the  central 
program,  PREP.  Typical  operations  performed  within  PREP  are  mesh 
refinement,  assignment  of  properties,  constraints  and  loads,  and 
geometry  plotting.  In  the  following  subsections,  each  modelling 
operation  available  in  the  PREP  preprocessor  is  described,  and 
typical  sequences  of  operations  are  described  which  make  the  most 
effective  use  of  the  program's  capabilities. 

The  general  mode  of  operation  of  the  PREP  module  involves  the 
use  of  model  data  files,  with  up  to  13  files  being  active  at  any 
one  time.  Each  file  may  be  named  for  later  reference.  The 
typical  PREP  operation  then  consists  of  the  sequence 

INPUT  MODEL - >  (OPERATION) - >  OUTPUT  MODEL 

FILE ( S )  FILE ( S ) 

For  example,  the  user  may  wish  to  refine  the  finite  element  mesh 
in  a  certain  area  of  a  model.  Upon  entering  the  command  'REFINE', 
he  is  requested  to  enter  the  name  of  an  existing  model  whose 
element  mesh  is  to  be  refined.  PREP  then  requests  additional 
input  to  define  those  elements  to  be  subdivided  and  the  degree 
of  refinement,  and  generates  a  new  model  file  (which  is  assigned 
a  new  name)  with  the  requested  element  refinement. 

The  input/output  files  for  the  PREP  preprocessor  are  described 
in  Section  6.1  and  in  Appendix  A.  PREP  is  organized  such  that 


data  files  may  be  manipulated  in  the  program,  saved  on  peripheral 
storage,  and  reentered  into  PREP  at  a  later  time.  This  capability 
permits  the  generation  of  a  model  to  be  accomplished  in  several 
sessions  if  necessary,  and  allows  separate  segments  of  the  finite 
element  model  to  be  prepared  separately  and  combined  when 
compl e  te . 

The  various  modelling  functions  available  in  PREP  are 
described  in  Subsections  4.1  through  4.7;  those  utilities  which 
are  useful  in  managing  data  files  within  the  program  are  also 
discussed  in  Subsection  4.1.  Subsection  4.8  outlines  some 
suggested  sequences  of  PREP  operations  which  have  proved  to  be 
most  effective  in  practice.  Subsection  4.8  also  contains  an 
alphabetized  listing  of  all  PREP  commands  for  quick  reference. 
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4.1  MODIFICATION  /  REFINEMENT  OF  MODEL  GEOMETRY 

PREP  operations  which  are  provided  for  manipulation  of  the 
geometry  of  a  finite  element  model  are: 

-  CREATE  -  build  model  from  'scratch'  by  user  input; 

-  EDIT  -  alter  coordinates  and  connectivity; 

-  FILL  -  nenerate  midside  and  interior  node  points; 

-  MAS?'  -  renove  midside  and  interior  node  points; 

-  MERGE  -  combine  two  models  with  compatible  geometries; 

-  REFINE  -  subdivide  finite  elements  to  form  a  finer  mesh; 

-  REFLECT  -  form  the  mirror  inane  of  a  model  about  a  plane; 

-  ROTATE  -  rotate  an  entire  model  in  three  dimensions;  and 

-  TRANSLATE  -  translate  a  model  in  three  dimensions. 

The  CREATE  and  EDIT  commands  permit  input,  modification  or  dele¬ 
tion  of  individual  items  of  modeling  data.  The  remaining  command 
are  geometric  operations,  which  alter  either  the  coordinate  data 
or  the  element  connection  data  for  the  entire  model  in  question. 
The  effect  of  each  of  these  commands,  as  well  as  the  available 
options  for  data  input,  are  summarized  in  Tables  4.1  throunh  4.7. 
An  example  of  each  operation  is  given  in  the  accompanying  Finures 
(4.1  throunh  4.4). 

It  is  qenerally  advisable  to  perform  FILL  and/or  MASK  opera¬ 
tions  at  the  beginning  of  a  modelling  session,  particularly  in 
cases  where  an  irregular  number  of  nodes  has  been  used  in  the 
initial  generation  of  a  model.  An  example  in  which  REFINE  has 
been  used  prior  to  FILL  with  such  an  irrenular  mesh,  is  shown  in 
Figure  4.5;  here  the  initial  geometry  has  been  defined  using 
elements  with  only  one  midside  node,  and  the  effect  of  the  REFINE 
operation  is  disastrous.  The  remainino  functions  (MERGE,  REFINE, 
REFLECT,  ROTATE  and  TRANSLATE)  can  normally  be  invoked  in  any 
order  as  circumstances  dictate. 


TABLE  4.1  -  The  FILL  Command 


Command  :  FILL 

Effect  :  Generates  midside  and/or  interior  nodes  to 

create  higher-order  finite  elements. 

Number  of  Input  Files  :  1 

Number  of  Output  Files:  1 

Options  :  (1)  All  3-D  elements  FILLed  to  16,  20,  26  or  27  nodes 

(2)  All  2-D  elements  FILLed  to  8  or  9  nodes 

Related  commands  /  options  :  MASK  is  an  ' inverse'  operation  to  FILL 

:  This  function  does  not  preserve  boundary  condition  and 
loading  data. 


Notes 


TABLE  4.2 


The  MASK  Command 


Command  :  MASK 


Effect  :  Removes  midside  and  /  or  interior  nodes  from  elements. 
Number  of  Input  Files  :  1 

Number  of  Output  Files  :  1 

Options  :  (1)  All  3-D  elements  MASKed  to  8,  16,  20,  26  or  27  nodes 
(2)  All  2-D  elements  MASKed  to  4,  8  or  9  nodes 


Related  commands  /  options  s  FILL  is  an  ’inverse'  function  to  MASK 


(1)  FILL  does  not  preserve  constraint  or  loading  data. 

(2)  Care  should  be  taken  in  using  MASK  with  curved 
elements,  since  straight-sided  elements  will  result 
in  many  cases. 
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TABLE  4.3  -  The  MERGE . Command 


Command  :  MERGe 

Effect  :  Combines  two  models  to  form  a  single  (output)  model,  with 
coincident  nodes  equivalenced  and  resequenced. 

Number  of  Input  Files  :  2 

Number  of  Output  Files  :  1 

Options  :  (None) 

Related  commands  /  options  s  The  TOLErance  command  can  be  used  to 

modify  the  distance  tolerance  used  in 
detecting  coincident  nodes. 

Notes  :  MERGE  does  not  preserve  constraint  or  loading  data. 


TABLE  4.4  -  The  REFINE  Command 


Command  :  REFIne 

Effect  :  Subdivides  selected  elements  of  a  model  to  produce  a 
finer  mesh  for  analysis. 

Number  of  Input  Files  :  1 

Number  of  Output  Files  :  1 

Options  :  (1)  Up  to  150  elements  may  be  refined  at  one  time. 

The  selection  options  include 

-  random  element  numbers 

-  ranges  of  elements  ( first,  last,  increment) . 

(2)  Elements  must  be  refined  in  one  direction  at  a 
time.  The  refinement  direction  follows  the 
natural  coordinate  directions  of  an  element 
(see  Figure  4.3).  The  PLOT  option  ORIENT  should 
be  used  to  determine  the  proper  direction  for 
element  refinement. 

(3)  From  one  to  five  'cuts'  may  be  made  through  each 
element  specified,  resulting  in  from  two  to  six 
elements  as  output.  The  default  option  is  equal 
spacing  of  these  cuts?  by  overriding  the  default, 
cuts  may  be  spaced  at  irregular  intervals  across 
the  elements. 

Related  commands  /  options  :  The  PLOT  utility  can  be  used  to 

identify  the  correct  element  numbers  for  input  to 
REFINE  (select  the  LABEL  option),  as  well  as  the 
correct  local  directions  (use  the  ORIENT  option). 

Notes  :  REFINE  does  not  preserve  boundary  condition  or  loading  data. 


TABLE  4.5 


The  REFLECT  Command 


Command  :  REFLect 

Effect  :  Forms  the  mirror  image  of  a  model  with  respect  to  a 
specified  plane  in  three-dimensional  space. 

Number  of  Input  Files  :  1 

Number  of  Output  Files  :  1 

Options  :  The  plane  of  reflection  may  be  defined  in  two  ways: 

-  Give  the  coefficients  of  Ax  +  By  +  Cz  =  D?  for 
example,  the  plane  x  =  5  has  A  =  1 ,  B  =  0,  C  =  0, 

D  =  5 . 

-  Define  three  non-collinear  points  which  lie  in  the 
plane. 

Related  commands  /  options  :  REFLECT  is  often  used  in  conjunction 

with  MERGE  to  create  a  full  model 
from  a  symmetric  one. 


Notes  :  (None) 


■■  TV 


TABLE  4.6 


The  ROTATE  Command 


Command  :  ROTAte 


Effect  :  Performs  a  sequence  of  three  rotations  of  a  model  in 
space. 

Number  of  Input  Files  :  1 

Number  of  Output  Files  :  1 

Options  :  (None) 

Related  commands  /  options  :  MERGE  may  be  used  following  ROTATE 
to  generate  a  full  model  from  one  for  which  a  single 
sector  (e.g.,  in  polar  coordinates)  has  been  defined 

s  Model  rotations  are  always  specified  in  degrees,  and 
are  always  performed  in  the  orders  X-rotation,  then 
Y- rotation,  then  Z- rotation. 


Notes 


TABLE  4.7 


The  TRANSLATE  Command 


Command  :  TRANs late 

Effect  :  Performs  a  rigid-body  translation  of  a  model  in  space 
Number  of  Input  Files  :  1 

Number  of  Output  Files  :  1 

Options  :  (None) 

Related  commands  /  options  :  MERGE  is  useful  in  conjunction  with 
the  TRANSLATE  function  when  the  geometry  to  be 
generated  is  periodic. 


Notes 


:  (None) 


In  addition  to  the  above  qeonetric  nodellinq  functions,  two 
inportant  auxiliary  operations  are  included  in  PREP  to  ensure  the 
correct  qeneration  of  finite  elenent  data.  They  are 

-  SIFT  -  renove  unused  nodes  fron  the  nodel ,  and 

-  TIDY  -  conbine  all  coincident  nodes  into  sinqle  node  points 

with  a  unique  node  nunber. 

For  models  qenerated  usinq  the  present  preprocessor,  these  two 
utilities  should  be  transparent:  SIFT  is  called  automatically 
followinq  any  MASK  operation,  and  TIDY  is  performed  whenever  the 
FILL,  PERCE,  or  REF  IDF  commands  are  executed.  However,  qeonetric 
data  fron  other  sources  may  contain  unused  and  /  or  duplicate 
nodes,  in  which  case  SIFT  and  TIDY  may  be  executed  as  separate 
operations.  Two  options  are  available  in  the  TIDY  operation: 

0PT1  will  qenerate  correct  results  for  any  situation,  but  always 
requires  more  execution  tine;  C'>T2  uses  a  qeonetric  partitioninq 
of  the  nodel  to  reduce  the  complexity  (and  execution  tine)  of  the 
operation,  but  nay  fail  in  some  instances  due  to  storaqe  limita¬ 
tions.  If  the  faster  TIDY  option  fails,  the  0PT1  version  is  in¬ 
voked  and  the  user  is  informed  of  the  situation. 

It  should  be  noted  that  any  qeonetric  model  may  be  listed 
and  /  or  plotted  at  any  tine  durinq  the  session  usinq  the  ”RI"T 
and  PLOT  commands  (Subsections  4.6  and  4.7).  The  PLOT  opti  ' 
•-articularlv  useful  before  a  refine  operation,  to  i-ien  tif”  t'.e 
local  directions  (R,F,T)  within  an  elenent  which  are  needed  as 
data  for  REF  IMF  (see  the  ORIENT  command,  Subsection  4.6). 

Additional  auxiliary  operations  which  are  useful  in  mananinn 
nodel  data  files  durinq  REFINEnent,  MERCinq  and  other  operations 
include 

-  COPY  -  creates  a  copy  of  a  model  file,  under  a  different  nane 

-  delete  -  deletes  the  nane  of  a  model  fron  the  active  list; 


-  HELn  -  nives  a  summary  of  all  available  connands,  and  further 

information  on  selected  commands  if  requested; 

-  NAME  -  assinns  a  user-defined  name  to  a  model  file; 

-  LIST  -  lists  the  nanes  of  all  active  model  files; 

-  TIME  -  prints  the  CPH  time  elapsed  since  siqn-on;  and 

-  TOLERANCE  -  allous  redefinition  of  the  tolerances  used  in 

detectinq  coincident  nodes  in  FILL,  NERCE  and  RTF  I ME . 

The  HELP  and  LIST  commands  are  particularly  useful  as  reminders  of 
operations  which  are  available  and  of  model  filer,  already  created 
durinq  the  session. 

Features  are  provided  in  PREP  to  allow  the  creation  of  models 
from  'scratch'.  The  CREATE  command  provides  a  means  for  input  of 
nodal  coordinates  and  elenent  connectivitv  to  define  a  model.  mhe 
model  must  be  composed  of  3-D  solid  elements  containinq  at  least 
13  nodes  per  element.  The  EnIT  feature  allows  the  user  to  alter 
nodal  coordinates  or  connectivity  for  any  model  in  PEEP.  "hose 
commands  are  fullv  interactive  and  reouire  all  data  to  be  entered 
in  3-P  coordinates. 


4.2  ELEMENT  PROPERTIES  SELECTION 


Three  facilities  exist  in  PREP  for  the  specification  of  finite 
element  properties: 


-  PROPerty  -  define  element  types,  material  codes  and 

numerical  integration  orders; 

-  SHELl  -  identify  selected  elements  as  thin  shells;  and 

-  CONTact  -  create  surface  contact  elements  from  certain 

surfaces  of  existing  three-dimensional  elements. 


The  PROPERTY  directive  is  the  usual  means  of  defining  the 
pertinent  property  information  for  existing  elements  of  the  model. 
Elements  may  be  assigned  properties  in  groups  (all  2-D  or  all  3-D 
elements,  random  lists  of  elements,  or  element  ranges),  with  a 
single  set  of  property  specifications  applying  to  all  elements 
in  the  group.  Elements  may  be  defined  as  being: 


Type  0 
Type  1 
Type  2 
Type  3 
Type  4 
Type  5 
Type  6 


Default  (solid,  plane  stress,  beam,  or  bar). 
Thin  shells. 

Plane  stress  (membrane)  elements, 

Plane  strain  elements, 

Axisymmetric  solids, 

Shear  panels,  or 
Surface  contact  elements. 


If  no  type  is  declared  for  an  element,  or  if  the  number  of  nodes 
is  inappropriate  for  the  type  assigned,  the  default  element  type 
for  the  existing  nodal  pattern  will  be  assigned  in  the  REFMT 
module  (see  Subsection  5.1).  Integration  orders  are  treated 
similarly,  and  generally  the  default  values  assigned  in  REFMT 
are  acceptable. 
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Material  property  codes  may  be  assigned  within  the  PROPerty 
utility,  using  material  codes  appearing  in  the  material  propertie 
library  described  in  Appendix  F.  Materials  data  can  also  be 
assigned  within  the  REFMT  processor,  and  existing  material 
coefficients  may  be  edited  as  well  (see  Subsection  5.1). 

The  SHELL  command  performs  a  similar  function  to  PROPERTY. 
Selected  elements  may  be  redefined  as  being  thin  shells  or  three 
dimensional  solids,  without  specifying  integration  orders  and 
material  property  codes. 

The  function  of  CONTact  is  complementary  to  that  of  PROPERTY: 
while  the  PROPERTY  command  allows  existing  elements  to  be  defined 
as  special  elements  to  be  used  in  surface  contact  analysis  (see 
Reference  [4]),  CONTACT  provides  a  means  of  creating  new  elements 
which  are  immediately  declared  to  be  surface  contact  elements. 

The  new  contact  elements  are  created  from  specified  faces  of 
existing  three-dimensional  solids. 

One  other  function  of  the  CONTACT  command  is  to  create 
elements  which  lie  on  the  faces  of  three-dimensional  solids, 
whose  properties  are  subsequently  redefined  using  the  PROPERTY 
command.  This  device  can  be  useful  when  thin  coatings  (e.g., 
protective  or  damping  layers)  are  to  be  applied  over  the  surface 
of  shell  or  solid  elements,  with  the  coating  layer(s)  to  be 
treated  as  membrane  (plane  stress)  finite  elements. 


4.3  boundary  conditions  add  constraints 

Modal  constraints  nay  be  applied  to  a  finite  elenent  nodel  in 
PRPP  us  inn  the  BOUNDS  and  LINEAR  commands.  Althouqh  boundary  con¬ 
ditions  nav  be  defined  at  any  tine  durinn  the  construction  of  a 
nodel ,  it  is  advisable  to  complete  most  of  the  neonetric  model  inn 
operations,  includinq  node  reorderinn  ,  before  the  constraints  are 
specified  (see  Subsection  4.8  for  typical  sequences  of  operations). 

L’ith  the  BOUNDS  operation,  honoqeneous  nodal  constraints  (i.e., 
displacement  =  0)  nav  be  specified  for  nrouns  of  nodes  selected  bv 
one  of  the  followinn  methods: 

-  select  all  nodes; 

-  select  all  nodes  on  a  plane  defined  by  Ax  +  By  +  Cz  =  P;  or 

-  select  all  nodes  in  a  niven  ranqe  (first,  last,  increment). 

Each  constraint  definition  benins  with  the  selection  of  one  of  the 
above  options  (or  selection  of  "no  more  boundary  conditions"). 

PREP  then  permits  all  of  the  selected  points  to  be  fixed  in  any 
combination  of  the  X,  Y  and  2  directions.  This  cycle  (select  nodes, 
define  constraint  direction)  continues  until  the  selection  option 
"no  nore  boundary  conditions"  is  specified.  BOUNDS  nay  be  used  as 
many  tines  as  required  to  define  all  simple  boundary  conditions  for 
the  nodel. 

The  LINEAR  function  allows  the  specification  of  simple  forms 
of  linear  constraints  at  individual  nodes.  This  form  of  constraint 
is  required  when  the  displacement  at  a  node  is  constrained  in  a 
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direction  which  does  not  coincide  with  the  global  x,  Y,  Z  axis  dir¬ 
ections.  Linear  constraints  are  defined  for  one  node  at  a  tine? 
first  a  node  nunber  is  selected  ,  and  then  the  direction  of  constraint 
is  defined.  The  direction  of  constraint  is  always  defined  by  nivinq 
the  plane  in  which  notion  is  permitted  at  the  node  (i.e.,  the  con¬ 
straint  is  applied  in  the  direction  normal  to  this  plane) .  "he 
plane  in  which  notion  occurs  can  be  defined  either  by  giving  the 
equation  of  the  plane  (in  the  form  Ax  +  Dy  +  Cs  ■  P ) ,  or  bv  listing 
the  coordinates  of  three  points  lyinq  in  the  plane.  As  with  nounns, 
the  LINEAR  operation  nay  be  performed  any  nunber  of  tines  to  define 
the  required  constraints. 


4.4  APPLIED  LOADING 


Applied  loadings  of  two  types  may  be  defined  in  PREP: 


-  concentrated  nodal  forces;  and 

-  distributed  element  loads,  including  body  forces,  surface 

pressures  and  edge  loads. 


Applied  loads  input  is  initiated  using  the  LOADS  command,  which 
may  be  given  as  many  times  as  necessary  in  the  construction  of 
a  model.  Five  basic  options  are  available  in  the  LOADS  function. 


They  are: 


-  (N)odal  loads  input; 

-  (E)lement  loads  input; 

-  (L)ist  currently  defined  loads  data; 

-  (H)elp  (print  a  list  of  available  options);  and 

-  (S)top  loading  data  input. 

These  options  may  be  executed  in  any  order  within  the  LOADS 
function . 

When  the  nodal  loads  option  (N)  is  selected,  PREP  will  accept 
data  defining  concentrated  forces  acting  at  existing  node  points 
of  the  model.  The  specification  of  nodal  loads  includes  the  set 
of  nodes  at  which  the  load  is  to  act,  a  loading  case  number,  and 
the  three  components  of  force  (X,  Y,  Z)  which  act  at  the  nodes 
selected.  Several  options  exist  for  specifying  the  node  numbers 
at  which  the  force  acts,  including: 

-  a  single  node; 

-  a  range  of  nodes  (first,  last,  increment);  and 

-  all  nodes  on  a  plane  defined  by  Ax  +  By  +  Cz  *  D. 

The  vector  force  (Fx,  Fy,  Fz)  input  is  assumed  to  act  at  each  node 
selected  by  one  of  the  above  methods.  It  should  also  be  noted 
that  the  'load  case'  number,  a  term  usually  associated  with 
linear  finite  element  analysis,  can  be  used  to  distinguish  sets 
of  loads  which  act  non-proportional ly  in  a  nonlinear  solution 
(see  REFMT,  Subsection  5.1). 

Element  loads  may  be  specified  for  a  single  element  at  a  time, 
a  range  of  elements  (defined  by  first  and  last  element  numbers  and 
an  element  number  increment),  or  all  elements  in  the  model.  For 
each  series  of  elements  selected,  a  single  load  specification  can 
be  defined  by  a  loading  case  number,  a  loading  type,  and  a  load 
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magnitude.  Loading  types  are  defined  as  follows: 

-1,  -2,  -3  :  denotes  a  body  force  (force  per  unit  volume), 

acting  in  the  global  X,  Y,  Z  directions, 
respectively; 

-  1,2,3, 4,5,6  :  for  3-D  elements,  signifies  a  distributed 

surface  pressure  acting  on  the  surface 
(-R, -S, -T, +R,+S,+T) ,  respectively; 

1,2, 3, 4  :  for  2-D  elements,  describes  an  edge  force 

(force  per  unit  length)  on  edges  1,  2,  3  or 
4,  respectively;  and 

1,2, 3, 4  :  for  axisymmetric  elements,  defines  a  surface 

pressure  on  edge  1,  2,  3  or  4  of  the  element. 

Surface  numbers  (for  pressure  loading  on  3-D  elements)  and  edge 
numbers  (for  line  loads  or  axisymmetric  pressures)  are  shown  in 
Figures  4.6  and  4.7.  The  correct  surface  or  edge  number  can  be 
identified  by  recalling  the  numbering  sequence  (  -R,  -S,  -T,  +R, 
+S,  +T  for  3-D,  and  -R,  -S,  +P.,  +S  for  2-D  )  ,  and  using  the  PLOT 
option  ORIENT  (see  Subsection  4.6)  to  display  the  local  coordinate 
directions  for  the  elements  in  question. 

4.5  NODE  POINT  REORDERING 

i 

Finite  element  models  generated  using  interactive  modelling 
techniques  are  typically  completed  with  most  of  the  nodal  and 
element  data  being  numbered  in  an  almost  random  fashion.  For 
efficiency  in  the  actual  finite  element  solution,  reordering  of 
either  node  points  or  elements  is  usually  advisable,  since  the 
number  of  operations  performed  in  the  solution,  as  well  as  the 
storage  requirements,  can  be  affected  to  a  significant  degree. 
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SURFACE 

LOCATION 

CORNER  NODES 

1 

R  *  -1 

2  -  3  -  7  -  6 

2 

S  *  -1 

3  -  4  -  8  -  7 

3 

T  3  -1 

5  -  6  -  7  -  8 

4 

R  *  +1 

1  -  4  -  8  -  5 

5 

S  3  +1 

1  -  5  -  6  -  2 

6 

T  3  +1 

1  -  2  -  3  -  4 

Figure  4.6.  Surface  Numbers  for  3-D  Pressure  Loads 


An  example  of  the  effect  of  reordering  a  small  finite  element 
mesh  upon  the  storage  requirement  for  solution  is  depicted  in 
Figure  4.8. 

Node  point  reordering  is  performed  in  PREP  using  the  RENUMBER 
command ,  which  requires  no  input  from  the  user  other  than  an  input 
file  label.  The  RENUMBER  option  performs  a  renumbering  of  the 
nodes  which  is  not  optimum,  but  which  is  usually  sufficient  to 
reduce  storage  and  solution  time  requirements  to  a  very  reasonable 
level.  In  particular,  the  renumbering  scheme  is  generally  quite 
effective  when  the  equation  solver  to  be  used  employs  a  variable 
bandwidth  (profile,  or  envelope)  storage  mode  [4], 

RENUMBER  displays  the  initial  and  final  'nodal  bandwidths'  as 
a  measure  of  how  effective  the  node  reordering  has  been.  The 
nodal  bandwidth  is  the  maximum  bandwidth  of  a  system  of  equations, 
based  on  the  mesh  in  question,  having  a  single  degree  of  freedom 
per  node;  the  actual  maximum  bandwidth  for  the  model  is  the  nodal 
bandwidth  multiplied  by  the  number  of  freedoms  per  node  (usually 
2,  3  or  6)  . 

The  RENUMBER  utility  is  quite  fast  in  execution,  and  can  be 
executed  more  than  once  to  attempt  further  reduction  in  the  system 
bandwidth.  Two  passes  through  RENUMBER  are  recommended,  to  verify 
that  the  model  bandwidth  has  been  reduced  as  much  as  possible. 


4.6  PLOTTING  UTILITIES 


Geometry  plotting  may  be  performed  in  PREP  at  any  time  by 
entering  the  PLOT  command.  PREP  requests  the  name  of  an  active 
data  file  to  be  plotted ,  and  enters  the  plot  command  mode 
(identified  by  the  prompt  symbol  . 
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In  the  PLOT  mode,  plotting  commands  may  be  entered  to  control 
the  appearance  of  plots  to  be  made.  For  example,  the  EYE  command 
sets  up  a  viewer  position  for  use  in  plotting,  but  the  plot  is  not 
actually  drawn  until  a  DRAW  command  is  issued.  Most  of  the  PLOT 
mode  commands  control  one  or  more  plotting  parameters,  which  may 
be  set  and  reset  as  often  as  desired.  These  parameter- related 
commands  include: 


-  AXES  : 

-  CLIP  : 

-  CUBE  : 

-  DEFAULT  : 

-  ELEMENT  : 

-  EYE  : 

-  LABEL  : 

-  ORIENT  : 

-  PROJECT  : 

-  REFLECT  : 

-  ROTATE  : 

-  SCALE  : 

-  SHRINK  : 

-  TRANSLATE: 

-  VERTICAL  : 

-  ZOOM  : 


control  plotting  of  X,  Y,  Z  axis  lines 

define  a  "clipping  plane"  which  limits  drawn  lines 

to  a  specified  minimum  distance  from  the  viewer 

define  minimum  and  maximum  plotting  coordinates 

restore  all  plotting  parameters  to  default  values 

select  a  range  of  elements  to  be  plotted 

define  the  viewer  eye  position 

select  node  and/or  element  numbering 

control  plotting  of  local  element  orientations 

needed  in  REFINE  and  LOADS  operations 

select  orthogonal  or  perspective  projection 

select  reflection  of  geometry  about  a  certain 

pi  ane 

define  rotations  (for  plotting  purposes  onxy) 

scale  or  unscale  plot 

select  exploded  view  to  be  plotted 

define  translations  (for  plotting  only) 

select  a  vertical  axis  direction 

select  an  area  of  the  screen  for  close-up  views 


Sufficient  default  values  are  set  upon  entry  to  PLOT  that  the 
model  may  be  drawn  without  setting  any  of  the  above  parameters, 
and  often  this  is  a  useful  procedure  for  "getting  oriented"  with 
a  new  model.  Once  any  of  the  above  parameters  has  been  set, 
however,  the  new  plotting  parameters  will  remain  active  even 
when  PLOT  is  terminated  and  reentered;  this  measure  eliminates 
the  need  for  entering  the  same  sequence  of  commands  every  time 
PLOT  is  executed. 
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The  remaining  commands  available  in  the  PLOT  module  are 
control  functions,  namely: 


-  BOUNDS 

Plot  symbols  indicated  boundary  conditions  at  node) 

-  DRAW 

Plot  the  model  using  all  current  parameter  values 

-  HELP 

Print  a  list  of  available  commands  and  options 

-  MAIN 

Return  control  to  PREP  following  plotting 

-  SITE 

Redefine  center  of  plotting  region 

-  SUMMARY 

Print  a  summary  of  all  plotting  parameter  values 

-  TIME 

Print  the  CPU  time  used  since  sign-on 

All  of  the  control  functions  request  some  form  of  immediate  action, 
with  no  further  input  required  from  the  user. 

Most  of  the  plotting  options  described  above  are  demonstrated 
in  the  modelling  examples  in  Section  7. 

4.7  DATA  LISTING  UTILITIES 

The  functions  available  for  listing  of  data  within  PREP  are  of 
two  types.  The  LIST  command  produces  a  table  of  active  model 
files,  with  their  corresponding  file  numbers.  The  PRINT  utility 
allows  data  from  a  specific  model  file  to  be  printed  at  the 
terminal  or  on  a  formatted  file  suitable  for  line  printer  output. 

A  typical  LIST  output  might  appear  as  follows: 

LISTING  OF  DATA  FILES  AVAILABLE 

#10 - PANEL 

#11 - PANELREF 

#13 - PANEL  +  BEAM 

Here  the  names  *  PANEL* ,  'PANELREF',  and  'PANEL  +  BEAM'  are  all 
names  of  active  data  files  which  may  be  specified  as  input  for 


any  PREP  operation.  The  corresponding  file  numbers  (10,  11,  13) 
are  important  only  on  CDC  computers,  where  they  correspond  to 
the  local  file  names  TAPE10,  TAPEll  and  TAPE 13 ,  respectively. 

Upon  leaving  PREP,  the  model  'PANEL  +  BEAM'  would  be  saved  by 
cataloging  the  local  file  TAPE13  as  a  permanent  disk  file,  for 
example.  In  the  VAX  11/780  version  of  PREP,  data  files  may  be 
saved  by  name  prior  to  leaving  the  program. 

With  the  PRINT  command,  the  contents  of  a  single  model  may  be 
listed  at  the  terminal  and/or  written  to  a  formatted  file  for 
off-line  printing.  Printing  of  header  datar  nodal  coordinates, 
element  data,  boundary  conditions,  and  nodal  or  element  loads 
may  each  be  selected  as  needed.  For  the  nodal  and  element  data, 
which  can  become  quite  lengthy,  ranges  of  nodes  or  elements  may 
be  specified  for  printing. 

4.8  TYPICAL  SEQUENCES  OF  OPERATIONS 

A  complete  listing  of  the  commands  and  options  available  in 
PREP  appears  in  Table  4.8.  These  commands  can  be  grouped 
logically  into  four  categories, 

-  utility  commands, 

-  geometric  modelling  functions, 

-  boundary  conditions  and  loads,  and 

-  property  definition. 

The  utility  commands  are  those  which  perform  data  management 
or  informative  functions;  generally,  they  may  be  executed  at  any 
time,  in  any  sequence.  The  commands  which  fall  into  this 
category  include: 


TABLE  4.8 


Alphabetized  Listing  of  PREP  Commands 


COMMAND 


DESCRIPTION 


BOUNDS 

CONTACT 

COPY 

CREATE 

DELETE 

EDIT 

FILL 

HELP 

LINEAR 

LIST 

LOADS 

MASK 

MERGE 

NAME 

PLOT 

PRINT 

PROPERTY 

REFINE 

REFLECT 

RENUMBER 

ROTATE 

SHELL 

SIFT 

STOP 

TRANSLATE 

TIDY 

TIME 

TOLERANCE 


Specify  nodal  constraints 

Create  contact  elements  from  surfaces  of  3-D 
solids 

Make  a  copy  of  a  specified  model  file 

Begin  a  new  model  file  by  keyboard  entry  of  data 

Delete  a  model  file 

Make  minor  editing  changes  to  an  existing  model 

Generate  midside  and/or  interior  nodes 

List  available  commands  and  options 

Generate  linear  constraint  conditions  at  nodes 

List  all  active  model  files  and  file  numbers 

Define  nodal  forces  and/or  distributed  element 

loads 

Eliminate  midside  and/or  interior  nodes 
Combine  two  model  files 

Assign  a  descriptive  name  to  a  file  number 
Enter  plotting  mode 

Print  all  or  part(s)  of  a  model  data  file 
Assign  element  property  data 
Subdivide  elements  to  create  a  finer  mesh 
Form  the  mirror  image  of  an  entire  model 
Reorder  node  points  to  reduce  solution  time 
Rotate  model  geometry  in  space 
Select  certain  elements  to  be  thin  shells 
Eliminate  unused  node  points 
Exit  from  PREP 

Translate  model  geometry  in  space 

Eliminate  duplicate  node  points 

Print  CPU  time  since  start  of  session 

Reset  tolerance  for  detecting  coincident  nodes 


COPY 

DELETE 

HELP 

LIST 


NAME 

PLOT 

PRINT 

TIME 


Property  definiton  functions,  which  are  also  relatively  benign 
with  respect  to  their  execution  sequence,  include 

CONTACT 

PROPERTY 

SHELL 


However,  it  may  be  advantageous  in  some  applications  to  segregate 
elements  of  one  type  or  property  set  from  those  of  another  early 
in  the  model  generation,  when  fewer  elements  are  involved. 
Definition  of  many  different  sets  of  properties  in  a  complicated 
model  presents  a  potential  source  of  errors  and  oversights. 

With  the  remaining  two  groups  of  commands,  esp*  ir  ience  has 
shown  that  some  sequences  of  commands  result  in  very  fast  and 
effective  modelling,  while  others  can  lead  to  long  and  tedious 
sessions  at  the  terminal.  It  is  generally  recommended  that  the 
definition  of  loads  and  boundary  condition  data  be  performed  as 
the  LAST  step  in  model  generation,  after  such  operations  as 
REFINE,  RENUMBER,  etc.  The  geometric  modelling  operations  also 
have  a  preferred  order,  which  can  best  be  described  as  three 
'stages'  of  the  geometric  model  generation.  Commands  associated 
with  these  modelling  stages  are  summarized  below: 

Stage  1  :  FILL  MASK 

Stage  2  :  MERGE  REFINE  REFLECT 

ROTATE  TOLERANCE  TRANSLATE 
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Stage  3  :  RENUMBER 


In  Stage  1,  irregular  elements  entered  in  CREATE,  and/or 
higher  order  elements  which  are  to  be  reduced  to  lower-degree 
nodal  patterns  are  adjusted  to  contain  the  final  number  of  nodes 
per  element  which  is  to  be  used  throughout  the  mesh.  Once  these 
steps  have  been  accompl ished ,  other  operations  such  as  REFINE 
and  MERGE  should  proceed  smoothly. 

In  Stage  2,  the  actual  mesh  division  to  be  used  for  analysis 
is  created,  by  REFIning  a  coarse  mesh,  performing  TRANslations 
and  ROTAtions  to  reproduce  periodic  sections  of  the  grid,  and 
MERGing  disjoint  sections  of  the  model  together.  Since  most  of 
the  operations  in  this  phase  of  the  modelling  will  affect  the 
number  and/or  ordering  of  nodal  coordinates  and  elements, 
RENUmbering,  BOUNdary  conditions  and  LOADs  are  best  left  until 
this  stage  is  complete. 

In  Stage  3,  the  geometry  of  the  mesh  (number  of  nodes  and 
elements  and  their  relationship  to  one  another)  is  complete,  and 
RENUmber  can  be  used  to  minimize  the  effects  of  the  almost 
random  numbering  of  the  finite  element  model. 

As  the  final  step  in  generating  the  model,  boundary  condi¬ 
tions  and  loading  data  are  entered  in  PREP.  At  this  point,  the 
geometry  of  the  model,  as  well  as  the  numbering  of  nodes  and 
elements,  is  fixed;  PLOT  should  be  used  to  obtain  the  node  and 
element  numbers  needed  for  specifying  boundary  conditions  and 
loading  . 
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SECTION  5 


MODEL  DATA  TRANSLATION 


The  preprocessor  contains  two  facilities  for  producing 
formatted  output  from  PREP-generated  finite  element  models. 

These  are  the  REFMT  and  NEUTRAL  utilities,  which  are  described 
in  Subsections  5.1  and  5.2,  respectively. 

REFMT  provides  output  in  a  format  compatible  with  the  MAGMA 
nonlinear  finite  element  program  [4].  Additional  data  required 
by  MAGNA  which  is  not  prepared  in  PREP  can  also  be  generated  in 
REFMT  interactively.  In  most  instances,  the  REFMT-generated 
data  deck  is  ready  for  submission  directly  to  MAGNA  for  analysis. 
When  special  options  or  minor  data  changes  are  necessary,  the 
formatted  data  file  can  be  modified  'by  hand'  using  the  resident 
computer's  editing  facilities. 

The  NEUTRAL  utility  translates  the  PREP- format  data  file  into 
a  formatted  form  which  is  listable  on  the  line  printer  and  may  be 
transferred  between  computer  systems  with  a  minimum  of  difficulty. 
NEUTRAL  also  performs  the  inverse  translation  function  (back  to 
PREP  internal  format)  needed  to  re-introduce  the  finite  element 
data  into  the  preprocessor  on  another  computer. 


5.1  REFMT  -  MAGNA  Data  Translation 


The  REFMT  translator  accepts  as  input  the  PREP  internal  format 
data  file,  and  generates  a  formatted  data  deck  suitable  for  use 
with  the  finite  element  analysis  program  MAGNA  [4].  All  of  the 
finite  element  types,  and  most  of  the  special  options,  available 
in  MAGNA  are  supported  by  REFMT. 

In  addition  to  reformatting  the  preprocessor-generated  data, 
REFMT  requests  information  from  the  analyst  to  complete  the 
solution  data  needed  in  MAGNA.  This  additional  information 
includes : 

-  the  type  of  solution  (linear  or  nonlinear,  static  or  dynamic), 

-  information  concerning  restart  files  for  nonlinear  analysis, 

-  postprocessing  file  options,  and 

-  nonlinear  solution  strategy  (iteration  method,  etc.). 

The  data  required  for  multiple  loading  conditions,  contact 
analysis  and  other  options  is  generated  automatically  from  the 
preprocessor  data. 

Additional  editing  of  the  model  data  is  possible  in  REFMT  for 
material  properties  and  loading  data.  Material  properties  may  be 
selected  from  a  library  of  stored  properties  (see  Appendix  F) ,  or 
entered  directly  at  the  keyboard.  When  the  element  properties  are 
complete  for  all  elements  in  the  model,  materials  data  is  listed 
and  can  be  modified  selectively  as  the  need  arises. 


Applied  loading  data  may  take  on  different  meanings  in  linear 
and  nonlinear  analysis,  and  for  transient  dynamic  and  nonlinear 
solutions,  the  time  history  of  the  loading  must  also  be  specified. 
REFMT  allows  a  moderate  amount  of  editing  of  the  loading  data  to 
suit  the  type  of  analysis  to  be  performed.  The  options  include: 

-  options  to  combine  multiple  load  cases,  or  analyze  separately; 

-  generation  of  time-history  curves  for  use  in  nonlinear  and/or 
dynamic  analysis  (several  built-in  curves  may  be  selected  and 
scaled  as  needed);  and 

-  specification  of  non-proportional  loading  sets  in  nonlinear 
and  transient  analysis,  with  each  set  scaled  independently. 

Most  of  the  options  available  within  REFMT  are  demonstrated  in 
the  sample  preprocessing  sessions  included  in  Section  7.  It  should 
also  be  noted  that  the  MAGNA  data  deck  (whether  output  from  REFMT_ 
or  hand-generated)  may  be  translated  back  into  the  preprocessor 
file  formats  using  the  TRNSFR  utility  described  in  Subsection  3.1. 


5.2  NEUTRAL  -  Neutral  File  Translation  Program 

The  NEUTRAL  utility  performs  two  data  translation  functions: 

-  preprocessor  internal  format  to  neutral  data  file,  and 

-  neutral  data  file  to  preprocessor  internal  format. 

These  two  functions  are  inverses  of  one  another,  with  all 
preprocessor  data  types  supported. 


The  two  data  files  (internal  format  and  neutral  format)  are 
described  in  Subsections  6.1  and  6.2,  respectively.  The  neutral 
file  format  exists  solely  for  the  purpose  of  storing  the  prepro¬ 
cessor  data  in  a  readable  form  which  is  easily  transferred 
between  computer  systems.  Thus,  finite  element  data  prepared 
using  PREP  on  one  system  may  be  translated  to  the  neutral  format 
using  NEUTRAL,  transmitted  to  a  completely  different  machine 
type,  returned  to  the  internal  format  using  NEUTRAL,  and  used 
within  PREP  on  the  second  machine.  The  translation  to  neutral 
format  is  also  useful  for  archiving  of  data,  either  on  punched 
cards  or  in  text  library  format  (using  the  LIBRARY  utility  on 
VAX  11/780,  or  UPDATE  on  CDC  machines) . 

Execution  of  NEUTRAL  requires  only  the  declaration  of  the 
files  involved,  and  the  selection  of  the  type  of  translation  to 
be  performed.  On  CDC  machines,  the  internal- format  data  file  is 
local  file  UNFMT,  and  the  formatted  (neutral)  file  is  the  local 
file  FMTDAT.  With  the  VAX  version  of  NEUTRAL,  the  necessary 
file  names  are  requested  as  input  at  the  beginning  of  execution. 


SECTION  6 


MODEL  DATA  FILES 


Optimal  data  handling  of  information  associated  with  model 
geometry  and  finite  element  analysis  is  highly  dependent  on  how 
that  data  is  accessed  and  stored.  The  preprocessing  system 
utilizes  several  types  and  formats  for  data  files  to  take 
advantage  of  specific  features  offered  by  each.  Data  that  is 
transferred  between  the  various  preprocessing  modules  in  the 
stages  of  model  development,  modification,  anc7  assignment  of 
analysis  parameters  is  stored  in  an  internal,  unformatted  data 
file.  This  binary  data  file  is  only  readable  by  programs  that 
are  properly  equipped  to  read  the  particular  file  in  question. 

Two  different  data  files  are  employed  in  the  internal,  unformatted 
data  type;  one  for  the  output  data  generated  by  CORGEN,  SPATCH, 
or  AGRID,  and  a  second  data  type  generated  by  PREP,  IJKGEN, 

CREATE,  or  NEUTRAL.  The  second  type  of  data  file  used  in  the 
preprocessor  system  is  the  external,  formatted  data  file.  This 
data  file  type  is  used  for  transferring  model  data  between 
preprocessing  programs,  and  for  archival  of  data  in  a  readable 
form.  The  external,  formatted  data  file  is  a  character 
text  file  and  may  therefore  be  read  without  any  need 
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for  special  programs.  General  descriptions  of  the  form  and 
purpose  of  each  of  these  data  file  types  are  given  in  the 
remainder  of  this  Section.  Detailed  descriptions  of  the  file 
formats  for  each  type  are  presented  in  Appendix  A. 

6.1  Preprocessor  Internal  Data  Structure 

The  preprocessing  system  defines  two  functionally  similar 
data  structures  for  storing  model  geometry  data.  Execution  of 
any  of  the  data  entry  and  interface  programs  discussed  in 
Sections  2  and  3  results  in  the  creation  of  either  a  two- 
dimensional  surface  model  geometry  file  (see  CORGEN,  AGRID,  and 
SPATCH)  or  a  three-dimensional  solid  model  geometry  file  (see 
CREATE,  IJKGEN,  NEUTRAL,  TRNSFR,  and  IMPRINT) .  The  surface 
model  geometry  file  is  an  intermediate  model  file  designed  to 
be  input  to  the  EXPAND  module  to  convert  the  nine  node  surface 
element  model  to  a  twenty-seven  node  solid  element  model.  The 
shell  surfa.ce  file  and  the  standard  (PREP- format)  3-D  geometry 
file  are  the  subjects  of  Paragraphs  6.1.1  and  6.1.2,  respectively. 

6.1.1  Shell  Surface  Geometry  File 

Certain  modelling  requirements  for  windshield  canopies  and 
related  structures  fall  into  a  class  where  the  structures  may  be 
represented  by  2-D  quadratic  surface  elements  located  in  3-D  space. 
If  a  thickness  value  is  associated  with  the  surface  then  it 
becomes  a  straightforward  task  to  mathematically  expand  the 


6.2 


«r* 


surface  to  a  3-D  solid  representation.  The  CORGEN,  AGRID  and 
SPATCH  programs  generate  simplified  surface  models  which  can  be 
readily  represented  by  the  data  file  described  in  Appendix  A, 

Table  A.l.  This  data  file  contains  coordinate  data  with  an 
associated  thickness  value  for  each  node,  and  element  connectivity 
for  2-D  bi-quadratic  surface  elements.  The  surface  elements  are 
defined  by  nine  nodes:  four  corner  nodes,  four  midside  nodes  and 
a  center  node.  The  expand  module  will  convert  the  nine  node 
surface  element  model  to  three-dimensional  solid  finite  elements 
stored  in  the  three  dimensional  model  file  format  described 
below.  This  2-D  file  is  characterized  by  the  implicit  node  numbering 

of  coordinates  and  the  implicit  numbering  of  the  elements. 

6.1.2  Three-dimensional  Model  File 

The  requirements  of  maintaining  data  neccessary  for  model 
geometry  as  well  as  for  numerical  analysis,  determine  the  structure 
of  the  three-dimensional  model  file.  Table  A. 2  shows  the  details 
of  this  file  structure.  This  model  file  is  the  workhorse  of  the 
preprocessing  system  in  that  all  models  that  utilize  the  PREP  and 
REFMT  modules  must  be  converted  to  this  data  structure.  Particular 
features  of  this  model  file  include  the  arrangement  of  element 
loads,  nodal  forces  and  boundary  conditions  in  such  a  manner  that 
this  data  need  not  be  provided  nor  must  file  space  be  allocated 
to  provide  for  it  later.  While  this  file  type  is  similar  to 
the  surface  model  data  structure  the  user  should  be  aware  of  the 
explicit  node  and  element  numbering  that  is  included  with  this 
data  structure  along  with  the  loss  of  the  thickness  data  required 


for  EXPAND  and  the  addition  of  element  codes  that  regulate  and 
affect  the  numerical  solution  options.  This  file  is  used  as 
input  to  several  of  the  preprocessing  modules  as  illustrated  in 
Figure  1.1.  Once  all  preprocessing  functions  have  been  executed 
for  model  definition  this  file  is  input  to  the  REFMT  module  to 
convert  it  to  a  reformatted  character  file  with  the  correct 
options  for  input  to  the  appropriate  analysis  program. 

6.2  Neutral  Data  File  Structure 

The  neutral  data  file  is  an  external,  formatted  data  file 
intended  for  use  in  archiving  finite  element  data  and  transmitting 
model  data  from  one  computer  system  to  another.  The  neutral  data 
file  is  read  and  written  only,  by  the  NEUTRAL  utility  (Subsection 
5.2),  which  performs  the  translation  between  neutral  files  and  the 
internal  format  of  the  preprocessor  (Subse.ction  6.1). 

The  neutral  file  contains  exactly  the  same  information  as  the 
standard  preprocessor  data  file  (see  Appendix  A).  However,  the 
fact  that  the  neutral  file  is  a  formatted  (text)  file  has  several 
useful  consequences: 

-  data  may  be  archived  in  neutral  file  form  using  the  LIBRARY 
functions  provided  by  VAX/VMS,  or  UPDATE  on  CDC  machines; 

-  alternatively,  model  data  may  be  punched  directly  on  cards 
for  later  use;  and 

-  transmission  of  data  via  RJE  terminals  can  be  accomplished 
with  a  minimum  of  difficulty. 


The  data  sequence  and  formats  for  the  neutral  data  file  are  listed 
in  Appendix  A  of  this  report. 


6.3  MAGNA  Input  Data  File 

The  MAGNA  [4]  finite  element  analysis  program  requires  the 
user  to  provide  data  for  an  analysis  in  a  particular  formatted 
form.  While  the  data  stored  in  the  internal  unformatted  form 
is  adequate  for  preprocessing  functions,  this  data  must  be 
converted  to  the  appropriate  structure  for  input  to  MAGNA.  The 
MAGNA  input  data  file  is  created  by  the  REFMT  program  in  response 
to  user  selected  solution  options.  This  file  contains  the 
complete  information  concerning  solution  options,  nodal  coor¬ 
dinates,  element  connectivity,  boundary  conditions,  loading 
requirements  and  all  associated  parameters  for  element  types, 
material  properties,  time  increments,  etc.  All  data  on  this 
file  is  formatted  according  to  the  structure  outlined  in  detail 
in  Reference  [4].  Due  to  the  formatted  nature  of  the  MAGNA 
data  file,  the  user  may  readily  edit  it  directly  to  alter  any 
values  that  may  have  been  improperly  specified  in  the  REFMT 
program.  In  addition,  the  user  is  encouraged  to  make  changes 
to  the  file  directly,  in  the  event  minimal  changes  are  necessary 
for  altering  solution  options.  Reference  [4]  will  be  necessary 
for  the  user  wishing  to  alter  data  on  this  file. 


SECTION  7 


DEMONSTRATION  PROBLEMS 


Several  modelling  examples  are  presented  in  this  Section  which 
illustrate  the  capabilities  of  the  preprocessor  and  the  types  of 
procedures  involved  in  generating  a  finite  element  model.  Three 
of  the  data  entry  programs  (CREATE,  CORGEN,  and  IJKGEN) ,  which 
represent  the  most  common  points  of  entry  into  the  system,  are 
used  for  the  creation  of  the  original  models.  In  the  case  of 
CORGEN,  which  generates  a  surface  model  ..file,  the  use  of  EXPAND 
to  create  a  fully  three-dimensional  model  is  also  demonstrated. 

In  each  example,  further  modelling  operations  are  performed 
within  PREP,  and  the  final  data  is  translated  into  the  MAGNA 
input  format  using  REFMT. 

7.1  EXECUTION  OF  CREATE  /  PREP  /  REFMT  SEQUENCE 

The  following  example  illustrates  the  generation  of  a  two 
dimensional  finite  element  model  using  the  CREATE  data  entry 
program.  The  model  is  completed  using  PREP,  and  formatted  data 
for  analysis  with  MAGNA  is  generated  using  REFMT. 


The  model  generated  below  involves  a  thin,  flat  rectangular 
plate  with  a  circular  hole  at  the  center;  only  one  quadrant  of 
the  structure  is  represented  due  to  symmetry.  Of  particular 
interest  is  the  technique  of  using  only  selected  midside  nodes 
for  input  in  CREATE  to  capture  the  curvature  of  the  boundary 
along  the  hole.  This  method  allows  a  minimum  of  hand-generated 
input,  but  requires  the  use  of  FILL  (or  MASK)  within  PREP  to 
generate  a  properly  connected  mesh  of  elements. 


7.1.1  Summary  of  Modelling  Procedure 

The  basic  geometry  of  the  plate  is  defined  using  CREATE,  as 
detailed  in  the  sample  terminal  session..of  Section  7.1.2.  Ten 
node  points  and  three  elements  are  used  to  define  the  overall 
geometry.  Note  that  only  the  corner  nodes  are  used,  except  where 
curved  boundaries  are  to  be  defined.  It' is  important  to  observe 
that,  if  only  selected  midside  nodes  are  used,  the  ordering  of 
element  connectivity  may  contain  zeroes  in  positions  corresponding 
to  unused  nodal  positions.  This  situation  occurs  in  the  first 
two  elements  of  the  input  mesh. 

The  primary  operations  performed  in  PREP  are  the  FILL  option, 
in  which  additional  node  points  are  generated  to  fill  out  the 
midside  and  element  centroid  positions,  and  a  series  of  REFIne 
operations  which  subdivide  each  of  the  three  major  portions  of 
the  mesh  to  obtain  a  final  model  containing  80  nine-node  elements. 


The  REFMT  step  specifies  a  linear,  static  analysis.  Material 
properties,  which  were  not  defined  within  PREP,  are  also  selected 
from  the  materials  library  during  execution  of  REFMT.  Following 
the  execution  of  REFMT,  the  MAGNA  input  file  (local  file  FDATA) 
is  saved  as  the  permanent  file  PLATEP  on  an  indirect  file  library. 


7.1.2  CREATE  Execution 


COMMAND-  ATTACH .  P, PREPROCESSORPROC,  ID*iROC<HAN,MR*l  . 

AT  OY*  998  SN-AFFDL 
COMMAND-  8EGIN.CRCATE.P. 


*************  TSS5555555S555S55S55S5u555u 
BEGIN  CREATE  -  DATA  ENTRY  AND  EDITING  OF 
COARSE  GRID  FINITE  ELEMENT  MODEL  DATA 
**  **  *  tsnmmtxtmumnxitiuxnm 

DO  YOU  HAUE  A  PREVIOUSLY  GENERATED  DATA 
FILE  FROM  CREATE  TO  RE-EDIT?  (V.N) . <  N 


tXXtXSSSSXSSSX-  INITIAL  NOK  POINT  INPUT  XXXXXXX******* 

IN  THIS  PART  OF  THE  INPUT,  COORDINATE  DATA  IS  ENTERED  AT  THE 
KEYBOARD  FOR  EACH  NODE  OF  THE  MODEL.  THE  DATA  IS  CHECKED  FOR 
CONSISTENCY  AND  NAY  BE  EDITED  LATER.  THE  MAXIMUM  ALLOUABLE 
NODE  POINT  NUMBER  IS  SBB. 

ENTER  THE  DIMENSIONALITY  OF  THE  MESH  TO 
BE  DEFINED  (  2  OR  3  )  . 2 

txxxs- BEGIN  DIRECT  NODAL  POINT  INPUT  - 
ENTER  NODE  NO.  AMO  2  COORDINATES  AT  EACH  NODE 
(  ENTER  ALL  ZEROES  TO  TERMINATE  INPUT  > 


ENTER  NOK. 

X. 

V 

-  1 

B. 

S. 

ENTER  NOK. 

X. 

Y 

-  2 

1.87 

4.6 

ENTER  NOK, 

X, 

Y 

-  3 

3.SS 

3.SS 

ENTER  NOK, 

X. 

Y 

-  4 

4.6 

1.87 

ENTER  NOK, 

X, 

Y 

-  5 

S. 

B. 

ENTER  NOK. 

X. 

Y 

-  6 

IB. 

B. 

ENTER  NOK. 

X. 

Y 

-  7 

19. 

IB. 

ENTER  NOK, 

X. 

Y 

-  8 

B. 

IB. 

ENTER  NOK. 

X, 

Y 

-  9 

B. 

29. 

ENTER  NOK. 

X. 

Y 

-IB 

IB. 

29. 

ENTER  NOK. 

X. 

V 

“B.B.B 

EDIT  THE  NODAL  POINT  DATA  <  Y  ,  N  ) . N 

*Xt*XX*X*XXXX*X  INITIAL  ELEMENT  INPUT  XHXXXXXSXXXXXSX 

THIS  PORTION  OF  THE  INPUT  GIVES  THE  CONNECTION  DATA  DEFINING 
FINITE  ELEMENTS  IN  TERMS  OF  THE  NODES.  CONNECTIVITY  FOR  THE 
ELEMENTS  FOLLOWS  THE  CONVENTIONS  FOR  TWO  AND  THREE  DIMENSIONAL 
ELEMENTS  IN  MAGNA  - 


......  TWO-DIMENSIONAL  . . THREE-DIMENSIONAL . 

NODE  POINTS  LOCATIONS  NODE  POINTS  LOCATIONS 

1-4  CORNERS  1  -  B  CORNERS 

S  -  t  HIDSIDCS  9  -  2S  MIDSIDES 

9  CENTROID  21-26  NIDFACES 

27  CENTROID 

THE  MAXIMUM  ALLOWABLE  ELEMENT  NUMBER  IS  IBB. 

FOR  EACH  ELEMENT,  ENTER  (DELEMENT  NUMBER,  (2)NUMBER  LOCAL 
NODES  ( I .E. .THE  LENGTH  OF  THE  CONNECTIVITY  LIST)  AND  (3)THE 
LIST  OF  COWCCTED  NODES.  ENTER  ELEHENT-S  TO  TERMINATE  INPUT 

FOR  2-D  CLEMENTS.  THE  MAX.  LOCAL  NODE  NUMBER  SHOULD  BE  BETWEEN  2  AND  9 

ENTER  ELEMENT  NO. -I 
MAX.  LOCAL  MOOE  -4 
CONNECTIVITY  LIST-1. 1.3,7,B,2 

ENTER  ELEMENT  NO. -2 
MAX.  LOCAL  NOOE  -6 
CdOCCTIVITY  LIST-7,  3,6,6, 9,4 

ENTER  ELEMENT  HO. -3 
MAX.  LOCAL  NOK  -4 
CONNECTIVITY  LIST-10.9.S.7 

ENTER  ELEMENT  NO.-B 
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7.1.3  PREP  Execution 


co' hand  files 

—  uOCAL  FILES— 

*°  BINPUT  JOOTPUT  UHFfIT 

COMMAND-  REWIND,  (FHT 
COMMAND-  COPYBF.UXFMT.TAPEie 
EOI  ENCOUNTERED  AFTER  COPY  OF  FILE 
e.  RECORD  1 
COMMAND-  RETURN. UNF NT 
COMMAND-  BEGIN. PREP. P. 


nnttmtttnimiuimtmuimmn 
nmnntnt  begin  prep  mmimm 

BBSS  MAGNA  PREPROCESSING  UTILITIES  BBBBB 


PLOTTING  TERMINAL  (Y.N)  ?nY 
TEKTRONIX  OR  HP  <  T  ,  H  )  ?  ,jT 

TEKTRONIX  TERMINAL  TYPES  - 

S.  444S-1 

1.  4414  /  4412  /  4413  /  4452 

2.  4414  /  4415 

3.  4414  /  4415  <  ENM.  GR.  NOD.  ) 

4.  4114 

ENTER  TERMINAL  TYPE  (4-4  >tj3 
ENTER  CHARACTERS  PER  SECOND  |jl24 
NUMBER  OF  INPUT  DATA  FILES  ?  Jjl 
ENTER  FILE  B  jjlS 

ENTER  LABEL  ,, PLATE 

TYPE  'HELP'  FOR  LIST  OF  COMMANDS 

»l  PLOT 

ENTER  INPUT  FILE  41  LABEL  PLATE 

INITIALIZATION  OF  PLOTTING  ROUTINE  EJ5  SES  SCREEN. 
READY  (Y.N)  ?  ....I  Y 


(HELP 

OOMNANO  lESCRIPTION 

(tutu  ltll«B»»I 

AXES  AXES  DRAW  AND  LABEL 

BOUNDS  PLOT  BOUNDARY  CONDITIONS 

CLIP  CLIP  PLANE  POSITION 

CUBE  SET  MINIMA  AND  MAXIMA 

DEFAULT  SET  DEFAULT  UaLUES 

KMU  DRAU  MODEL 

ELEMENTS  PLOTS  ALL  OR  SELECTED  ELEMENTS 

EYE  EYE  POSITION 

HELP  LIST  ALL  COMMANDS 

LABELS  LABELS  ELEMENTS  ANOsOR  NODES 

MAIN  RETURNS  CONTROL  TO  MAIN  PROGRAM 

ORIENT  PLOTS  R.S.T  ORIENTATION  AXES 

PROJECT!®??  PROJECTION  TYPE 
REFLECT  REFLECT  A  PLANE 

ROTATE  ROTATE  MODEL  ABOUT  AXES 

SCALE  SCALE  PLOT 

SHRINK  SHRINK  ELEMENTS 

SUMMARY  LIST  ALL  P  ARAIK  TEA  UALUES 

TIME  PRINT  CPU  TIME  SINCE  START  OF  SESSION 

TRANSLATE  TRANSLATE  MflJEL  FROM  ORIGIN 

UERT1CAL  VERTICAL  AXIS 

ZOOM  ZOOM  ON  THE  MODEL 

SITE  SITE  POSITION 

SUES 

PLOT  AND  LABEL  THE  AXEST(V.N)  . >  Y 

BLASE 

LABEL  THE  ELEMENTS?  (Y.N)  . «  Y 

LABEL  THE  NOOES?(V,N>  . . 

LABEL  ALL  THE  SURF  ACES?  (Y.N  I  . I  Y 

SQRIENT 

PLOT  ORIENTATION  AXEST(Y.N)  .  V 


AT  END  OF  PLOT,  ENTER 

READY  (Y.N)  ?  ....<  V 


CHARACTER  TO  CONTINUE. 


•  i  FIU. 

ENTER  INPUT  FILE  *1  LABEL  j  j PLATE 

ENTER  OUTPUT  FILE  LABEL  jj PLATE 

CREATING  NEW  DATA  FILE 
FILL  OPTIONS  - 

1.  3-D  ELEMENTS  ONLY 
E.  S-D  ELEMENTS  ONLY 
ENTER  FILL  OPTION  <l.S)....;jE 
MAX  NODE  POSITION  B  IN  FINAL  ELEMENTS? 

<E,9>...j|9 

FILL  COMPLETE. 

PERFORM  SUBSEQUENT  TIDY  ?  <V#N>  ->Y 

TIDY  OPTION  ONE  OR  TWO  ?  (OPT1,  OPTS,  HELP)  ->OPTS 

TIDY  COMPLETE 
2  DUPLICATE  NODES  DELETED 

««  PLOT 

ENTER  INPUT  FILE  SI  LABEL  j  j PLATE 

SDRAM 

AT  END  OF  PLOT,  ENTER  CHARACTER  TO  CONTINUE. 
READY  (V.N)  ?  ....I  Y 


«'  <1 


::  .I-T 

‘LUTING  OF  DATA  FILES  AVAILABLE 

•1» - PLATE 

•it  —  PLAT2 

•«  DELE 

ENTER  LABEL  PLATE 

it  REFINE 

ENTER  INPUT  FILE  B1  LABEL  jjPLATS 

ENTER  OUTPUT  FILE  LABEL  j  jREFINl 

CREATING  NEW  DATA  FIU 

BEGIN  REFINE  ROUTINE... 

t  -  RANDOM  ELEMENT  INPUT 

2  -  RANGE  OF  ELEMENTS 

3  -  EXECUTE  REFINEMENT 

ENTER  OPTION - >1 

MAXIMUM  NUMBER  OF  ELEMENTS 

SPECIFIABLE - >  1SB 

NUMBER  STILL  AVAILABLE  - >  15* 

ENTER  ELEMENT  NUMBER  AND  <RETURtO  FOR  EACH 
ELEMENT  TO  BE  REFINED.  ELEMENTS  NEED  NOT  BE 
ENTERED  IN  ASCENDING  ORDER.  ELEMENT  NUMBER  • 
TERMINATES  INPUT. 

ENTER  ELEMENT  B - >1 

ENTER  ELEMENT  • - >2 

ENTER  ELEMENT  t - >• 

1  -  RANDOM  ELEMENT  INPUT 

2  -  RANGE  OF  ELEMENTS 

3  -  EXECUTE  REFINEMENT 

ENTER  OPTION - >3 

ENTER  CUT  DIRECTION  (R.S.T)  >R 


MAXIMUM  CUTS  PER  ELEMENT 
ENTER  CUTS  PER  ELEMENT  — 


->  S 
->S 


EQUALLY  SPACED  CUTS  ?  <Y,N>  >N 

ENTER  CUT  POSITIONING  IN 

ASCENDING  ORDER  (O-ISO)  — >1Q  22  35  SS  7S 

REFINE  COMPLETE. 

PERFORM  SUBSEQUENT  TIDY  ?  <Y,N>  ->Y 

TIDY  OPTION  ONE  OR  TUO  ?  (OPT1,  0PT»..HE$P>  ->0PT2 


TII'V  COMPLETE 
It  DUPLICATE  NODES  DELETED 

««  PLOT 

ENTER  INPUT  FILE  *1  LABEL  REF  INI 

(DRAW 

AT  END  OF  PLOT,  ENTER  CHARACTER  TO  CONTINUE. 
READY  <Y,N>  ?  ....«  V 


1  1.1 1 N 


It  REFINE 

ENTER  INPUT  FILE  *1  LABEL  j jREFINt 

ENTER  OUTPUT  FILE  LABEL  jjREFIN2 

CREATIN'!  NEU  DATA  FILE 

BEGIN  REFINE  ROUTINE... 

1  -  RANDOM  ELEMENT  INPUT 

2  -  RANGE  OF  ELEMENTS 

3  -  EXECUTE  REFINEMENT 

ENTER  OPTION - >2 

MAXIMUM  NUMBER  OP  ELEMENTS 

SPECIFIABLE - >  IS* 

NUMBER  STILL  AVAILABLE  - >  1SB 

ENTER  BEGINNING  ELEMENT  NUMBER.  ENDING  ELEMENT 
NUMBER.  AND  INCREMENT  FOR  RANGE. 

ENTER  NBEG.NEND.  I  NCR - >1.13.1 

1  -  RANDOM  ELEMENT  INPUT 

2  -  RANGE  OF  ELEMENTS 

3  -  EXECUTE  REFINEMENT 

ENTER  OPTION  - >3 

ENTER  CUT  DIRECTION  (R.S.T)  >S 

MAXIMUM  CUTS  PER  ELEMENT  - >  S 

ENTER  CUTS  PER  ELEMENT  - >4 

EQUALLY  SPACED  CUTS  ?  (Y.N)  >Y 

REFINE  COMPLETE. 

PERFORM  SUBSEQUENT  TIDY  ?  (Y.N)  ->Y 
TIDY  OPTION  ONE  OR  TUO  ?  COPTI,  OPT2.  HELP)  ->0PT8 
TIDY  COMPLETE 

119  DUPLICATE  NODES  DELETED 


II  PLOT 

ENTER  INPUT  FILE  SI  LABEL  if REFINE 

•ORIENT 

PLOT  ORIENTATION  AXEST(Y.N)  .  I 

SLABCL 

LABEL  THE  ELEMENTS?  (Y.N)  . . 

LABEL  THE  NODES? (Y.N)  . I  I. 


AT  END  OF  PLOT.  ENTER  CHARACTER  TO  CONTINUE. 
READY  (Y.N)  ?  ....I  Y 


L2 


M.ilN 


::  LIST 

LISTING  OF  DATA  FILES  AUAILABLE 

•  10 - REF  INI 

til - PLATE 

f  12  —  REFXN2 

it  DELE 

ENTER  LABEL  jjREFINI 

it  DELE 

ENTER  LABEL  j jPLAT2 

it  REFINE 

ENTER  INPUT  FILE  tl  LABEL  i j REF INS 

ENTER  OUTPUT  FILE  LABEL  | jREFINS 

CREATING  NEU  DATA  FILE 

BEGIN  REFINE  ROUTINE... 

1  -  RANDON  ELENENT  INPUT 

2  -  RANGE  OF  ELEMENTS 

3  -  EXECUTE  REFINEMENT 

ENTER  OPTION - >8 

MAXIMUM  NUMBER  OF  ELEMENTS 

SPECIFIABLE - >  15* 

NUMBER  STILL  AVAILABLE  - >  1SB 

ENTER  BEGINNING  ELEMENT  NUMBER.  ENDING  ELEMENT 
NUMBER,  AND  INCREMENT  FOR  RANGE. 

ENTER  NBEG.NEND,  I  NCR - >1.5,1 

1  -  RANDOM  ELEMENT  INPUT 
8  -  RANGE  OF  ELEMENTS 


->V 

0PT8,  HELP)  ->0PT8 


it  PLOT 

ENTER  INPUT  FILE  tl  LABEL  jjREFINS 

ttM 

INVALID  COMMAND.  ENTER  'HELP'  FOR  LIST  OF  COMMANDS 
tDRAU 

BTEMQF  PLOT.  ENTER  CHARACTER  TO  CONTINUE. 

ii"A2>  (  I  »  ...  I  V 


3  -  EXECUTE  REFINEMENT 
ENTER  OPTION 


->3 


ENTER  CUT  DIRECTION  (R.S.T)  >R 


MAXIMUM  CUTS  PER  ELEMENT 
ENTER  CUTS  PER  ELEMENT  — 


->  S 
->4 


EQUALLY  SPACED  CUTS  ?  (Y»N>  >Y 

REFINE  COMPLETE. 

PERFORM  SUBSEQUENT  TIDY  ?  (Y.N) 

TIDY  OPTION  ONE  OR  TUO  ?  (OPTl. 

TIDY  COMPLETE 
SD  DUPLICATE  NODES  DELETED 
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LISTING  OF  DATA  FUES  AVAILABLE 


•  10 - REFIN3 

•  12 - REFIN2 

««  DELE 

ENTER  LABEL  j  >REFIN2 

tt  RENUN 

ENTER  INPUT  FILE  SI  LABEL  j  jREFIN3 

ENTER  OUTPUT  FILE  LABEL  } j PLATE —4 

CREATING  NEU  DATA  FILE 

NODE  POINT  REORDERING  COHPLETE 

NODAL  BANDWIDTH  (OLD)  •  9? 

NODAL  BANDWIDTH  (NEU)  •  9? 


II  DELE 

ENTER  LABEL  jjR£FlN3 

if  LIST 

LISTING  OF  DATA  FILES  AVAILABLE 
•11  -  PLATE-4 


it  BOUNDS 

ENTER  INPUT  FILE  SI  LABEL 
ENTER  OUTPUT  FILE  LABEL 
CREATING  NEU  DATA  FILE 


> I PLATE -4 
j j PLATf+BCS 


BOUNDARY  CONDITIONS  NAY  BE  ENTERED  FOR  RANDOM 
NODES.  RANGES  OF  NODES,  ALL  NODES  ON  A  PLANE. 
OR  THE  ENTIRE  MODEL. 


FOR  EACH  GROUP  OF  NODES,  CONSTRAINTS  HAY  BE  SET 
IN  ANY  COMBINATION  OF  THE  X.  V,  t  Z  DIRECTIONS. 
(NOTE  THAT  2-D  MODELS  MUST  BE  FIXED  IN  Z) 


CONSTRAINT  DIRECTION  CODE  CONSTRUCTED  AS  FOLLOWS! 
ENTRY  1  - 
ENTRY  R  - 


'1' 

'O' 


FOR  CONSTRAINED  IN  X  DIRECTION 
FOR  UNCONSTRAINED  IN  X  DIRECTION 


SAME  FOR  Y 
SAME  FOR  Z 


->  'l.S.i' 


FOR  CONSTRAINTS  IN  X  AND  Z 
DIRECTIONS  BUT  NOT  IN  V. 


1  -  RANDOM  NODES 
R  -  RANGE  OF  NODES 

3  -  SPECIFIED  PLANE 

4  -  ALL  NODES 

5  -  EXIT 


ENTER  NODS  SELECTION  OPTION  <1.8.. ..5) 


->4 


•  i  e 


rt‘  -0ssT  *n  c:re>:  ion  c<  :e  mp-  f.  v  iue:;  - 

1  RANCOR  node; 
a  RANGE  OF  NCiES 
3  -  SPECIFIED  PLANE 
A  -  ALL  NODES 
S  -  EXIT 

ENTER  NODE  SELECTION  OPTION  (1.3... .5)  - >3 

ENTER  CONSTRAINT  DIRECTION  CODE  (THREE  VALUES )  ->1  ( 

- >1  9 


PLANE  IS  DEFINED  BY  AX  ♦  BY  ♦  CZ  •  D. 
ENTER  COEFFICIENTS  (A.B.C.D) 


->3 


34  NODES  FOUND 

1  -  RANDOM  NODES 
3  -  RANGE  Of  NODES 

3  -  SPECIFIED  PLANE 

4  -  ALL  NODES 

5  -  EXIT 

ENTER  NODE  SELECTION  OPTION  (1.3.... S)  - 

ENTER  CONSTRAINT  DIRECTION  CODE  (THREE  VALUES)  ->•  1 

PLANE  IS  DEFINED  BY  AX  ♦  BY  ♦  C2  "  D. 

ENTER  COEFFICIENTS  (A.B.C.D)  - 

14  NODES  FOUND 

1  -  RANDOM  NODES 
3  -  RANGE  OF  NODES 

3  -  SPECIFIED  PLANE 

4  -  ALL  NODES 

5  -  EXIT 

ENTER  NODE  SELECTION  OPTION  (1.3.. ..S)  - 


->•  1 


~>5 


»  BOUNDARY  CONDITIONS  ADDED 
3D  BOUNDARY  CONDITIONS  TOTAL 


n  LIST 

LISTING  OF  DATA  FILES  AVAILABLE 

»}•  -  PLATE  ♦BC3 

•11  -  PLATE-4 

II  dele 
enter  UttEL 


jj  PLATE -4 

ll  PLOT 

EHTKyINPUT  FILE  «i  LABEL  ijPLATE+DCS 

UNION  AXIS  IS  VERTICAL? 

ENTER  1  FOR  X.  3  FOR  V,  OR  3  FOR  Z . I  t 

•LYE 

ENTER  THE  EVE  POSITION, 

XEVE.  VEYE.  ZEYE . . . « 

•DRAG  * 

a¥*Kro  n 


•  • 


•  • 


1NM 
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V 


30  VGU  WAIT  n,,;  2r,;,n  FU  CT  ;<»*  (V.H  .  ?  . s  • 

-'-ALE  THE  LOON  AREA  (V.N)  ■> . .  .• 

DIGITIZE  THE  LOWER  LEFT  CORNER 
RNO-TMC-UFRER  RICmt  CORNER  of  the 
»OON  AREA 

RCjWV?<V.N>  . .  y 

•U  IE 

.'BEL  'HE  EU  lENTSXV.i)  . .  . 

U1EL  THE  NOI-:S?tV.N)  . I  .......t  v 

^J-^BEL  ALL  THE  SURF  ACES  ?<V,N) . .....I  •, 

tORAW^  W,fU  186  **E5,<Y'*>  . .  N 


SLSffS^n^r*  ENTt5  CH^CTER  TO  CONTINUE 

RCAOV  <V#H)  ?  m»«I  V 


j|PLATE*BCS 
I I FINAL 


DC  VOU  UAff‘  THE  ZOOM  FUNCTION  (V,N>  ?  .  N 

♦.'AIM 


it  LIST 

LISTING  OF  DATA  FILES  AVAILABLE 

•  10 - PLATE+BCS 

ll  LOAD 

ENTER  INPUT  FILE  B1  LABEL  jjPL 

ENTER  OUTPUT  FILE  LABEL  I I FI 

CREATING  NEU  DATA  FILE 

LOAD  SPECIFICATION  OPTIONS 
N  -  NODAL  LOAD  ENTRY 
E  -  ELEMENT  LOAD  ENTRY 
L  -  LIST  EXISTING  LOADS 
H  -  PRINTS  THIS  LIST 
S  -  STOP  LOAD  ENTRY 

ENTER  LOAD  SPECIFICATION  OPTION  - >N 

NODE  SPECIFICATION  OPTIONS 
S  -  SINGLE  NODE 
R  -  RANGE  OF  NODES 
P  -  ALL  NODES  ON  A  GIVEN  PLANE 
H  -  PRINTS  THIS  LIST 

E  -  EXIT  NODAL  LOAD  SPECIFICATION  SECTION 


ENTER  NODE  SPECIFICATION  OPTION  - >S 

ENTER  CASE  NUMBER  >1 

ENTER  FORCE  VECTOR  (FX.FY.FZ)  - >•  100  • 

NODE  NUNBER  TO  BE  LOADED  ?  - >3M 

ENTER  NODE  SPECIFICATION  OPTION  - >E 

LEAVING  NODAL  LOAD  SPECIFICATION  SECTION 

ENTER  LOAD  SPECIFICATION  OPTION  - >H 

LOAD  SPECIFICATION  OPTIONS 
N  -  NODAL  LOAD  ENTRY 
E  - 
L  - 
H  r 

S  -  STOP  LOAD  ENTRY 

ENTER  LOAD  SPECIFICATION  OPTION  - >S 

LEAVING  LOAD  SPECIFICATION  NODULE 
n  LIST 

LISTING  OF  DATA  FILES  AVAILABLE 

BIO  -  PLATE+BCS 

BU - FINAL 

ll  DELE 

ENTER  LABEL  J  »PLATE+BCS 


««  STOP 

LISTING  OF  DATA  FILES  AVAILABLE 


0S37S O  MAXIMUM  EXECUTION  FL. 

19. SOI  CP  SECONDS  EXECUTION  TIME. 

FILE  QUOTA  EXCEEDED 
COMMAND-  FILES 
—LOCAL  FILES — 

TAPE 14  SP  SIMPUT  SOUTPUT  TAPE1S 

TAPE1B  TAPE13  TAPE11  TAPE1S  TAPE1S 
TAPE! 7 

COMMAND-  RETURN, TAPE10»TAPEia,TAPEl3,TAPS14,TAPElS,TAPElS,TAPEl? 


MCSM  MN2MM CXECUTION  Ft. 

cornua*  *r?i  ft  stC0NM  b*cution  n«. 

LCCAL  I  £S-  - 

*  iscin  *IHPUT  i* output  jni  rr 


FDATA 


C0NHAM9-  RETURN, JNf AT,  ISCSI 
CONNAH0-  REWIND. FDATA 
COWHAND-  REQUEST. DATA. SPF. 

COWHAND-  COPY *F, FDATA, DATA 

COWHAND-  CATALOG. DATA, SHEETUITHHOl£,CY*l»RP*QA#, IDsBROCKHAN 
INITIAL  CATALOG 

CT  IS-  FROCWIAH  PFN»SHEETWITHHOLE 
CT  CY*  HI  SH-AFFDL  MMAACfl  1 1  WORDS. 
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7.2  EXECUTION  OF  CORGEN  /  EXPAND  /  PREP  /  REFMT  SEQUENCE 

This  Subsection  contains  an  example  of  the  use  of  CORGEN  to 
generate  a  finite  element  mesh  from  coordinate  data  in  lofting 
form,  which  in  this  instance  is  stored  on  a  data  file.  The 
EXPAND  utility  is  used  to  translate  the  surface  model  output  by 
CORGEN  to  three-dimensional  form.  PREP  is  used  to  refine  the 
mesh  and  specify  constraints.  REFMT  generates  the  formatted 
data  deck  ready  for  analysis. 

The  geometric  data  used  in  preparing  this  model  is  lofting 
data  describing  the  forward  canopy  of  the  T-38  aircraft.  Near 
the  forward  edge,  the  canopy  rests  on  a  metallic  arch  structure 
which  does  not  provide  for  any  bolted  attachment;  the  final 
model  will  contain  a  rather  small  row  of  elements  along  this 
boundary,  which  could  be  used  to  generate  contact  elements  for 
detailed  analysis  of  the  true  support  condition. 

7.2.1  Summary  of  Modelling  Procedure 

The  modelling  data  provided  for  the  T-38  airplane  canopy 
consists  of  nine  lofting  planes  with  seven  points  defined  on 
each  plane.  This  data  was  entered  to  CORGEN  as  a  file  struc¬ 
tured  so  that  each  lofting  line  corresponded  to  one  group  of 
data.  The  groups  were  input  such  that  the  first  group  repre¬ 
sented  the  first  lofting  line  at  the  fore  point  on  the  canopy 
and  the  ninth  group  represented  the  aft  area  of  the  canopy. 


The  CORGEN  program  was  executed  to  convert  the  lofting 
data  into  the  preprocessor  internal  format  description.  The 
program  execution  was  uninvolved  for  this  model  as  the 
default  options  for  group  formulation  and  element  generation 
were  utilized.  The  sequence  of  commands  to  create  the  model 
are  identified  below: 

1  -  NODE?  FILE?  READ?  EXIT  -  node  input  from  file, 

2  -  ELEM?  GENE?  EXIT  -  element  generation, 

3  -  FINI  -  write  model  data  to  file. 

These  operations  resulted  in  a  2-D  surface  model  consisting  of 
twelve  nine  node  elements. 

The  EXPAND  program  was  executed  to  translate  the  surface 
model  to  a  three-dimensional  solid  element  model.  There  exists 
only  the  option  of  selecting  the  direction  to  expand  the  surface 
in  this  program. 

The  PREP  program  was  executed  to  perform  several  functions: 

-  verification  of  the  lofting  model  data, 

-  mask  to  16  node  elements  (from  27  node  elements), 

-  refine  the  model, 

-  optimize  the  node  numbering  for  the  analysis, 

-  select  14  point  integration  rule  for  analysis,  and 

-  specify  boundary  conditions. 

The  sequence  of  steps  to  execute  PREP  are  outlined  below  and 
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listed  in  the  following  paragraphs  of  this  section. 


STEP  1 


STEP  2 


STEP  3 


STEP  4 


STEP  5 


STEP  6s 


:  Verification  of  CORGEN/EXPAND  model 

-  PLOT  (plot  original  canopy  model) 

-  PLOT  (find  local  element  orientation  axes) 

s  Create  16-node  solid  elements  from  27  node  model 

-  MASK  (eliminate  extraneous  nodes) 

s  Refine  the  model 

-  REFINE  (execute  2  times  for  R  direction) 

-  PLOT  (verify  refinement,  check  new  element 

numbering  for  next  refinement) 

-  REFINE  (only  selected  elements  in  S  direction) 


Optimize  node  numbering 

-  RENUMBER  (repeated  to  insure  optimal  numbering) 

-  PLOT  (verify  model  geometry,  and  identify  node 

numbers  for  boundary  conditions) 


Select  element  properties 

-  PROPS  (select  14  point  integration  rule, 

remainder  of  properties  given  in  REFMT ) 

Specify  boundary  conditions 

-  PLOT  (ZOOM  to  find  nodes  to  be  constrained) 

-  BOUNDS  (constrain  nodal  degrees  of  freedom) 


A  natural  frequencies  analysis  will  be  performed,  so  no  loads 
data  are  specified  at  this  time. 


The  REFMT  program  was  executed  to  prepare  a  data  file 
for  MAGNA  analysis.  This  section  details  data  specification 
for  a  natural  frequency  analysis  of  the  model.  The  REFMT 
program  has  interpreted  the  16  node  elements  as  MAGNA  type  8 
element  of  which  there  are  45  elements.  The  sample  run 
illustrates  the  definition  of  undefined  material  property 
types.  A  final  file  is  generated  for  a  MAGNA  analysis. 

Sections  7.2.2  through  7.2.5  list  the  program  executions 
of  the  CORGEN  /  EXPAND  /  PREP  /REFMT  sequence  for  completing 
the  T-38  aircraft  windshield  model. 
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CORGEN  Execution 


CONNAHD-  ATTACH.P.PREPROCESSORPROC,  ID*  WOCICnAH,f1R*l 
AT  CV.  993  SN'AFFDL 

COfWrtMD-  ATTACH, TAPCU  ,  UINDSHIELDDATA,  CY*  1 . 

CONHAND-  FILES 
—LOCAL  FILES— 

•INPUT  •OUTPUT  «P 
COMMAND-  KGlH.COftGEM.P. 


STARCH 


7.2.3  EXPAND  Execution 


CONNAHD-  FILES 
—LOCAL  FILES— 

•INPUT  SOUTPUT  » 
CONNAND-  ICG IN, EXPAND. P. 


TAPEIO 


•COIN  EXPAND 

THICKNESS  EXPANSION  PREPROCESSOR 


THICKNESS  EXPANSION  OPTIONS! 
t  -  EXPAND  FROH  SOTTON  -  UP 
8  -  EXPAND  ADOUT  NIDPLANE 
3  -  EXPAND  FRON  TOP  -  DOWN 
ENTER  OPTION  (1.8.3) . t  3 


DO  VO U  WISH  A  LISTING  OP 
EXPANSION  DATA  (V,N)?...<  N 
STOP 

034700  NAXIHUH  EXECUTION  FL. 

0.357  CP  SECONOS  EXECUTION  TINE. 
COHNAND-  FILES 
—LOCAL  FILES— 

TAPE12  SINPUT  SOUTPUT  SP 
TAPE 11 

RETURN.  TAPE10.TAPE18 


TAPE10 
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7.2.4  PREP  Execution 


COWMNO-  FILES 
--LOCAL  FILES" 

•INPUT  SOUTPUT  tP 
COMMHD-  BEGIN, PREP. P 


TAPClt 


tnnmraummituunmmmtm 
muwmtl  BEGIN  PREP  »umuu« 
tttt  HAGHA  PREPROCESSING  UTILITIES  S»M 


4ui 


■*■<**» 


JoitJ 


it  LIST 

LIST! NO  OF  DATA  HUS  AUAXUWU 

lit  -  ACFIN3 

Ml  -  RCHUNR4 

M3  -  KEFIRS 

it  SCUTE 

ENTER  UWCL 

j  j REFINE 

it  SCUTE 

ENTER  UWCL 

,,RCFIN3 

if  RCMUNICR 

ENTER  INTUT  FIU  M  UWCL 
ENTER  OUTPUT  FIU  UWCL 
CREATIMO  NCU  DATA  FIU 

jjRENUWM 

llRCNURRS 

NOOC  POINT  RCORBCRXNtt  CONFUTE 

MORAL  BAHBUIOTH  (OLS)  • 
MORAL  RANWIIPTII  (NCU)  ■ 

TO 

TO 

ft  SCUTE 

INTER  UWCL 

llACNUHBR 

t*  PLOT 

Bfgg^IHPUT  FIU  M  UWCL 

f  jRCNUNM 

AT  CNR  OF  PLOT.  OTTER  CHARACTER  TO  CONTINUE. 
MART  (V.N)  . . . 

Wm 
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KAKl 

LABEL  THE  EL£NENTS7<Y,N>  . . 

LABEL  THE  NODES»(V.N>  . . . V 

LABEL  ALL  THE  SURFACES?! Y,N>  . «  N 

ENTER  THE  SURFACE  NO.  <0-0,  #  PLOTS  ALL)  ....I  8 
•HELP 

COWHAND  DESCRIPTION 

it»m  ttntnttst 

AXES  AXES  DRAM  AND  LABEL 

BOUNDS  PLOT  BOUNDARY  CONDITIONS 

CLIP  CLIP  PLANE  POSITION 

CUBE  SET  NINIHA  AND  MAXIMA 

DEFAULT  SET  MFAULT  VALUES 

DRAW  DRAU  MODEL 

ELEMENTS  PLOTS  ALL  OR  SELECTED  ELEMENTS 

EYE  EYE  POSITION 

HELP  LIST  ALL  COMMANDS 

LABELS  LABELS  ELEMENTS  AND/OR  NODES 

MAIN  RETURNS  CONTROL  TO  MAIN  PROGRAM 

ORIENT  PLOTS  R.S.T  ORIENTATION  AXES 

PROJECTION  PROJECTION  TYPE 

REFLECT  REFLECT  A  PLANE 

ROTATE  ROTATE  MODEL  ABOUT  AXES 

SCALE  SCALE  PLOT 

SHRINK  SHRINK  ELEMENTS 

SUMMARY  LIST  AU.  PARAMETER  VALUES 

TIME  PRINT  CPU  TIME  SINCE  START  OF  SESSION 

TRANSLATE  TRANSLATE  MODEL  FROM  ORIGIN 

VERTICAL  VERTICAL  AXIS 

ZOOM  ZOOM  ON  THE  MODEL 

SITE  SITE  POSITION 

SEYE 

ENTER  THE  EYE  POSITION. 

XEVE.  VEVE.  ZEVE .  ISO  200  ISO 

BELEM 

PLOT  ALL  THE  ELEHENTSTtV.N)  . . 

ELEMENT  INPUT  OPTIONS 

OPTION  11  RANDOM  ELEMENTS (4. 13. 64....  ) 

WHERE  FIRST  DIGIT  •  NO.  ELEMENTS 
OPTION  3»  RANGE  OF  ELEMENTSU.8.8  FOR  1-9. INC  OF  3) 

ENTER  THE  INPUT  0PTI0NU.3)  . «  1 

ENTER  THE  ELEMENTS . *  S  7,34,36,40.41 


ELEMENTS 

EYE 

HELP 

LABELS 

MAIN 

ORIENT 

PROJECTION 

REFLECT 

ROTATE 

SCALE 

SHRINK 


TRANSLATE 

VERTICAL 


NT  END  OP  PLOT,  ENTER  CHARACTER  TO  CONTINUE 
HEADY  <V,N!  T  .  V 


DO  YOU  UANT  TV*  ZOOM  FUNCTION  (V.N) 

SCAIC  THC  ZOON  MICA  (V.N)  ? . 

OIQITIZZ  THC  LOUOt  LEFT  CORNC* 

AND  THC  UFFC*  RIGHT  COANCR  OF  THC 
ZOON  ARCA 

RCAOV?(Y.N)  . 


SMAU 


AT  DO  OF  W.OT,  CMTCR  CMAWACTCH  TO  CONTXNUC 
RCAOV  (V.N)  ?  ....<  V 


t  N 
.1  N 


DO  YOU  WANT  THE  ZOOM  FUNCTION  (V.N)  ?  . 

CUN 

PLOT  ALL  THE  EUNEHTS?<Y,N)  . 

CLCnCNT  INPUT  OPTIONS 

OPTION  It  RANDON  ELENENTS(4.13,S4, . . .  )  _ 

UHERE  FIRST  DIGIT  -  NO.  ELENEHTS _ 

OPTION  at  RANGE  OF  EUHENTSU.9.8  FOR  1-9. INC  OF  8) 

ENTER  THE  INPUT  0PTI0NC1.8)  . '  1  _  _ 

ENTER  THE  CLEHENTS  . . . •  *  7.34.35.40.41,48 

LAK 

LASCL  THE  ELENEHTS?<Y.N)  . . »  N 

LABEL  THE  N0DES7<V,N)  .  Y 

LABEL  ALL  THE  SURFACES? (Y.N)  . I  N 

ENTER  THE  SURFACE  NO.  CO-S.  S  PLOTS  ALU  ....t  3 


AT  END  OF  PLOT,  ENTER  CHARACTER  TO  CONTINUE. 
READY  (V,N>  . . . 


i  v 
t  N 


00  VOU  UAHT  THE  200f1  FUNCTION  CY.N)  ? . 

SCALE  THE  ZOON  AREA  (V.N)  ?  . 

DIGITIZE  THE  LOWER  LEFT  CORNER 
AND  THE  UPPER  RIGHT  CORNER  OF  THE 
200N  AREA 

READY?! V.N >  . t  V 


(DRAW 

AT  END  OF  PLOT,  ENTER  CHARACTER  TO  CONTINUE. 
READV  (V.N)  ?  ....I  V 


c‘; 


•  rLEM 

PLOT  ALL  THE  ELEMENTS?(Y.N>  .  N 

ELEMENT  INPUT  OPTIONS 

OPTION  II  RANDOM  ELEMENTS (4. 13, 64, . . . ) 

WHERE  FIRST  DIOIT  •  NO.  ELEMENTS 
OPTION  St  RANGE  OF  ELEMENTS! 1,8,3  FOR  1-9, INC  OF  8) 


ENTER  THE  INPUT  OPTION<1.2>  . i  8 

ENTER  THE  ELEMENTS . i  1,9.1 

•ZOOM 

DO  VOU  WANT  THE  ZOOM  FUNCTION  <V,N>  ? . i  N 

SORMU 


AT  END  OF  PLOT.  ENTER  CHARACTER  TO  CONTINUE. 
READV  <V,N)  . . V 


t: 

« 
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•  . . -  - 


-•  -  *  -  •  - 


ll  STOP 

LISTING  OF  DPT*  FILES  AMAILAiLE 


•IS - - 

•11 - R£NUH»4 

tia - BOUNDS 


STOP 

NNN  NAXINUN  EXECUTION  FI. 

17. 031  CP  SECONDS  EXECUTION  TI1C. 

FILE  QUOTA  EXCEEDED 
COWHAND-  FILES 
—LOCAL  FILES— 

•INPUT  SOUTPUT  SP  TAPE11  TAPE!# 

TAPE I 4  TAPE13  TAPE13  TAPE17  TAPE1S 

TAPC1S  TAPE1S 
FILE  OUOTA  EXCEEDED 

COHHANO-  RETURN, TAPE 1*. TAPE 11, TAPE 13, TAPE 14, TAPE IS 
CONNAND-  RETURN.TAPE18.TAPE17.TAPE1* 

CONHAUD-  FILES 
—LOCAL  FILES— 

•INPUT  •OUTPUT  SP  TAPE18 

COWHAND-  REWIND. TAPC12 
COWHAND-  C0PYDF,TAPE13.UNFHT 
EOI  ENCOUNTEMS  AFTER  COPY  OF  FILE 


•.  RECORD  1 
RETURN. TAPEia 
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7.2.5  REFMT  Execution 


command-  PILES 
--LOOAL  PILES— 

•INPUT  BOUTPUT  *P 
COMMAND-  BEGIN, PEP NT, P. 


UNFNT 


BEGIN  REFNT  —  MAGNA  INPUT  GENERATOR 


«  SUMMARY  OP  INITIAL  SCAN  OP  DATA  PILE  XX 


OP  NODAL  POINTS . 

_  OP  ELEMENTS  (TOTAL)  . 

NUMBER  OP  CONSTRAINT  RECORDS 
NUMBER  OP  LINEAR  CONSTRAINTS 

NUMBER  OF  NODAL  LOADS  . 

NUMBER  OP  ELEMENT  LOADS  .... 
DISTINCT  LOAD  CASES  /  GROUPS 


33S 

45 

4# 

B 

A 

B 

B 


ELEMENT  TYPE 


NUMBER  OP 
ELEMENTS 


ELEMENTS  WITH 
UNSPEC.  HATRL 


ENTER  A  TMREE-L1ME  PROBLEM  TITLE  (UP  TO  SB  CHARACTERS  PER  LINE) 

. CANy|syR^IC,TOa,^siw,»UPPOT‘-iiGENU^>SOLufi0Mt 

GENERATED  IE  JAM.  18E3 


MAJOR  SOLUTION  OPTIONS  A NO  PARAMETERS 


ENTER  ANALYSIS  TYPE  (UNBAR, NONLINEAR)  -LINEAR 

ENTER  ANALYSIS  SUBTYPE  < STATIC. DYNAMIC ) - DYNAMIC 

DYNAMIC  SOLUTION  TYPE  (EIQEN  .TRANSIENT)  -EIQEN 
INTER  NUMBER  OP  FREQUENCIES  TO  IE  COMPUTED  -4 
CONSISTENT  NASS  MATRICES  TO  BE  USED  (Y/N)  —V 
POSTPROCESSOR  PILE  TO  BE  WRITTEN  (Y/N)  - Y 


OR)  OP  OPTIONS  SPECIFICATIONS 


INDIVIDUAL  ELEMENTS  IN  THE  NOBEL  CONTAINUNDEPINED  PROPERTIES 

ELEMENT  TYPE  S 

NUMBER  OP  ELEMENTS*  4S 

PLEASE  DEFINE  A  DEFAULT  PROPERTY  COSE  FOR  THIS  ELEMENT  TYPE, 
OR  ENTER  MATERIALS  DATA  OIRECTLY  BELOU. 
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L  ■  (DIST  CURRENT  PROPERTIES  TABLE 
C  •  <C)MANG£  AN  ENTRY  XN  THE  TABU 
S  •  (S)TOP  EDITING 

ENTER  OPTION  I  l  ,  C  ,  S  )  - 


CONNAND-  CATALOG, DATA, UINDSHIELDDATAEIGEN* ID*BROCKNAN,RP«B#G 
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7.3  EXECUTION  OF  IJKGEN  /  PREP  /  REFMT  SEQUENCE 


This  Subsection  contains  a  preprocessing  example  demon¬ 
strating  the  use  of  the  default  options  provided  in  IJKGEN 
(Subsection  2.2)  to  generate  a  three-dimensional  model  of  a 
folded  plate-type  structure.  The  use  of  user-written  sub¬ 
routines  within  IJKGEN  is  described  separately  in  Appendix  D. 

The  component  to  be  modelled  in  this  hypothetical  example 
involves  a  folded  plate  containing  both  flat  and  cylindrical 
geometry.  The  plate  is  twenty  inches  wide,  with  symmetry 
assumed  along  its  centerline.  Its  folded  edges  are  repre¬ 
sented  by  90-degree  cylindrical  segments  of  mean  radius  5.1 
inches,  and  ten  inches  long  (half  the  width  of  the  plate) . 

A  flat  segment  between  the  folded  edges  is  40  inches  in  length, 
and  symmetry  is  assumed  along  the  center  of  this  panel  as  well. 
The  thickness  of  the  entire  panel  is  uniform,  at  0.20  inches. 

7.3.1  Summary  of  Modelling  Procedure 

The  basic  geometric  components  of  the  plate  are  the  flat 
section  and  its  cylindrical  ends.  These  are  generated  in  IJKGEN 
separately,  using  the  default  coordinate  system  and  mesh  grading 
options : 

-  Rectangular  segment  s  Cartesian  coordinates,  graded  mesh 

-  Cylindrical  segment  :  Cylindrical  coordinates,  uniform  mesh 


The  rectangular  (flat)  portion  of  the  mesh  has  been  refined 
slightly  near  the  symmetric  center  section,  where  a  band  of 
pressure  will  be  applied  in  the  analysis.  The  row  of  elements 
to  be  joined  to  the  cylindrical  end  has  been  made  four  times  as 
long  as  the  row  at  the  center,  and  the  grid  is  uniform  in  the 
remaining  two  directions. 

The  two  substructures  created  in  IJKGEN  are  assembled  with 
the  use  of  PREP,  after  using  the  TRANslate  and  ROTAte  operations 
to  position  them  properly.  Observe  that  the  flat  section  of  the 
plate  could  have  been  positioned  properly  in  IJKGEN  by  an 
appropriate  selection  of  the  minima  and  maxima;  however,  the 
repositioning  has  been  performed  here  to  demonstrate  more  of  the 
features  of  PREP. 

The  sequence  of  operations  performed  in  PREP  is  summarized 
below  for  easy  reference: 

Step  1.  Verification  of  I JKGEN-Generated  Models 

-  PLOT  (Plot  original  cylinder  data) 

-  MASK  (Convert  cylinder  to  16-node  elements) 

-  PLOT  (Plot  original  plate  data) 

-  MASK  (Convert  plate  to  16-node  elements) 

Step  2.  Repositioning  of  Plate  Model  Segment 

-  ROTAte  (Two  90-degree  rotations  of  plate  model) 


-  PLOT  (Verify  that  rotations  are  correct) 

-  TRANslate  (Move  plate  next  to  cylindrical  segment) 

-  PLOT  (Verify  that  plate  is  positioned) 

Step  3.  Combination  of  Flat  Plate  and  Cylindrical  Sections 

(Combine  the  two  models) 

(Verify  that  it  has  been  done  correctly) 
Step  4.  Finish  Geometric  Model 

-  RENUmber  (Reorder  nodes  to  reduce  bandwidth) 

-  RENUmber  (Verify  that  a  second  pass  not  needed) 


-  MERGe 

-  PLOT 


Step  5.  Apply  Boundary  Conditions  and  Loads 


BOUNds 


-  PLOT 


-  LOADS 


(Apply  symmetry  conditions  on  the  edges 
Z=0  and  Y=-10»;  clamp  the  free  end  of 
the  cylindrical  segment,  X=0.) 

,(Use  ORIENT  plot  option  to  determine 
the  correct  surface  numbers  for 
pressures) 

(Apply  inward  pressure  on  the  upper 
faces  of  elements  41  through  50) 


Step  6.  Finish  Model  by  Assigning  Properties 


-  PROPerty 


(Assign  all  elements  as  solids,  to  be 
integrated  using  a  14-point  quadrature 


LIST 


rule;  leave  the  specification  of 
materials  data  for  the  REFMT  step) 
(Make  sure  of  the  local  file  number 
which  contains  the  completed  model; 
in  this  case,  local  file  TAPE12  must 
be  saved) 

-  STOP  (Leave  PREP) 


The  above  sequence  of  operations  is  typical  for  modelling 
data  generated  with  IJKGEN,  in  which  the  original  model  can  be 
generated  with  the  necessary  degree  of  refinement.  It  is  worth 
noting  that  the  specification  of  properties,  which  is  quite 
simple  in  this  case,  can  become  more  tedious  when  several 
materials  are  involved  and  a  number  of  REFIne  operations  are 
needed.  In  such  cases,  it  is  a  good  practice  to  assign  the 
properties  at  an  early  stage,  where  the  element  numbering 
sequence  is  generally  more  regular. 

The  final  phase  of  model  generation  is  the  REFMT  step,  in 
which  a  complete  data  deck  for  the  MAGNA  program  is  generated 
in  this  example.  The  initial  output  from  REFMT  indicates  that 
the  data  has  been  read  correctly,  and  that  REFMT  has  properly 
interpreted  the  16-node  elements  as  being  Type  8  elements  in 
MAGNA.  The  following  options  have  been  selected  through 
additional  user  inputs: 


-  nonlinear,  static  solution; 

-  10-increment  solution,  using  a  constant  time  step; 

-  combined  full  and  modified  Newton-Raphson  iterations; 

-  printing  at  every  other  increment,  and  writing  of  a  post¬ 
processor  file  at  every  five  increments; 

-  restart  file  to  be  written  at  every  other  increment,  with 
the  restart  file  label  FPLT. 

The  next  segment  of  REFMT  execution  deals  with  material 
properties,  since  none  were  specified  within  PREP.  Since  al 
elements  in  the  model  contain,  the.  same  material  code  (zero), 
the  user  is  requested  to  define  a  default  material  code  for 
all  Type  8  elements  having  this  material  code.  Entries  for 
carbon  steels  are  listed  from  the  materials  library,  and  a 
default  material  code  of  120  (G10500  steel)  is  selected.  No 
editing  of  the  materials  data  is  performed,  although  entries 
may  be  altered  selectively  as  many  times  as  necessary. 

Finally,  information  is  requested  concerning  the  variation 
of  the  applied  loads  for  the  nonlinear  solution.  Since  the 
solution  is  static,  a  ramp  function  is  selected  (option  3), 
causing  the  loads  to  applied  incrementally  in  equal  steps. 

The  pressure  magnitude  input  in  PREP  is  retained  as  the  peak 
loading  value,  and  all  pressures  are  defined  as  being  'live 
loads' ,  to  be  updated  according  to  the  current  surface  orien¬ 
tation  continuously  during  the  nonlinear  ana .ysis. 


Once  the  data  generation  is  complete,  the  local  file  FDATA, 
which  contains  the  complete  input  deck  for  MAGNA,  is  cataloged 
as  a  permanent  disk  file  with  the  name  FOLDEDPLATEDATA.  All 
that  remains  to  perform  the  nonlinear  analysis  is  the  submission 
of  the  proper  job  control  stream  for  MAGNA,  as  described  in 
Reference  [4]. 

The  complete  modelling  session,  including  all  printed  output 
and  plotting,  is  reproduced  in  Paragraphs  7.3.2  through  7.3.4 
below. 


IJKGEN  Execution 


COWI*®-  ATTACH.P,PRCPR0CCSS0RPR0C.ID>BR0CKMAH.NR>1. 

AT  CY-  90S  SH-AFFDL 
COMMAND-  BCGIN.IJKGEN.P. 


ninniiiHHitmtimnmmnHiiii 
numm-  begin  ijkcen 


OP  GEOMETRIC  NCSH  DATA  FOR  SOHO. 
rnnSniSfrS Ei252I  USING  AN  INTEGER  - 

ta\  )  -  DEFINE  NCSH  GEONETRV 

tiFnSrrniSMZXi: V»g >  *  COORD.  THANSFORHATIOH 
.ySiJS*  RECTANGULAR,  CYLINDRICAL  OR  SPHERICAL 
COORDINATES.  AND  UNIFORM  OR  PROPORTIONALLY  GRADED  NESH  SPACING 


(3) 


uimir  -  USER  PARAMETER  INPUT  ROUTINE  (INITIALIZE  data 
IN  BLANK  CANNON) 


USER  SUBROUTINE  'SURF AC'  NOT  QIUCH 


BUILT-IN  NESH  DIUISION  OPTIONS  ARE  AS  FOLLOWS  - 
<1 >  “  UNIFORM  NESH  IN  EACH  DIRECTION  FOfcU*" 

(R)  -  GRADED  NESH  (SPECIFY  RATIO  OF  FIRST/LAST  ELEMENT  SIZE 


ENTER  OPTION  (1,2) . . . g 

ENTER  THE  RATIO  OF  FIRST  /  LAST  ELEMENT  LENGTHS  FOR  EACH 
COORDINATE  DIRECTION  (ALPHA.  BETA.  ZETA)  (R-l  FOR  UNIFORM) 


ENTER  LENGTH  RATIOS  (R1.R2.R3) . 4..  1..  1. 

USER  SUBROUTINE  'CRDTRN'  NOT  GIVEN 


BUILT-IN  COORDINATE  SYSTEM  TRANSFORMATION  OPTIONS  ARC  - 

(l)  RECTANGULAR.  (2)  CYLINDRICAL.  (3)  SPHERICAL 


ENTER  COORDINATE  SYSTEM  OPTION  (1.2.3)  -1 

SS  PLEASE  NOTE  THE  FOLLOWING  CONVENTIONS  FOR  RECTANGULAR  SYSTEM  XX 

ZETA'  •  Z 


'ALPHA'  -  X  'BETA'  •  Y 
A  RIGHT-HANDED  SYSTEM  IS 


ENTER  UNITING  SURFACE  COORDINATE  VALUES  - 


1.  ALPNA(NIN) 

3.  BETA  (HIM) 
S.  ZETA  (HIM) 


2.  ALPHA (MAX) 

4.  BETA  (MAX) 

5.  ZETA  (MAX) 


S.S.  20.0,  S.S.  IS.,  S.S,  0.2 


ENTER  THE  NUMBER  OF  ELEMPTTS  TO  IE  GEN¬ 
ERATED  IN  THE  ALPHA,  BETA  AND  ZETA  CO¬ 
ORDINATE  DIRECTIONS.  RESPECTIVELY  ....  S,  S.  1 


■COIN 


ION 


NUMBER  OF  NODES  TO  IE  QEHDWTSd  • 
NUMBER  OF  ELEMENTS  TO  BE  GENERATED" 


THE  THICKNESS  DIRECTION  OF  THE  MODEL  (IF  ONE  EXISTS.  AS  IN  A 
THICK  SMELL)  RUST  IE  IDENTIFIED  TO  ORIENT  ELEMENTS  PROPERLY. 


OPTIONS  ARE  (DALFHA,  (S)SCTA,  (3)ZETA,  OR  (4 )<UN IMPORT ANT>. 
ENTER  THICKNESS  DIRECTION  CODE  (1.2, 3. 4)  -3 


DATA  OEMDIATIQN  COMPLETE 
1JX0BN  TERMINATED 


USA  MAXIMUM  EXECUTION  FL. 

S.17S  CP  SCCONOO  EXECUTION  TIME. 
'  REWIND, UNFMT 
wnnr  COPVBF.UNFWT, TAPE  1 1 
EDI  ENCOUNTERED  AFTER  COPY  OP  FILE 


S.  RECORD  1 
RETURN, UMFNT. 


COMMAND-  BECIN.IJKGEN.P. 


mmmiiummutnunnmmmt 
tUJIiWt  BEGIN  IJKGEN  lllllllttl 
nnxnotmtixtttunxitimmnmti 

1  JUEN  -  GENERATION  OF  GEOMETRIC  MESH  DATA  FOR  SOLID, 
THICK  SHELL  OR  PLATE  FINITE  ELEMENT  MODELS,  USING  AN  INTEGER  - 
COORDINATE  INDEXING  SCHEME.  OPTIONAL  USER  ROUTINES  ARE  - 
(t)  SURF AC  (1,J,K, ALPHA, BETA, ZETA)  -  DEFINE  MESH  GEOMETRY 

(2)  CRDTRN  (ALPHA, BETA, ZETA, X.Y.Z )  -  COORD.  TRANSFORMATION 
BUILT-IN  OPTIONS  INCLUDE  RECTANGULAR.  CYLINDRICAL  OR  SPHERICAL 
COORDINATES,  AND  UNIFORM  OR  PROPORTIONALLY  GRADED  MESH  SPACING 

(3)  UINPUT  -  USER  PARAMETER  INPUT  ROUTINE  (INITIALIZE  DATA 

BLANK  COMMON) 


Wntmt  USER  SUBROUTINE  'SURF AC'  NOT  CIUEN  XXXXXXXXXX 

BUILT-IN  MESH  DIVISION  OPTIONS  ARE  AS  FOLLOWS  - 

(1)  -  UNIFORM  MESH  IN  EACH  DIRECTION 

(2)  -  GRADED  MESH  (SPECIFY  RATIO  OF  FIRST/LAST  ELEMENT  SIZE 

ENTER  OPTION  (1,2)  . 1 

mmmx  user  subroutine  'CRDtrh'  not  given  «mn«n 

BUILT-IN  COORDINATE  SYSTEM  TRANSFORMATION  OPTIONS  ARE  - 

(1)  RECTANGULAR.  (2)  CYLINDRICAL,  (3)  SPHERICAL 

ENTER  COORDINATE  SYSTEM  OPTION  (1.2.3)  -2 

«  PLEASE  NOTE  THE  FOLLOWING  CONVENTIONS  FOR  CYLINDRICAL  SYSTEM  XX 

'ALPHA'  -  RADIUS  'BETA'  •  ANGLE  'ZETA'  •  AXIAL 
SYSTEM  IS  RIGHT-HANDED.  WITH  ALL  ANGUS  MEASURED  IN  DEGREES. 

ENTER  LIMITING  SURFACE  COORDINATE  VALUES  - 

1.  ALPHA(MIN)  2.  ALPNA(HAX) 

3.  BETA  (MIN)  4.  BETA  (MAX) 

S.  ZETA  (MIN)  S.  ZETA  (MAX) 

G.O,  S.2,  O.O,  00.0,  0,0,  10.0 

ENTER  THE  NUMBER  OF  ELEMENTS  TO  BE  GEN¬ 
ERATED  IN  THE  ALPHA,  BETA  AND  ZETA  CO¬ 
ORDINATE  DIRECTIONS,  RESPECTIVELY  ....  1,4,S 

XXX  BEGIN  GENERATION  PHASE  XXX 

NUMBER  OF  NODES  TO  BE  GENERATED  •  297 
NUMBER  OF  ELEMENTS  TO  BE  GENERATED-  20 

THE  THICKNESS  DIRECTION  OF  THE  MODEL  (IF  ONE  EXISTS,  AS  IN  A 


OPTIONS  ARE  (1 (ALPHA,  (B)BETA,  (3>ZETA,  OR  (4 ><UN IMPORT ANT> 
ENTER  THICKNESS  DIRECTION  CODE  (1.2, 3,4)  -  1 


XXX  DATA  GENERATION  COMPLETE  XXX 
XXX XXX  IJKGEN  TERMINATED  XXXXXX 


020100  MAXIMUM  EXECUTION  FL. 

0.141  CP  SECONOS  EXECUTION  TINE. 
COMMAND-  REWIND. UNF NT. 

COMMAND-  COPVBF.UNFMT.TAPEIO. 

E01  ENCOUNTERED  AFTER  COPY  OF  F1U 


0,  RECORD  1 
RETURN, UNFMT . 


AD-A129  923  MAGNA  (MATERIALLV  AMD  GEOHETRICfN.LV  NONLINEAR  ANALVSIS)  2/1 
'  PART  II  PREPROCES.  .  (U)  DAVTON  UNIV  OH  RESEARCH  INST 

T  S  BRUNER  ET  AL.  DEC  82  UDR-TR-82-U2 
UNCLASSIFIED  AFHAL-TR-82-2898-PT-2  F22615-80-C-2403  .  F/G  9/2  NL 


pa 

PREP  Execution 


COMNAHO-  BEGIH.PREP.P. 


imumnn  begin  mcp  mwnwTii 
nu  NAOHA  PREPROCESSING  UTILITIES  Util 


PLOTTING  TERMINAL  <  V  .  N  )  ?,,V 
TEKTRONIX  OR  HP  (  T  ,  H  )?  ,jT 

TEKTRONIX  TERMINAL  TYPES  - 

0.  RM-I 

1.  401#  /  4#12  /  4#13  ✓  4#52 
3,  4#14  /  4#lS 

3.  4#  14  /  4#lS  (  ENH.  GR.  NOD.  ) 

4,  4114 

ENTER  TERMINAL  TYPE  <  •  -  4  )»,3 
ENTER  CHARACTERS  PER  SECOND  tjlM 
NUMBER  OF  INPUT  DATA  FILES  ?  «|3 
ENTER  FILE  *  IjlS 

ENTER  LABEL  l|CVL-l 

ENTER  FILE  B  lilt 

ENTER  LABEL  jjPLT-1 


HELP'  FOR  LIST  OF  CONNAHM 
it  LIST 

LISTING  OF  DATA  FILES  AVAILABLE 


SIR 

Bll 


CYL-1 

PIT-1 


i«  PLOT 

ENTER  INPUT  FILE  SI  LABEL 


j jCVL-1 


INITIALIZATION  OF  PLOTTING  ROUTINE  ERASES  SCREEN, 
NEABV  (V.N)  ?  ....I  Y 


ENTER  INPUT  FILX  *1  LABEL 
ENTER  OUTPUT  FILX  LABEL 
CREATING  NCU  DATA  FILX 
MSK  OPTIONS  - 

1.  3-D  ELEMENTS  ONLY 
X.  2-0  ELEMENTS  ONLY 
ENTER  OPTION... ill 

MX  NOOC  POSITION  S  IN  PINAL  ELCNCNTS? 
(s.is.2s.as.a7)  uis 

MSK  COMPLETE. 

PERFORM  SUBSEQUENT  SIFT  ?  (Y.N)  ->Y 

SIFT  COHPtXTE 
139  UNUSED  NODES  DELETED 

><  LIST 

USTIN8  OP  DATA  FILES  AUAXLADLE 

SIS  -  CYL-1 

til  -  PLT-1 

112  -  CVL-1SM 

«»  PLOT 

ENTER  INPUT  FIU  Si  LABEL 


IjCVL-l 

ilCVL-lSN 


I, PLT-1 


AT  CMS  OF  PLOT*  ENTER  CHARACTER  TO  CONTINUE 

■Easy  (y.n>  ?  ....»  v 


«»  l«T 

LISTING  Of  MTS  FILES  MANILAS L£ 


SIS 

•12 

•13 


ALT-GOT 

CYL-16N 

PLT-1GN 


t*  TRANSLATE 

ENTER  INPUT  FILE  *1  LABEL  >j ALT -ROT 

ENTER  OUTRUT  FILE  LABEL  jjPLT-RDV 

CREATING  NEW  DATA  FILE 

ENTER  TRANSLATION  F ACTORS ( XFAC. YF AC. 2F AC) - 1  SGG.  G. 


It  DELETE 
ENTER  LABEL 

s< 


,, ALT-ROT 


ENTER  INPUT  FILE  Bi  LABEL 
ENTER  INPUT  FILE  SR  LABEL 
ENTER  OUTPUT  FILE  LABEL 
CREATING  NEU  ORTA  FILE 


jjPLT-RDV 
IlCVL-lSN 
j jSOTH 


TIDY  OPTION  ONE  OR  TUO  ?  (0PT1.  OPTS.  HELP)  ->OPTR 

TIDY  COMPLETE 
11  DUPLICATE  NODES  DELETES 


n  PLOT 

ENTER  INPUT  FILE  SI  LABEL  ijDOTH 

SBRAU 


0.  1G.G 


AT  ENG  OF  PLOT.  ENTER  CHARACTER  TO  CONTINUE 
HEADY  (V.N)  ?  ....<  V 


LISTING  or  DATA  FILES  AVAILABLE 


lt« - DOTH 


><  RENUMBER 

ENTER  INPUT  FILE  11  LAtEL 
ENTER  OUTPUT  FILE  LAtEL 
CREATING  NOI  DATA  FILE 


j jBOTH 
titOTH/R 


NODE  POINT  REORDERING  CONPLETE 


NODAL  BANDWIDTH  (OLD) 
NODAL  BANDWIDTH  (NEW) 


59 

59 


tt  DELE 
ENTER  LABEL 


jjBOTH/R 


It 


ENTER  INPUT  FILE  B1  LABEL 
ENTER  OUTPUT  FILE  LABEL 
CREATING  NEW  DATA  FILE 


1 1  BOTH 

jjCEOM  ♦  DCS 


BOUNDARY  CONDITIONS  NAY  BE  ENTERED  FOR  RANDON 
NODES,  RANGES  OF  NODES,  ALL  NODES  ON  A  PLANE, 
OR  THE  ENTIRE  NOBEL. 


rOR  EACH  GROUP  OT  NODES.  CONSTRAINTS  NAY  BE  SET 
IN  ANY  COMBINATION  OT  THE  X.  Y,  &  2  DIRECTIONS. 
(NOTE  THAT  2-D  NOBELS  RUST  BE  FIXED  IN  2) 


CONSTRAINT  DIRECTION  CODE  CONSTRUCTED  AS  FOLLOWS  l 


1  -  '1' 
_ .  'O' 

ENTRY  2  - 

3  - 


FOR  CONSTRAINED  IN  X  DIRECTION 
FOR  UNCONSTRAINED  IN  X  DIRECTION 
SANE  FOR  V 
FOR  2 


-> 


'1.0.1'  FOR  CONSTRAINTS  IN  X  AND  2 
DIRECTIONS  BUT  NOT  IN  V. 


1  - 

S  -  RANGE  OF  NOSES 
3  -  SPECIFIED  PLANE 

S  "  EXIT 


SELECTION  OPTION  (1.E....S) 
ENTER  CONSTRAINT  DIRECTION  CODE  (1 


->3 

->•  0  1 


IS  DEFINED  BY  AX  ♦  BY  ♦  CZ  •  D. 
ENTER  COEFFICIENTS  (  A  ,  B  .  C  .  D)  - 


•>0  0  10 


43  NODES  FOUND 


1  •  RANDOM  NODES 
3  -  RANGE  OF  NOOES 

3  -  SPECIFIED  PLANE 

4  -  ALL 

5  -  EXIT 


ENTER  NODS  SELECTION  OPTION  (1.E....S) 


->3 


ENTER  CONSTRAINT  DIRECTION  CODE  (THREE  VALUES)  ->1  1  0 

— - >1  0  0  0 

BE  NODES  FOUND 


PUNE  IS  OEFINED  BY  AX  ♦  BY  ♦  C2  «  D. 
ENTER  COEFFICIENTS  (  A  .  B  ,  C  ,  D) 


1  -  RAHOON  NODES 
3  -  RANGE  OT  NOOES 

3  -  SPECIFIED  PLANE 

4  -  AU.  NODES 


LISTING  OF  MTU  FILES  AVAILABLE 


*10  -  ALL  BUT  PROPS 

BH - FINAL 


BP 

TAPE14 


STOP 

tS41tt  WAX  I  NUN  EXECUTION  FL. 

1B.S70  CP  SECONDS  EXECUTION  TINE. 

FILE  QUOTA  EXCEEDED 
COWHAND-  FILES 
—LOCAL  FILES— 

TAPE1B  TAPE IB  SINPUT  SOUTPUT 

TAPE11  TAPE13  TAPC1S  TAPE13 

TAPE1S  TAPE17 

FILE  QUOTA  EXCEEDED 

COWHAND-  RETURN, TAPElt,TAPE12,TAP£13,TAPE14,TAPE15,TAPE16,TAPE17 
COWHAND-  FILES 
—LOCAL  FILES— 

TAPE1S  SINPUT  SOUTPUT  BP  TAPElt 

COWHAND*  RETURN, TAPE1S 
COWHAND*  REWIND. TAPE 11 
COWHAND*  COPVBF.TAPEll.UNFWT 
EOI  ENCOUNTERED  AFTER  COPY  OF  FILE 
S.  RECORD  1 
CONNAND-  RETURN, TAPE1 I 


MUM  SOLUTION  OPTIONS  MM  WMHTBI 


ENTER  ANALYSIS  TYPE  (UNBAR.  NONLINEAR)  -NONLINEAR 

ENTtt  ANALYSIS  SUBTYPE  (STATIC. DYNAMIC) - STATIC 

ENTER  TINS  STEP  OPTION  (CONST., VARIABLE)  —CONST 

SNTSK  THE  INITIAL  SOLUTION  TINE  STEP  - IS. 

ENVER  THE  NUMBER  OP  SOLUTION  TINE  STEPS  - IS 


EQUILIBRIUM  ITERATION  OPTIONS  ARE  AS  FOLLOWS 
•  •  NO  ITERATION 

1  *  MODIFIES  NCUTON  (CONST.  STIFFNESS) 

2  •  FULL  NCUTON- RAPNSON  ITERATION 

3  ■  CONSINED  FULLs MODIFIED  NCUTON 

ENTER  ITERATIVE  SOLUTION  OPTION  (S.  1.2.3)  —3 

POSTPROCESSOR  FILE  TO  BE  WRITTEN  (Y/N)  — Y 

ENTER  TME  FREQUENCY  (IN  INCREMENTS)  AT  UNION 
RESULTS  ASS  TO  SAUED  ON  POSTPROC.  FILE  - S 

ARC  RESTART  FILES  TO  IE  READ  (Y/N) . N 

ARE  RESTART  FIUS  TO  BE  WRITTEN  <Y/N> - — Y 


ENTER  NEW  RESTART  FILE  LABEL  (4  CHARS. )  — — FPLT 


XXXXXmmn 
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APPENDIX  A 


PREPROCESSOR  FILE  DESCRIPTIONS 


The  preprocessor  described  in  this  report  deals  with  two 
internal  forms  and  two  external  forms  of  finite  element  model 
data.  These  forms  are  discussed  briefly  in  Section  6  -  Model 
Data  Files.  Table  A.l  describes  the  unformatted  2-D  internal 
data  file  that  is  utilized  for  shell-like  surface  models 
defined  via  the  programs  CORGEN,  AGRID  and  SPATCH  described  in 
Section  2.  The  EXPAND  utility  converts  all  shell  surface  model 
files  to  the  standard  unformatted  3-D  internal  data  file  illus¬ 
trated  in  Table  A. 2.  This  file  type  is  accessed  by  the  prepro¬ 
cessor  and  the  data  translation  programs.  Once  the  user  has 
defined  the  model,  there  exist  two  external  file  types  for 
transferring  data  from  the  preprocessor  to  other  programs  or 
computer  systems.  Table  A. 3  describes  the  neutral  data  file 
format,  which  can  be  used  to  transfer  data  between  the  prepro¬ 
cessor  and  other  programs,  to  transmit  data  between  computer 
systems,  or  to  archive  modeling  data  in  a  form  which  can  be 
listed  on  a  terminal  or  line  printer.  The  MAGNA  input  file, 
which  is  input  to  TRNSFR  and  output  from  REFMT,  is  described 
in  detail  in  Reference  [4] . 


TABLE  A.  1 


SHELL  SURFACE  GEOMETRY  FILE 


File  description:  This  is  an  unformatted  (binary) ,  sequential  file 

of  surface  (2-D)  model  data  output  from  the  data 
entry  programs  AGRID,  CORGEN  and  SPATCH.  The 
numbering  of  node  points  and  elements  is  sequen¬ 
tial  with  numbering  implicit. 


Section 

#  Records 

Variables 

I 

1 

NUMNP 

NUMEL 

II 

NUMNP 

COORD ( 4 ) 

III 

NUMEL 

NC0N(9 ) 

Description 

Number  of  node  points 
Number  of  elements 

Cartesian  coordinates  and  thick¬ 
ness  for  node  point 

Nodal  connectivity  for  element 


TABLE  A. 2 


THREE-DIMENSIONAL  MODEL  DATA  FILE 


File  description?  This  is  an  unformatted  (binary) »  seauentlal  file 
of  three-diaensiohal  aodal  data  output  fro*  the 
data  entry  rrosraas  CREATE*  EXPAND  and  IJKGEN 
and  tha  intarfaca  routinas  (a. a.  NEUTRAL*  IMPRINT 
and  TRN8FR  . 

Section  ♦  Records  Variables  Description 


I 

1 

NUMNP 

NUMEL 

NUMBC 

II 

NUMNP 

NODNUM 
COORD (3) 

III 

NUMEL 

IELNUH 

I  TYPE 

MATRL 

IAXSET 

INTORD 

IDUH(3) 

NC0N(27) 

IV 

NUMBC 

I8E0 

IEND 

I  NCR 
IDIR(3) 

V 

1 

NUMNL 

VI 

NUMNL 

NBEO 

NENB 

I  NCR 

I  CASK 

F  (3) 

VII 

1 

NUMLE 

VIII 

NUMLE 

NVCS 

NCND 

I  NCR 
ISURF 
ICASC 
FORCE 

IX 

1 

NUMLC 

X 

NUHLC 

NODE 

XMULT 

YMULT 

ZHULT 

Nuabar  of  node  Points 

Number  of  elements 

Nuabar  of  boundary  conditions 

Node  point  number  (arbitrary) 
Cartesian  coordinates  (X*Y*Z) 

Eleaent  nuabar  (arbitrary) 
Eleaant  tap*  code 
Material  type  code 
Material  axis  code 
Inteeration  order 
Unused 

Nodal  connectivity  for  elements 
(variable  1-27  nodes  solid) 

First  node 
Last  node 
Node  increaant 

List  of  constrained  directions 

Nuabar  of  nodal  loads 

First  noda  . 

Last  node 
Nod*  increaent 
Case  number 

Nodal  force  vector  (X*Y*Z) 

Number  of  eleaent  loads 

First  eleaent 
Last  eleaent 
Eleaent  increaent 
Surface  or  edse  nuaber 
Case  nuaber 

Force  per  unit  area/lenath 
Nuaber  of  linear  constraints 
Nod*  nuaber 

Multiplier  for  X-displece**nt 
Multiplier  for  Y-displaceaent 
Multiplier  for  Z-displaceaent 


TABLE  A. 3 


NEUTRAL  FILE  DATA  FORMAT 


File  description!  This  is  a  forastted  (character ) »  seouential  fils 
of  three-dieensional  eodel  data  output  fros  the 
NEUTRAL  prodraa*  Rseords  ars  all  80  co loans* 


Saction  0  Rseords  Variable  Forsat  Description 


I 

1 

NUHNP 

IS 

NUMEL 

IS 

NUHBC 

IS 

II 

NUHNP 

MQBMUH 

IS 

COORD <  3 ) 

E13*i 

III 

NUMEL 

IELNUH 

IS 

ITYPE 

IS 

MATRL 

15 

IAX8ET 

IS 

INTORD 

IS 

IB(JM(3) 

IS 

NC0N(27) 

15 

IV 

NUHBC 

IBES 

IS 

I  END 

IS 

I  NCR 

IS 

IDZR<3> 

IS 

V 

1 

NUMNL 

IS 

VI 

NUMNL 

NBCO 

IS 

NEMO 

IS 

INCH 

IS 

I  CASE 

15 

F<3> 

3tlS.S 

VII 

1 

NUMLE 

IS 

VIII 

NUMLE 

NBEB 

IS 

NEMO 

IS 

I  NCR 

IS 

I  SURF 

IS 

ICASE 

IS 

FORCE 

CIS. 8 

IX 

1 

NUHLC 

IS 

X 

NUHLC 

NOSE 

IS 

XNULT 

C1S.S 

THULT 

£15. 8 

ZHULT 

ElS.i 

Nuabar  of  nods  saints 

Nuabar  of  alaaants 

Nuabar  of  boundary  conditions 

Noda  soint  nuabar  (arbitrary) 
Cartaaian  eoordinatas  (X»Y*Z) 

Elaaant  nuabar  (arbitrary) 
Elaaant  tuna  coda 
Material  type  coda 
Notarial  axis  coda 
Intadration  order 
Unusad 

Connectivity  for  alaaants 
(variable  1-27  nodes) 

First  node 
Last  node 
Noda  increaant 

List  of  constrained  directions 

Nuabar  of  nodal  loads 

First  node 
Last  Noda 
Node  increaent 
Casa  nuabar 

Nodal  force  vector  (XrT.Z) 

Nuabar  of  aleaent  loads 

First  elaaant 
Last  elaaant 
Eleeent  increaent 
Surface  or  adds  nuabar 
Case  nuabar 

Force  ear  unit  area/leneth 
Nuabar  of  linear  constraints 
Node  nuabar 

Multiplier  for  X-disPlsceeent 
Multiplier  for  Y-displaceaent 
Multiplier  for  Z-disPlaceeent 


APPENDIX  B 


POSSIBLE  DATA  PATHS 


Tables  B.l  and  B.2  contain  a  brief  summary  of  the  preprocessing 
programs  and  data  file  types  referred  to  in  this  report.  These,  in 
conjunction  with  the  flow  chart  in  Figure  1.1,  are  intended  to  aid  in 
determining  the  best  path  to  follow  through  the  system  in  preparing 
a  finite  element  analysis  model. 

For  a  particular  application,  the  type  of  data  available  for  use 
should  be  compared  with  the  data  entry  descriptions  given  in  Table 
B.l.  The  type  of  data  generated  by  any. of  the  data  entry  or  translatio 
utilities  can  then  be  determined  from  Table  B.2.  The  data  type 
generated  can  next  be  compared  with  the  input  file  types  in  Table  B.2 
to  determine  which  utility  (usually  EXPAND  or  PREP)  will  be  executed  #■'* 
next.  Once  the  model  data  is  stored  in  the  internal,  three-dimensionax 
file  format  used  by  PREP,  the  data  path  will  generally  take  the  form 
PREP  — >  REFMT  — >  MAGNA  for  analysis,  or  PREP  — >  NEUTRAL  for  the 
translation  to  the  format  required  by  other  programs. 


PREPROCESSOR  MODULE  FUNCTIONS 


INITIAL  DATA  GENERATION 


AGRID 


CORGEN 


CREATE 


EXPAND 


IJKGEN 


SPATCH 


Definition  of  surface  geometry  via  an 
arbitrarily  arranged  grid  of  points. 

Lofting  and/or  direct  data  input  of  arbitrary 
surface  geometry. 

Direct  input  and  editing  of  coarse  mesh 
geometry. 

Expansion  of  surface  geometry  description 
ii.to  three-dimensional  model  form. 

Mesh  generation  on  analytically-defined 
surfaces . 

Model  generation  from  surface  'patch* 
description. 


INTERFACING  MODULES 
IMPRINT  : 

NEUTRAL  s 

TRNSFR  : 


PREPROCESSING 

PREP 

REFMT 


Conversion  of  IMPRESS  preprocessor  file  to 
MAGNA  preprocessor  format. 

Translation  of  model  data  between  internal 
format  and  external  (text)  format. 
Conversion  of  MAGNA  input  to  preprocessor 
format. 


Model  refinement,  merging,  plotting, 
properties  definition,  etc. 

Generate  MAGNA  input  data  deck. 


TABLE  B .  2 

PREPROCESSOR  FILE  TYPES 


1 

FILE  TYPE 

1  TYPE  FILE 

1  OUTPUT  FROM  I  INPUT  TO 

1 

l 

1  Internal  3-D 

1  Unformatted 

1  CREATE 

1  PREP 

1 

1  Data  File 

!  Sequential 

|  TRNSFR 

1  REFMT 

1 

IJKGEN 

1  NEUTRAL 

1 

1 

1 

• 

1 

PREP 

EXPAND 

NEUTRAL 

1  IMPRINT 

1 

1 

1 

1 

2 

1  Internal  2-D 

Unformatted 

1  CORGEN 

1  EXPAND 

1 

1 

1  Data  File 

1  Sequential 

1  SPATCH 

1  AGRID 

1 

J 

1 

3 

*  Material 

Formatted 

1  — 

|  PREP 

1 

1 

1  Properties 

1  Library  File 

!  Sequential 

1  REFMT 

1 

! 

1 

4 

I  Neutral  Data  File 

Formatted 

Sequential 

1  NEUTRAL 

|  NEUTRAL 

1 

1 

S 

1  MAGNA  Input  Deck 

1 

Formatted 

|  REFMT 

1  MAGNA 

1 

Sequential 

|  TRNSFR 

_  APPENDIX  C 


ACCESS  PROCEDURES  FOR  PREPROCESSOR  MODULES 


Procedures  for  accessing  the  individual  preprocessor  modules  are 
summarized  in  Tables  C.l  and  C.2  for  quick  reference.  On  CDC  machines 
all  of  the  programs  may  be  accessed  through  a  single  CCL  (CYBER  Contro 
Language)  command  procedure,  PREPROCESSORPROC?  execution  is  initiated 
by  entering  a  single  BEGIN  command  as  described  in  Table  C.l.  On  the 
VAX  11/780,  the  DCL  (Digital  Command  Language)  procedure  [MAGNA. RAB] 
CONTROL  may  be  used  to  control  execution,  or  a  RUN  command  may  be 
issued  directly. 


s 

( 

a 

■ 

• 

■ 

U 

• 

II 

1 

a 

1 

1 

f 

1 

1 

1 

1 

1 

1 

1 

1 

( 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

i 

i 

i 

i 

i 

i 

i 

i 

1 

N 

1 

1 

1 

1 

» 

1 

1 

1 

ft 

ft 

1 

1 

i 

i 

• 

B 

1 

1 

1 

1 

0 

1 

i 

s 

■ 

1 

1 

1 

m 

1 

0 

1 

i 

o 

R 

1 

1 

1 

1 

w 

1 

i 

1 

1 

1 

1 

k* 

R 

1 

1 

1 

1 

V 

1 

i 

o 

R 

1 

1 

1 

1 

ft 

1 

i 

V 

a 

1 

1 

1 

1 

*» 

1 

i 

« ■ 

R 

1 

1 

1 

1 

£ 

1 

i 

• 

■ 

1 

1 

1 

l 

ft 

1 

i 

• 

R 

R 

1 

I 

t 

1 

1 

1 

1 

1 

3 

e 

1 

1 

i 

i 

R 

I 

1 

1 

l 

9 

1 

i 

• 

■ 

1 

1 

1 

1 

•  w 

1 

i 

*  - 

X 

a 

1 

pn 

1 

m 

1 

1 

1 

i 

o 

a 

• 

1 

1 

• 

1 

ft 

1 

i 

• 

M 

ft 

1 

ft  ft 

1 

ft 

• 

ft 

n  • 

•  • 

1 

0 

1 

0  ♦ 

1 

• 

i 

•  m 

•» 

a 

ft 

1 

ft  ft 

1 

ft 

ft 

ft  ft 

ft 

ft  ft 

1 

ft 

1 

ft 

1 

ft 

ft  ft 

i 

ft  » 

«» 

a 

ft 

1 

ft  ft 

1 

ft 

ft 

ft  ft 

ft  ft 

ft  ft 

1 

ft  ft 

1 

ft 

1 

ft  ft 

ft  ft 

i 

ft  ft 

■. 

• 

» 

c 

1 

£  C 

1 

C 

ft 

ft  £ 

ft  ft 

ft  ft 

1 

ft  ft 

1 

•  M 

1 

ft  ft 

ft  c 

i 

•  • 

a 

ft 

1 

ft  ft 

1 

ft 

c 

£  ft 

ft  C 

£  C 

1 

ft  c 

1 

£  N  • 

1 

ft  £ 

C  ft 

i 

-  £ 

m 

a 

mm 

1 

m 4  m m 

1 

mm 

ft 

ft  0 

S  ft 

•  ft 

1 

C  ft 

1 

ft  0  ft 

1 

£  ft 

ft  mm 

i 

ft  ft 

• 

mm 

a 

•Vi 

1 

mt  mm 

1 

mm 

mm 

k  0 

1 

ft  0 

1 

0  0  ft 

1 

ft  0 

mm  mm 

» 

—  mrn 

. 

a  ft- 

W 

1 

0  0 

1 

0 

m» 

mt  0 

0*0 

1 

0  •  0 

1 

0  0  •  •  ft 

1 

mm  mm 

mm  0 

i 

mm  k 

j 

X 

■  ft 

* 

1 

k  •  k 

1 

k  • 

0 

V  .  * 

0*0 

0  0  0 

1 

0  0  0 

1 

0  0  •  0  £ 

1 

mm  k. 

0  •  k 

i 

0B0 

4 

u 

■  ft 

•  U 

1 

X0  o 

1 

O  0 

m 

*  0  0 

0  0k 

k  »  m 

1 

0  •  * 

1 

•  •  •  ft  k  ft 

1 

0  •  » 

k  0  U 

i 

k  0  k 

w 

mm 

a  u 

0  0 

1 

X  »  k 

1 

0  k 

0 

O  *  0 

»  *  0 

O  0  0 

1 

•  ml  0 

1 

1 

k  0  0 

©  k  0 

i 

0  k  tk 

X 

a  o 

k  0 

1 

0X0 

1 

X  0 

X 

0  0  0 

m  z  x 

0X0 

0  0  0 

1 

1 

0  »  0 

»  x  0 

i 

0  X  X 

• 

3 

a  0 

00 

1 

00  0 

1 

0  0 

0 

0X0 

9  0  0 

0  0  0 

1 

X  X  0 

1 

0  >UI«€Uk 

1 

X  0  0 

0  U  X 

i 

Qb.li. 

-  “ 

• 

0  9 

000 

1 

X  0 

z 

0  0  0 

a  o  x 

0X0 

1 

0  0  a 

1 

0  0  0  k  0  0  » 

1 

0  x  a 

00  9 

i 

0  a  x 

X 

» 

a 

0  0 

1 

0  0  0 

1 

9  0 

9 

0  00 

9X9 

0  0  0 

1 

X  9  0 

1 

00  0  0  •  0  0 

1 

9  0  0 

0  0  a 

i 

X  x  9 

ui 

• 

* 

®  * 

1 

1 

-  0 

a« 

m-  X  m 

»  -a  • 

0  X  k 

1 

9  0k 

1 

0  0  0  -00  k 

1 

9  0k 

0  0k 

i 

0  X  - 

0 

mm 

a 

0  0 

1 

X  U  CD 

1 

X  U 

o 

•  0  © 

z  0  a 

•  0  a 

1 

•  x  a 

1 

—  0  »  k  k  »  Q 

1 

*  x  a 

k  a  o 

i 

k  0  a 

OB 

mm 

a 

•  O 

1 

U  •  O 

1 

u  » 

o 

X  -  o 

U  k  o 

x  •  o 

1 

X  mO 

1 

x  k  x  a  a  w  o 

1 

X  k  o 

X  •  o 

i 

X  k  o 

. 

► 

0 

a 

X  -J 

1 

0X  mi 

1 

0  X 

0 

u  X  0 

0X0 

U  X  mJ 

1 

0  X  mi 

1 

u  x  x  0  x  a  0 

1 

«J  Z  -J 

U  X  mJ 

i 

«J  X  -J 

M 

a 

0  0 

1 

0  0  0 

1 

0  mm 

0 

0  0  0 

0  0  0 

0  0  0 

1 

0  0  0 

1 

0  0  9  9  0  0  0 

l 

0  k.  0 

0  «  0 

i 

«H0 

• 

a 

0  N» 

1 

0  O  0 

1 

0  CD 

0 

0  a  0 

0  a  0 

0  e  0 

1 

0  a  0 

1 

0000X00 

) 

0  a  0 

0  o  0 

i 

0  a  0 

k 

■ 

0  0 

1 

0  0  0 

1 

0  0 

0 

0  0  0 

0  0® 

0  0  0 

1 

0  0  0 

1 

1 

0  0*0 

0  0  0 

i 

t-  ua  « 

0 

3 

a 

0  U 

1 

u  a  u 

1 

0  0 

u 

«  a  u 

U  ft  U 

0  a  u 

1 

0  a  u 

1 

0  a  x  x  x  o  u 

1 

0  a  u 

0  a  u 

i 

«  m  u 

Z 

V 

a 

1 

1 

1 

1 

1 

i 

1 

M 

ft 

a 

1 

1 

1 

1 

1 

i 

0 

O 

a  «►  • 

—  w 

1 

■»  •  • 

1 

—  mm 

• 

•  • 

mm  ^^m  mm 

aa  •  • 

1 

w  «■  a 

1 

1 

mm  mm  mm 

i 

-  *■  — 

0 

0 

a  • 

1 

1 

1 

1 

1 

0 

i 

.» 

0 

k 

a  • 

1 

1 

1 

1 

a 

i 

u 

ft 

a  c 

1 

1 

1 

1 

1 

o 

i 

Q 

a 

1 

1 

m  0 

1 

1 

1 

0 

i 

*  1 

0 

0 

a  «• 

1 

1 

U  X 

1 

1 

I 

u 

i 

0 

u 

a  9 

ft 

1 

ft 

1 

ft 

• 

0  X 

ft 

1 

ft 

1 

• 

1 

ft 

a  x 

i 

ft 

-■ 

0 

u 

a  o 

c 

1 

c 

1 

c 

£ 

0  0  0 

£ 

1 

c 

1 

£ 

» 

c 

9  »  x 

t 

£ 

0 

a  k  k 

0 

1 

0 

1 

o 

O 

9X0 

e 

1 

ft 

t 

O 

1 

0 

a  x  o 

i 

o 

■ 

0 

ft 

a  ft  X 

e 

1 

c 

1 

c 

C 

a  9  x 

£ 

1 

£ 

1 

£ 

1 

c 

3UI- 

i 

£ 

.* 

a  ft  9 

1 

1 

9ftU 

1 

1 

1 

0  0 

i 

1  • 

ft 

a  90 

1 

1 

mm 

1 

1 

1 

O  0 

i 

Mb 

mm 

0 

c  . 

a 

1 

1 

1 

1 

1 

9  9 

i 

f 

Z 

•k 

a 

1 

1 

1 

1 

1 

mm 

i 

4 

U 

a 

ft  • 

§  mm  mm 

—  ■— 

1 

•  •  • 

i 

• 

•  »  • 

~  mm  ~ 

M  M 

1 

«  »  mm 

1 

1 

«—  mm  ~ 

»  —  • 

i 

•  _  « 

< 

«■ 

a 

1 

i 

1 

» 

1 

i 

w  « 

0 

z 

ft  ■ 

a 

» 

l 

1 

1 

1 

1 

i 

,  • 

-I 

0 

a  *■— 

1 

1 

1 

1 

1 

i 

-  ■ 

A 

XX 

»  -ft  • 

(J 

1 

o 

l 

0 

1 

0 

1 

©  M 

1 

u 

i 

•  , 

« 

o 

ft  - 

a  ft  — 

0 

1 

•* 

i 

0 

0 

0 

0 

1 

0  0 

1 

—  M 

1 

0 

X 

i 

0 

'k  • 

0 

0 

3g 

a  4  ft 

0 

1 

0 

i 

X 

0 

X 

0 

1 

a  X 

1 

0  0 

1 

0 

0 

i 

X 

,  a  • 

i  ^ 

0 

1 

0 

l 

0 

0 

0 

0 

1 

0  0 

1 

0  0 

1 

0 

X 

i 

0 

»  ■* 

«-  0 

aft  — 

9 

1 

0 

i 

X 

0 

X 

0 

1 

X  X 

1 

00 

1 

a 

9 

i 

X 

1 

0 

X0 

R  O  Ik 

0 

1 

0 

i 

9 

0 

9 

0 

1 

0  9 

1 

0  0 

1 

0 

a 

9 

* 

0 

k*  0 

a 

1 

l 

1 

1 

1 

i 

0 

0 

■ 

1 

i 

1 

1 

1 

i 

U 

©x 

■  ■»  ■» 

1 

•  •  • 

I 

“  mm 

• 

•  •  • 

»  ~  » 

m*  mm  wm 

1 

«a»  <a*i  w 

1 

1 

k  »  a> 

•  •  • 

i 

mm  w  m 

u 

•  0 

a 

I 

1 

1 

1 

1 

i 

X 

ft  k 

1 

mm  | 

i 

mm 

1 

1 

1 

1 

i 

ft  0 

a 

1 

ft 

I 

ft 

1 

1 

1 

i 

.  ' 

•  0 

■  w 

X 

1 

£0© 

i 

£  O 

o 

© 

« 

0 

1 

©  M 

1 

© 

i 

0 

a 

U  X 

a  ft  • 

O 

1 

O  mm  « 

i 

O  0 

0 

0 

1 

0  0 

1 

0  N 

1 

0 

a 

i 

0 

- 

a 

U  X 

a  4  0 

0 

1 

•*•0  0 

I 

-k  X 

0 

ft 

0 

1 

x  a 

1 

00 

1 

X 

0 

i 

a 

u 

•  u 

R  C  -* 

0 

» 

«*  0  0 

I 

•»  0 

0 

£ 

0 

1 

0  0 

1 

X0 

1 

0 

0 

i 

0 

- 

o 

R  •*  0 

X 

1 

4  0  0 

1 

4  X 

0 

0 

0 

1 

X  X 

1 

00 

1 

X 

0 

i 

X 

ft  0 
00 

a 

a 

» 

1 

1 

00  0 

i 

I 

O  9 

e 

0 

1 

1 

90 

1 

1 

00 

1 

1 

9 

0 

i 

i 

0 

b  ' 

1 

l 

ft  0 

V 

1 

i 

1 

1 

i 

•  M 

■ 

ft 

1 

ft 

i 

ft 

ft 

• 

ft 

1 

ft 

1 

ft 

1 

ft 

ft 

i 

ft 

«  ' 

ft  - 

a 

> 

1 

> 

l 

> 

> 

> 

1 

> 

1 

> 

1 

> 

> 

» 

> 

► 

u 

a  a 

mt 

1 

mt 

l 

0 

0 

1 

0 

1 

mt 

1 

•0 

ft  o 

a  e 

*» 

1 

m* 

i 

m* 

«• 

m* 

1 

1 

m* 

1 

m» 

«• 

i 

m» 

•  0 

a  —  ft 

u 

1 

O 

l 

U 

(1 

O 

1 

«J 

1 

U 

1 

U 

U 

i 

U 

<=*  0 

■  -»a 

ft 

» 

ft 

i 

ft 

ft 

ft 

1 

ft 

1 

ft 

1 

ft 

ft 

i 

ft 

90 

a  •  o 

k 

1 

i 

k 

k 

k 

1 

k 

k 

k 

i 

k 

©o 

R  «-x 

ft 

1 

ft 

i 

ft 

0  U 

ft 

ft 

1 

ft 

1 

ft 

1 

ft 

ft 

i 

ft 

00 

a  ft 

k* 

1 

m* 

I 

0 

0 

m» 

m» 

<•» 

i 

k* 

ft  0 

a  4 

e 

1 

£ 

i 

c 

C  ft 

e 

£ 

1 

£ 

1 

£ 

1 

c 

£ 

i 

£ 

0 

r  e 

0 

1 

mm 

0  m 

0 

0 

1 

0 

1 

mm 

1 

M 

mm 

i 

0 

►*V 

fcO 

a 

1 

I 

1 

1 

1 

i 

oo 

R 

1 

1 

1 

1 

i 

•  0 

d  “  — 

•  • 

1 

*  — 

•  —  * 

•a  •  • 

•  •  • 

1 

mm  mm  mm 

1 

1 

k  •  • 

mm  ^m 

i 

•  •  • 

ftfc 

a 

1 

1 

1 

1 

i 

ft  0 

a 

1 

1 

1 

1 

i 

•> 

u  0 

a 

1 

1 

1 

1 

i 

04 

a  ft 

1 

1 

1 

1 

i 

“fl 

k  » 

a  ft 

1 

0 

mJ 

1 

1 

i 

44 

R  4  ft 

1 

X 

i 

0 

X 

X 

0 

1 

1 

X 

i 

X 

ft  « 

a  ft  ft 

0 

1 

0 

l 

0 

0 

X 

t 

1 

0 

u 

i 

0 

k  X 

a  o  ft 

0 

1 

© 

0 

o 

« 

l 

0 

1 

X 

0 

i 

a 

1U 

R  4  X 

0 

1 

0 

I 

0 

0 

X 

0 

9 

t 

0 

1 

0 

0 

i 

x 

0 

R  0 

0 

1 

o 

0 

-J 

X 

0 

l 

X 

1 

0 

0 

i 

X 

00 

0  0 

tl 

R 

0 

1 

I 

u 

l 

u 

M 

0 

1 

1 

X 

1 

1 

0 

1 

1 

X 

a 

i 

i 

0 

<0 

*1 

—  — 

1 

t 

-  —  - 

I 

— 

- 

— 

— 

_  «•  _ 

1 

1 

•  — .  — 

1 

1 

1 

1 

—  —  — 

— 

i 

i 

—  _  _ 

A. 10 


4 


APPENDIX  D 


USER  SUBROUTINES  FOR  IJKGEN 


The  user  subroutines  for  IJKGEN  (UINPUT,  SURFAC,  and  CRDTRN )  have 
the  following  general  forms: 

SUBROUTINE  UINPUT 
COMMON  <variable  list> 


<  user  -  written  FORTRAN  code  to  read  or  initialize 
variables  in  blank  COMMON  > 


RETURN  *•] 

END 


SUBROUTINE  SURFAC  (  I,  J,  K,  ALPHA,  BETA,  ZETA  ) 


COMMON  <variable 

list> 

COMMON  /  LIMITS  / 

IMAX  , 

JMAX  , 

KM  AX 

+ 

ALFALO  , 

ALFAHI  , 

BETALO 

+ 

BETAHI  , 

ZETALO  , 

ZETA HI 

<  user  -  written  FORTRAN  code  to  define  the  curvilinear 


coordinates  ALPHA,  BETA,  and  ZETA  as  functions  of  I, 
J,  and  K.  The  generated  coordinates  should  form  a 
right-handed  coordinate  system  but  are  otherwise 
arbitrary  > 


• 

RETURN 

END 

SUBROUTINE  CRDTRN 

(  ALPHA, 

BETA,  ZETA, 

X,  Y, 

COMMON  <variable 

list> 

COMMON  /  LIMITS  / 

IMAX  , 

JMAX  , 

KM  AX 

+ 

ALFALO  , 

ALFAHI  , 

BETALO 

+ 

BETAHI  ., 

ZETALO 

ZETAHI 

<  user  -  written  FORTRAN  code  to  define  the  transformation 
from  curvilinear  coordinates  ALPHA,  BETA,  and  ZETA  to 
right-handed  Cartesian  coordinates _X,  Y,  Z  > 


RETURN 

END 


Normally,  UINPUT  will  not  be  used  unless  at  least  one  of  the  other 
routines  is  also  supplied  to  IJKGEN.  When  UINPUT  is  used  to  set 
the  values  of  parameters  in  blank  COMMON,  these  values  will  not 
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be  modified  elsewhere  within  the  program  and  may  be  used  and/or 
modified  as  needed  within  SURFAC  and  CRDTRN. 

The  named  COMMON  block  /LIMITS/  contains  data  concerning  the 
limits  of  the  mesh  being  generated,  which  is  often  useful  within 
the  user  subroutines  SURFAC  and  CRDTRN.  The  three  integer  values 
IMAX,  JMAX,  KMAX  contain  the  maximum  values  of  the  generator 
indices  (their  lower  limit  is  always  one),  and  the  remaining 
parameters  describe  the  lower  and  upper  limits  of  each  of  the 
curvilinear  coordinates,  which  are  requested  as  keyboard  input 
during  execution.  While  the  data  in  COMMON  /LIMITS/  is  always 
accessible  to  the  user  routines,  the  data  values  should  not  be 
modified  to  ensure  proper  operation  of  the  program. 


APPENDIX  E 


USER  SUBROUTINES  FOR  SPATCH 

The  three  available  user-written  subroutines  for  SPATCH  are: 

-  UPATCH  -  define  surface  patch  coefficients  and  generation  dati 

-  UOPEN  -  open  and/or  position  input  file 

-  UCLOSE  -  close  input  file 

In  most  instances,  only  the  first  routine,  UPATCH,  will  be  needed  to 
perform  the  data  conversion. 


The  general  form  of  UPATCH  is  as  follows: 


C 

C 

C 

C 

C 

C 

C 

C 

C 

C 


SUBROUTINE  UPATCH  (  PATCH,  THICK,  NU,  NV,  ICOUNT,  I END  ) 
DIMENSION  PATCH  (  4,  4,  3  ) 


PATCH 


Output  array  of  patch  coefficients,  having  the  same  form 
as  in  equation  (2.3). 

PATCH  (i,j,l)  *  coefficient  array  (4x4)  for  X-coordinate 
PATCH  (i, j,2)  =  coefficient  array  (4x4)  for  Y-coordinate 
PATCH  ( i , j , 3 )  =  coefficient  array  (4x4)  for  Z-coordinate 


THICK  ■  Thickness  value  for  the  patch,  for  use  in  EXPAND 
NU  =  Number  of  quadratic  elements  to  be  generated  from  this 
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surface  patch,  along  the  u-coordinate  direction 


NV 


ICOUNT  = 


I  END 


Number  of  quadratic  elements  to  be  generated  from  this 
surface  patch,  along  the  v-coordinate  direction 


Sequence  number  for  this  patch.  ICOUNT  is  provided  in 
case  the  user  routine  must  keep  track  of  the  number  of 
patches  read.  When  UPATCH  is  entered  for  the  n-th  time 
the  value  of  ICOUNT  =  n,  unless  modified  by  the  user. 


End  of  data  flag,  to  be  set  in  UPATCH.  If  no  more  data 
exists,  I END  should  be  set  to  1  to  inform  SPATCH  that 
all  of  the  data  has  been  read.  If  IEND=0  on  exit  from 
UPATCH,  SPATCH  will  generate  surface  data  from  the 
parameters  PATCH,  THICK,. .NU  and  NV. 


<  user-written  code  > 


RETURN 

END 


The  UPATCH  routine  may  read  data  from  any  type  of  file  (random, 
sequential,  formatted,  or  unformatted) .  For  this  reason,  a  second 
user-written  routine,  UOPEN,  is  provided  to  permit  initialization 
of  the  file  to  be  read.  The  form  of  UOPEN  is  simply 


SUBROUTINE  UOPEN 


<  code  to  initialize  file  > 


RETURN 

END 


On  CDC  machines,  the  local  file  TAPE50  is  reserved  for  use  by  UPATCH, 
and  the  UOPEN  routine  need  only  be  provided  if  the  surface  patch  data 
is  stored  on  a  random  access  file;  in  this  instance,  UOPEN  may  be  used 
to  initialize  the  random  file,  by  calling.  OPENMS,  for  example.  Random 
file  index  keys  may  be  declared  in  blank  COMMON  if  desired,  since 
SPATCH  does  not  use  unlabelled  COMMON  blocks.  With  the  VAX  version  of 
SPATCH,  UOPEN  must  be  provided  to  open  the  file  containing  surface 
patch  data. 

UOPEN  is  called  at  the  beginning  of  execution,  before  any  data  is 
requested  from  UPATCH.  When  the  file  to  be  read  contains  other  data 
preceding  the  surface  patch  data,  UOPEN  can  be  used  to  position  the 
file  properly. 

The  third  user  routine,  UCLOSE,  is  provided  to  permit  closing  of 
the  user  data  file  following  execution.  Normally,  UCLOSE  will  not  be 
required  for  either  machine  version  (VAX,  CDC),  but  is  provided  for 
use  in  exceptional  circumstances.  UCLOSE  is  called  as  the  last  step 


APPENDIX  F 


MATERIAL  PROPERTIES  LIBRARY 


The  material  properties  library  used  by  PREP  (Section  4)  and  REFMT 


(Subsection  5 

.  1)  contains  materials  data  for  twelve  classes  of  materia 

ACRYL 

— 

Acrylics 

ALUMI 

- 

Aluminum  alloys 

CASTI 

- 

Cast  irons 

COPPR 

- 

Copper-based  alloys 

GLASS 

- 

Glasses 

MAGNS 

- 

Magnesium  Alloys 

NICKL 

- 

Nickel  alloys 

PLYMR 

- 

Polymeric  and  silicone  materials 

POLYC 

- 

Polycarbonates 

STEEL 

- 

Carbon  steels 

STSTL 

- 

Stainless  steels 

TITNM 

- 

Titanium 

Each  material 

in  the  library  is  described  by  an  accession  number  (five 

digits),  a  five- character  class  designation  (e.g.,  ALUMI),  a  text 
description,  and  a  list  of  numerical  property  values.  Stored  property 
information  consists  of: 

-  elastic  modulus? 

-  Poisson' s  ratio; 


elastic  shear  modulus; 
mass  density; 

coefficient  of  thermal  expansion;  and 
initial  yield  stress. 


At  present,  the  library  does  not  contain  material  data  for  the  plastic 
range  of  deformation  (e.g.,  strain-hardening  slope  or  full  stress  vs. 
strain  curves),  but  allowance  has  been  made  for  references  to  such 
tables,  stored  within  the  library.  Two  such  table  pointers  are  stored 
as  part  of  the  material  description  for  each  entry  in  the  library. 
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